[Dart-dev] DART/branches Revision: 11527

[dart at ucar.edu]

nancy at ucar.edu
2017-04-26 16:31:00 -0600 (Wed, 26 Apr 2017)
89
take a stab at adding the grid info to the output
files.  this has not been tested yet.




Modified: DART/branches/rma_output_files/models/ROMS/model_mod.f90
===================================================================
--- DART/branches/rma_output_files/models/ROMS/model_mod.f90	2017-04-26 22:22:28 UTC (rev 11526)
+++ DART/branches/rma_output_files/models/ROMS/model_mod.f90	2017-04-26 22:31:00 UTC (rev 11527)
@@ -615,26 +615,6 @@
                                    'nc_write_model_atts', 'inquire '//trim(filename))
 call nc_check(nf90_Redef(ncFileID),'nc_write_model_atts',   'redef '//trim(filename))
 
-! We need the dimension ID for the number of copies/ensemble members, and
-! we might as well check to make sure that Time is the Unlimited dimension.
-! Our job is create the 'model size' dimension.
-
-call nc_check(nf90_inq_dimid(ncid=ncFileID, name='copy', dimid=MemberDimID), &
-                           'nc_write_model_atts', 'copy dimid '//trim(filename))
-call nc_check(nf90_inq_dimid(ncid=ncFileID, name='time', dimid=  TimeDimID), &
-                           'nc_write_model_atts', 'time dimid '//trim(filename))
-
-if ( TimeDimID /= unlimitedDimId ) then
-   write(string1,*)'Time Dimension ID ',TimeDimID, &
-             ' should equal Unlimited Dimension ID',unlimitedDimID
-   call error_handler(E_ERR,'nc_write_model_atts:', string1, source, revision, revdate)
-endif
-
-! Define the model size / state variable dimension / whatever ...
-
-call nc_check(nf90_def_dim(ncid=ncFileID, name='StateVariable', len=model_size, &
-        dimid = StateVarDimID),'nc_write_model_atts', 'state def_dim '//trim(filename))
-
 ! Write Global Attributes
 
 call DATE_AND_TIME(crdate,crtime,crzone,values)
@@ -657,21 +637,6 @@
 ! to know the geometry, the output variables (PS,U,V,T,Q,...) etc. We're skipping
 ! complicated part.
 
-if ( output_state_vector ) then
-
-   ! Create a variable for the state vector
-   ! Define the actual (3D) state vector, which gets filled as time goes on ...
-
-   call nc_check(nf90_def_var(ncid=ncFileID, name='state', xtype=nf90_real, &
-                 dimids=(/StateVarDimID,MemberDimID,unlimitedDimID/),VarID=StateVarID),&
-                 'nc_write_model_atts','state def_var '//trim(filename))
-   call nc_check(nf90_put_att(ncFileID,StateVarID,'long_name','model state or fcopy'),&
-                 'nc_write_model_atts', 'state long_name '//trim(filename))
-
-   call nc_check(nf90_enddef(ncFileID),'nc_write_model_atts','state enddef '//trim(filename))
-
-else
-
    ! We need to output the prognostic variables.
    ! Define the new dimensions IDs
 
@@ -842,7 +807,6 @@
    call nc_check(nf90_put_var(ncFileID, VarID, WDEP ), &
                 'nc_write_model_atts', 'z_w put_var '//trim(filename))
 
-   endif
 
 ! Flush the buffer and leave netCDF file open
 call nc_check(nf90_sync(ncFileID), 'nc_write_model_atts', 'atts sync')
@@ -905,25 +869,6 @@
 
 write(filename,*) 'ncFileID', ncFileID
 
-! make sure ncFileID refers to an open netCDF file,
-
-call nc_check(nf90_inq_dimid(ncFileID, 'copy', dimid=CopyDimID), &
-            'nc_write_model_vars', 'inq_dimid copy '//trim(filename))
-
-call nc_check(nf90_inq_dimid(ncFileID, 'time', dimid=TimeDimID), &
-            'nc_write_model_vars', 'inq_dimid time '//trim(filename))
-
-if ( output_state_vector ) then
-
-   call nc_check(NF90_inq_varid(ncFileID, 'state', VarID), &
-                 'nc_write_model_vars', 'state inq_varid '//trim(filename))
-   call nc_check(NF90_put_var(ncFileID,VarID,state_vec,start=(/1,copyindex,timeindex/)),&
-                 'nc_write_model_vars', 'state put_var '//trim(filename))
-
-endif
-
-! Flush the buffer and leave netCDF file open
-
 call nc_check(nf90_sync(ncFileID), 'nc_write_model_vars', 'sync '//trim(filename))
 
 ierr = 0 ! If we got here, things went well.