Changeset 378 for codes/icosagcm/trunk/src

Timestamp:

05/25/16 09:57:29 (8 years ago)

Author:

ymipsl

Message:

Add vertical axis definition in output field (merged from heat/aquaplanet).

YM

File:

• codes/icosagcm/trunk/src/write_field.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/write_field.f90

@@ -1170 +1170 @@
       INTEGER :: nij
       CHARACTER(LEN=255) :: name
-      
+      INTEGER :: l,level_size, levId, dimlevId
             
       name=TRIM(ADJUSTL(name_in))
@@ -1222 +1222 @@
         status = NF90_DEF_DIM(ncid,'cell_i',ncell,ncellId)
         status = NF90_DEF_DIM(ncid,'nvert_i',nvert,nvertid)
-
+        level_size=0
         IF (Field(ind_b)%ndim==2)  THEN
           FieldVarId(NbField)%size=1
@@ -1230 +1230 @@
           ALLOCATE(FieldVarId(NbField)%nc_id(1))
           status = NF90_DEF_DIM(ncid,'lev',size(field(ind_b)%rval3d,2),dim3id)
+          level_size=size(field(ind_b)%rval3d,2)
         ELSE IF (Field(1)%ndim==4) THEN
           FieldVarId(NbField)%size=size(field(ind_b)%rval4d,3)
           ALLOCATE(FieldVarId(NbField)%nc_id(FieldVarId(NbField)%size))
           status = NF90_DEF_DIM(ncid,'lev',size(field(ind_b)%rval4d,2),dim3id)
+          level_size=size(field(ind_b)%rval4d,2)
 !          status = NF90_DEF_DIM(ncid,'Q',size(field(ind_b)%rval4d,3),dim4id)
         ENDIF
-      
+        PRINT*,"LEVEL_SIZE=",level_size
+
         status = NF90_DEF_DIM(ncid,'time_counter',NF90_UNLIMITED,timeId)
-      
+        IF (level_size>0) THEN
+          status = NF90_DEF_VAR(ncid,'lev',NF90_DOUBLE,(/ dim3id /),levId)
+          status = NF90_PUT_ATT(ncid,levId,"axis","Z")
+        ENDIF
+        
         status = NF90_DEF_VAR(ncid,'lon_i',NF90_DOUBLE,(/ ncellId /),lonId)
         status = NF90_PUT_ATT(ncid,lonId,"long_name","longitude")
@@ -1279 +1286 @@
       
         status = NF90_ENDDEF(ncid)      
+
+        if (level_size>0) status = NF90_PUT_VAR(ncid,levId,(/ (l,l=1,level_size) /))
          
          ALLOCATE(lon(ncell),lat(ncell),bounds_lon(0:nvert-1,ncell),bounds_lat(0:nvert-1,ncell))