Custom Query (964 matches)

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.

Added two new scripts ws_remove.sh and ws_remove.bash to remove trailing white spaces and trailing tabs. It also report all files that contain tabs (excluding the makefile) so can be removed by hand.

These are a nice scripts to clean source code. Many thanks to David Robertson for coding these scripts.

During initialization, the mask at PSi_points is recomputed in metrics.F. The range for this is:

      DO j=Jstr,Jend
        DO i=Istr,Iend
          IF ((rmask...   THEN
            pmask(i,j)=1.0_r8
          ...
        END DO
      END DO

But the pmask has a range from Istr:Iend+1, and Jstr:Jend+1. So I think this range should be Istr,Iend+1 and Jstr,Jend+1.