Custom Query (986 matches)

Corrected a problem in get_2dfld.F and get_2dfldr.F when reading point data. The Fmin and Fmax where not assigned, so the correct range of the processed data was incorrectly reported to standard output as zero. However, the correct data was used in ROMS internally.

Many thanks to Diego Narvaez for reporting this problem in the ​forum.

I also corrected a typo in the makefile. Many thanks to Kate hedstrom for reporting this typo.

Updated the Model Coupling Toolkit (MCT) library to its latest version 2.6.0 which was ​released on April 30, 2009.

Updated files get_ngfld.F and get_ngfldr.F to allow reading 3D variables with Jrec=1. The variable Jrec is used to differentiate reading between 3D and 2D arrays:

            IF ((Jrec.gt.1).or.                                         &
     &          ((Jrec.eq.1).and.(Iinfo(7,ifield,ng).gt.0))) THEN
              npts=(Iend-Istr+1)*Jrec
              CALL netcdf_get_fvar (ng, model, ncfile, Vname(1,ifield), &
     &                              A,                                  &
     &                              ncid = ncid,                        &
     &                              start = (/1,1,Trec/),               &
     &                              total = (/Iend-Istr+1,Jrec,1/))
            ELSE
              npts=Iend-Istr+1
              CALL netcdf_get_fvar (ng, model, ncfile, Vname(1,ifield), &
     &                              A,                                  &
     &                              ncid = ncid,                        &
     &                              start = (/1,Trec/),                 &
     &                              total = (/Iend-Istr+1,1/))
            END IF

Although this is not needed currently, it will be used in the future to read other fields.

Many thanks to John Warner for bringing this to my attention.