Changeset 1073 for trunk/LMDZ.MARS/util/localtime.F90

Timestamp:

Oct 17, 2013, 4:19:26 PM (11 years ago)

Author:

tnavarro

Message:

Better handling of the first date of the file : Read and write the controle field, if possible, for zrecast, localtime and concatnc + cosmetic change in lslin. This is retrocompatible with previous versions.

File:

• trunk/LMDZ.MARS/util/localtime.F90 (modified) (11 diffs)

trunk/LMDZ.MARS/util/localtime.F90

@@ -35 +35 @@
 integer :: varid
 ! varid: [netcdf] variable ID #
-real, dimension(:), allocatable:: lat,lon,alt,time
+real, dimension(:), allocatable:: lat,lon,alt,ctl,time
 ! lat(): array, stores latitude coordinates
 ! lon(): array, stores longitude coordinates
 ! alt(): array, stores altitude coordinates
+! ctl(): array, stores controle array
 ! time(): array, stores time coordinates
 integer :: nbvar,nbvarfile,ndim
@@ -45 +46 @@
 ! nbvarfile: total # of variables in an input file
 ! ndim: [netcdf] # (3 or 4) of dimensions (for variables)
-integer :: latdim,londim,altdim,timedim
+integer :: latdim,londim,altdim,ctldim,timedim
 ! latdim: [netcdf] "latitude" dim ID
 ! londim: [netcdf] "longitude" dim ID
 ! altdim: [netcdf] "altdim" dim ID
+! ctldim: [netcdf] "controle" dim ID
 ! timedim: [netcdf] "timedim" dim ID
-integer :: latvar,lonvar,altvar,timevar
+integer :: latvar,lonvar,altvar,ctlvar,timevar
 ! latvar: [netcdf] ID of "latitude" variable
 ! lonvar: [netcdf] ID of "longitude" variable
 ! altvar: [netcdf] ID of "altitude" variable
+! ctlvar: [netcdf] ID of "controle" variable
 ! timevar: [netcdf] ID of "Time" variable
-integer :: latlen,lonlen,altlen,timelen,timelen_lt,timelen_tot
+integer :: latlen,lonlen,altlen,ctllen,timelen,timelen_lt,timelen_tot
 integer :: ilat,ilon,ialt,it
 ! latlen: # of elements of lat() array
 ! lonlen: # of elements of lon() array
 ! altvar: # of elements of alt() array
+! ctlvar: # of elements of ctl() array
 ! timelen: # of elemnets of time() array
 ! timelen_tot:# =timelen or timelen+1 (if 1 more time to interpolate needed)
@@ -246 +250 @@
 !  write(*,*) "altlen: ",altlen
 
+   ierr=NF_INQ_DIMID(nid,"index",ctldim)
+   ierr=NF_INQ_VARID(nid,"controle",ctlvar)
+   if (ierr.NE.NF_NOERR) then
+      write(*,*) 'Field <controle> is missing in file'//file
+      ctllen=0
+      !stop ""
+   else
+      ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
+   endif
+!  write(*,*) "controle: ",controle
+
 !==============================================================================
 ! 2.3. Read (and check compatibility of) dimensions of
@@ -255 +270 @@
       allocate(lon(lonlen))
       allocate(alt(altlen))
+      allocate(ctl(ctllen))
 #ifdef NC_DOUBLE
       ierr = NF_GET_VAR_DOUBLE(nid,latvar,lat)
       ierr = NF_GET_VAR_DOUBLE(nid,lonvar,lon)
       ierr = NF_GET_VAR_DOUBLE(nid,altvar,alt)
+      if (ctllen .ne. -1) ierr = NF_GET_VAR_DOUBLE(nid,ctlvar,ctl)
 #else
       ierr = NF_GET_VAR_REAL(nid,latvar,lat)
       ierr = NF_GET_VAR_REAL(nid,lonvar,lon)
       ierr = NF_GET_VAR_REAL(nid,altvar,alt)
+      if (ctllen .ne. -1) ierr = NF_GET_VAR_REAL(nid,ctlvar,ctl)
 #endif
 !==============================================================================
@@ -351 +369 @@
 
    ! Initialize output file's lat,lon,alt and time dimensions
-      call initiate (filename,lat,lon,alt,nout,&
+      call initiate (filename,lat,lon,alt,ctl,nout,&
            latdimout,londimout,altdimout,timedimout,timevarout)
    ! Initialize output file's aps,bps and phisinit variables
@@ -561 +579 @@
 
 !******************************************************************************
-Subroutine initiate (filename,lat,lon,alt,&
+Subroutine initiate (filename,lat,lon,alt,ctl,&
                      nout,latdimout,londimout,altdimout,timedimout,timevarout)
 !==============================================================================
@@ -586 +604 @@
 real, dimension(:), intent(in):: alt
 ! alt(): altitude
+real, dimension(:), intent(in):: ctl
+! ctl(): controle
 integer, intent(out):: nout
 ! nout: [netcdf] file ID
@@ -604 +624 @@
 !integer :: latdim,londim,altdim,timedim
 integer :: nvarid,ierr
+integer :: ctldimout
 ! nvarid: [netcdf] ID of a variable
 ! ierr: [netcdf]  return error code (from called subroutines)
@@ -625 +646 @@
 ierr = NF_DEF_DIM(nout, "longitude", size(lon), londimout)
 ierr = NF_DEF_DIM(nout, "altitude", size(alt), altdimout)
+if (size(ctl).ne.0) ierr = NF_DEF_DIM(nout, "index", size(ctl), ctldimout)
 ierr = NF_DEF_DIM(nout, "Time", NF_UNLIMITED, timedimout)
 
@@ -664 +686 @@
 
 !==============================================================================
-! 4. Write "altitude" (data and attributes)
+! 5. Write "altitude" (data and attributes)
 !==============================================================================
 
@@ -688 +710 @@
 ierr = NF_PUT_VAR_REAL (nout,nvarid,alt)
 #endif 
+
+!==============================================================================
+! 6. Write "controle" (data and attributes)
+!==============================================================================
+
+if (size(ctl).ne.0) then
+   ! Switch to netcdf define mode
+   ierr = NF_REDEF (nout)
+
+   #ifdef NC_DOUBLE
+   ierr = NF_DEF_VAR (nout,"controle",NF_DOUBLE,1,ctldimout,nvarid)
+   #else
+   ierr = NF_DEF_VAR (nout,"controle",NF_FLOAT,1,ctldimout,nvarid)
+   #endif
+
+   ierr = NF_PUT_ATT_TEXT (nout,nvarid,"title",18,"Control parameters")
+
+   ! End netcdf define mode
+   ierr = NF_ENDDEF(nout)
+
+   #ifdef NC_DOUBLE
+   ierr = NF_PUT_VAR_DOUBLE (nout,nvarid,ctl)
+   #else
+   ierr = NF_PUT_VAR_REAL (nout,nvarid,ctl)
+   #endif
+endif
 
 end Subroutine initiate