Changeset 1073 for trunk/LMDZ.MARS/util/zrecast.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/zrecast.F90 (modified) (8 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -60 +60 @@
 !              constants (radius, R, etc.) are now read from file
 ! TN 01/2013 : Adapted for large output files with at least 2 variables > 2 GiB
+! TN 10/2013 : Read and write controle field, if available
 !
 implicit none
@@ -80 +81 @@
 integer nbvar4dinfile ! # of 4D (lon,lat,alt,time) variables in input file
 integer outfid ! NetCDF output file ID
-integer lon_dimid,lat_dimid,alt_dimid,time_dimid ! NetCDF dimension IDs
-integer lon_varid,lat_varid,alt_varid,time_varid
+integer lon_dimid,lat_dimid,alt_dimid,time_dimid,ctl_dimid ! NetCDF dimension IDs
+integer lon_varid,lat_varid,alt_varid,time_varid,ctl_varid
 integer gcm_layers_dimid ! NetCDF dimension ID for # of layers in GCM
 integer sigma_varid,aps_varid,bps_varid
@@ -103 +104 @@
 integer latlength ! # of grid points along latitude
 integer altlength ! # of grid point along altitude (of input datasets)
+real,dimension(:),allocatable :: ctl ! controle
+integer ctllength ! # of grid points along controle
 real,dimension(:),allocatable :: time ! time
 integer timelength ! # of points along time
@@ -439 +442 @@
   if (ierr.ne.NF_NOERR) then
       stop "Error: Failed to get altitude length"
+  endif
+endif
+
+! controle
+ierr=NF_INQ_DIMID(infid,"index",tmpdimid)
+if (ierr.ne.NF_NOERR) then
+  write(*,*) "Could not get controle dimension ID"
+  ctllength = 0
+else
+  ierr=NF_INQ_VARID(infid,"controle",tmpvarid)
+  if (ierr.ne.NF_NOERR) then
+    stop "Error: Failed to get controle ID"
+  else
+    ierr=NF_INQ_DIMLEN(infid,tmpdimid,ctllength)
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to get controle length"
+    else
+      allocate(ctl(ctllength),stat=ierr)
+      if (ierr.ne.0) then
+        write(*,*) "Error: Failed to allocate ctl(ctllength)"
+        write(*,*) "     ctllength=",ctllength
+        stop
+      endif
+      ierr=NF_GET_VAR_REAL(infid,tmpvarid,ctl)
+      if (ierr.ne.NF_NOERR) then
+        stop "Error: Failed reading controle"
+      endif
+    endif
   endif
 endif
@@ -1030 +1061 @@
 endif
 
+! controle
+if (ctllength .ne. 0) then
+  ierr=NF_DEF_DIM(outfid,"index",ctllength,ctl_dimid)
+  if (ierr.ne.NF_NOERR) then
+    write(*,*) "Error: Could not define controle dimension"
+    write(*,*) NF_STRERROR(ierr)
+    stop
+  endif
+endif
+
 ! GCM layers (for sigma or aps and bps)
 ierr=NF_DEF_DIM(outfid,"GCM_layers",altlength,gcm_layers_dimid)
@@ -1162 +1203 @@
   endif
 endif ! of if (have_sigma)
+
+! controle
+if (ctllength .ne. 0) then
+  ierr=NF_DEF_VAR(outfid,"controle",NF_REAL,1,ctl_dimid,ctl_varid)
+  if (ierr.ne.NF_NOERR) then
+    write(*,*) "Error: Could not define controle variable"
+    write(*,*) NF_STRERROR(ierr)
+    stop
+  endif
+
+  ! controle attributes
+  text='Control parameters'
+  ierr=NF_PUT_ATT_TEXT(outfid,ctl_varid,'title',len_trim(text),text)
+  if (ierr.ne.NF_NOERR) then
+    stop "Error: Problem writing title for controle"
+  endif
+endif
 
 ! GCM_layers
@@ -1509 +1567 @@
 endif
 
+! Write controle
+if (ctllength .ne. 0) then
+  ierr=NF_PUT_VAR_REAL(outfid,ctl_varid,ctl)
+  if (ierr.ne.NF_NOERR) then
+    write(*,*) "Error: Could not write controle data to output file"
+    write(*,*) NF_STRERROR(ierr)
+    stop
+  endif
+endif
+
 ! Write time
 ierr=NF_PUT_VARA_REAL(outfid,time_varid,1,timelength,time)
@@ -1771 +1839 @@
 subroutine init_planet_const(infid)
 ! initialize planetary constants using the "controle" array in the file
-! if "cointrole" array not found in file, look for it in "diagfi.nc"
+! if "controle" array not found in file, look for it in "diagfi.nc"
 use planet_const
 use netcdf