Opened 11 years ago

Closed 11 years ago

#580 closed upgrade (Done)

Metadata for all tracers climatology

Reported by: arango Owned by: arango
Priority: major Milestone: Release ROMS/TOMS 3.6
Component: Nonlinear Version: 3.6
Keywords: Cc:

Description

The climatology tracer (active and passive) metadata is now computed internally in mod_ncparam.F:

      DO i=1,MT
        varid=varid+1
        IF (varid.gt.MV) THEN
          WRITE (stdout,60) MV, varid
          STOP
        END IF
        idTclm(i)=varid
        DO ng=1,Ngrids
          Fscale(varid,ng)=1.0_r8
          Iinfo(1,varid,ng)=r3dvar
        END DO
        WRITE (Vname(1,varid),'(a)')                                    &
     &        TRIM(ADJUSTL(Vname(1,idTvar(i))))
        WRITE (Vname(2,varid),'(a,a)')                                  &
     &        TRIM(ADJUSTL(Vname(2,idTvar(i)))), ' climatology'
        WRITE (Vname(3,varid),'(a)')                                    &
     &        TRIM(ADJUSTL(Vname(3,idTvar(i))))
        WRITE (Vname(4,varid),'(a,a)')                                  &
     &        TRIM(Vname(1,varid)), ', scalar, series'
        WRITE (Vname(5,varid),'(a,a)')                                  &
     &        TRIM(ADJUSTL(Vname(1,idTvar(i)))), '_time'
      END DO

Notice that:

  • NetCDF variable name, Vname(1,:), is the same as the basic ROMS state tracer metadata.
  • Units, Vname(3,:), are the same as the basic ROMS state tracer data.
  • Default time attribute, Vname(5,:), is the same as the variable name but with the _time suffix. For example, temp_time, salt_time, and so on. Recall that other time-variables names are allowed provided that the input NetCDF variable has the "time" attribute with the appropriate value.
  • The following lines are removed from metadata file varinfo.dat because they are not longer needed:
    'temp'                                             ! Input
      'potential temperature climatology'
      'Celsius'                                        ! [Celsius]
      'temp, scalar, series'
      'temp_time'
      'idTclm(itemp)'
      'r3dvar'
      1.0d0
    
    'salt'                                             ! Input
      'salinity climatology'
      'nondimensional'                                 ! [PSU]
      'salt, scalar, series'
      'salt_time'
      'idTclm(isalt)'
      'r3dvar'
      1.0d0
    

I forgot to update the metadata when I updated the tracer climatology logic in src:ticket:569. Many thanks to Jiang-peng Liang for bringing this to my attention.

I corrected the format number in various ecosystem *_inp.h files when reporting the LtracerCLM switches. Thanks again to Jiang-peng Liang for reporting this typo.

Change History (1)

comment:1 by arango, 11 years ago

Resolution: Done
Status: newclosed
Note: See TracTickets for help on using tickets.