Changeset 378 for codes/icosagcm
- Timestamp:
- 05/25/16 09:57:29 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/write_field.f90
r351 r378 1170 1170 INTEGER :: nij 1171 1171 CHARACTER(LEN=255) :: name 1172 1172 INTEGER :: l,level_size, levId, dimlevId 1173 1173 1174 1174 name=TRIM(ADJUSTL(name_in)) … … 1222 1222 status = NF90_DEF_DIM(ncid,'cell_i',ncell,ncellId) 1223 1223 status = NF90_DEF_DIM(ncid,'nvert_i',nvert,nvertid) 1224 1224 level_size=0 1225 1225 IF (Field(ind_b)%ndim==2) THEN 1226 1226 FieldVarId(NbField)%size=1 … … 1230 1230 ALLOCATE(FieldVarId(NbField)%nc_id(1)) 1231 1231 status = NF90_DEF_DIM(ncid,'lev',size(field(ind_b)%rval3d,2),dim3id) 1232 level_size=size(field(ind_b)%rval3d,2) 1232 1233 ELSE IF (Field(1)%ndim==4) THEN 1233 1234 FieldVarId(NbField)%size=size(field(ind_b)%rval4d,3) 1234 1235 ALLOCATE(FieldVarId(NbField)%nc_id(FieldVarId(NbField)%size)) 1235 1236 status = NF90_DEF_DIM(ncid,'lev',size(field(ind_b)%rval4d,2),dim3id) 1237 level_size=size(field(ind_b)%rval4d,2) 1236 1238 ! status = NF90_DEF_DIM(ncid,'Q',size(field(ind_b)%rval4d,3),dim4id) 1237 1239 ENDIF 1238 1240 PRINT*,"LEVEL_SIZE=",level_size 1241 1239 1242 status = NF90_DEF_DIM(ncid,'time_counter',NF90_UNLIMITED,timeId) 1240 1243 IF (level_size>0) THEN 1244 status = NF90_DEF_VAR(ncid,'lev',NF90_DOUBLE,(/ dim3id /),levId) 1245 status = NF90_PUT_ATT(ncid,levId,"axis","Z") 1246 ENDIF 1247 1241 1248 status = NF90_DEF_VAR(ncid,'lon_i',NF90_DOUBLE,(/ ncellId /),lonId) 1242 1249 status = NF90_PUT_ATT(ncid,lonId,"long_name","longitude") … … 1279 1286 1280 1287 status = NF90_ENDDEF(ncid) 1288 1289 if (level_size>0) status = NF90_PUT_VAR(ncid,levId,(/ (l,l=1,level_size) /)) 1281 1290 1282 1291 ALLOCATE(lon(ncell),lat(ncell),bounds_lon(0:nvert-1,ncell),bounds_lat(0:nvert-1,ncell))
Note: See TracChangeset
for help on using the changeset viewer.