Changeset 436 for codes/icosagcm/trunk/src/observable.f90

Timestamp:

06/16/16 06:29:34 (8 years ago)

Author:

dubos

Message:

Correct output for Lagrangian vertical coordinate

File:

• codes/icosagcm/trunk/src/observable.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/observable.f90

@@ -101 +101 @@
     ELSE
        CALL output_field("temp",f_buf_i)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,85000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,85000.)
        CALL output_field("t850",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,50000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,50000.)
        CALL output_field("t500",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,preff)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,preff)
        CALL output_field("SST",f_buf_s)       
     END IF
@@ -135 +135 @@
        !       CALL output_field("pv",f_qv)
        
-       CALL vertical_interp(f_ps,f_buf_ulon,f_buf_s,85000.)
+       CALL vertical_interp(f_pmid,f_buf_ulon,f_buf_s,85000.)
        CALL output_field("u850",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_ulon,f_buf_s,50000.)
+       CALL vertical_interp(f_pmid,f_buf_ulon,f_buf_s,50000.)
        CALL output_field("u500",f_buf_s)
        
-       CALL vertical_interp(f_ps,f_buf_ulat,f_buf_s,85000.)
+       CALL vertical_interp(f_pmid,f_buf_ulat,f_buf_s,85000.)
        CALL output_field("v850",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_ulat,f_buf_s,50000.)
+       CALL vertical_interp(f_pmid,f_buf_ulat,f_buf_s,50000.)
        CALL output_field("v500",f_buf_s)
 
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,85000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,85000.)
        CALL output_field("w850",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,50000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,50000.)
        CALL output_field("w500",f_buf_s)    
 
        CALL w_omega(f_ps, f_u, f_buf_i)
        CALL output_field("omega",f_buf_i)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,85000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,85000.)
        CALL output_field("omega850",f_buf_s)
-       CALL vertical_interp(f_ps,f_buf_i,f_buf_s,50000.)
+       CALL vertical_interp(f_pmid,f_buf_i,f_buf_s,50000.)
        CALL output_field("omega500",f_buf_s)
     END IF