Changeset 927 for codes/icosagcm/devel/src/output
- Timestamp:
- 06/21/19 00:55:59 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/output/write_field.f90
r883 r927 292 292 FieldIndex(NbField)=1 293 293 294 PRINT *, 'Creating header for writefield : ', TRIM(FieldName(NbField)), ndomain, ind_b, ind_e ! FIXME295 PRINT *, mesh_glo%ndomain, SIZE(mesh_glo%primal_own)294 ! PRINT *, 'Creating header for writefield : ', TRIM(FieldName(NbField)), ndomain, ind_b, ind_e ! FIXME 295 ! PRINT *, mesh_glo%ndomain, SIZE(mesh_glo%primal_own) 296 296 SELECT CASE(Field(ind_b)%field_type) 297 297 CASE(field_T) 298 PRINT *, ' field_type == field_T' ! FIXME298 ! PRINT *, ' field_type == field_T' ! FIXME 299 299 IF(single) THEN ! include halos 300 300 cells => mesh_glo%primal_all … … 303 303 END IF 304 304 CASE(field_Z) 305 PRINT *, ' field_type == field_Z' ! FIXME305 ! PRINT *, ' field_type == field_Z' ! FIXME 306 306 IF(single) THEN ! include halos 307 307 cells => mesh_glo%dual_all … … 311 311 END SELECT 312 312 313 PRINT *, 'writefield : ind_b, ind_e :', ind_b, ind_e313 ! PRINT *, 'writefield : ind_b, ind_e :', ind_b, ind_e 314 314 315 315 ncell=0 … … 318 318 ncell = ncell + cells(ind)%ncell 319 319 END DO 320 PRINT *, 'writefield : found ',ncell, ' cells.' ! FIXME320 ! PRINT *, 'writefield : found ',ncell, ' cells.' ! FIXME 321 321 322 322 status = NF90_CREATE(TRIM(ADJUSTL(name))//'.nc', NF90_CLOBBER, ncid) … … 416 416 DO ind=ind_b,ind_e 417 417 n_end = n_begin + cells(ind)%ncell 418 PRINT *, 'create_header_gen ', n_begin, n_end, SHAPE(cells(ind)%bnds_lon), SHAPE(bounds_lon)418 ! PRINT *, 'create_header_gen ', n_begin, n_end, SHAPE(cells(ind)%bnds_lon), SHAPE(bounds_lon) 419 419 lat(n_begin+1:n_end) = cells(ind)%lat(:) *180./Pi 420 420 lon(n_begin+1:n_end) = cells(ind)%lon(:) *180./Pi
Note: See TracChangeset
for help on using the changeset viewer.