Changeset 860 for trunk/LMDZ.MARS/util/zrecast.F90

Timestamp:

Jan 9, 2013, 5:49:56 PM (12 years ago)

Author:

tnavarro

Message:

Modification to create large files, that is to say with at least 2 variables > 2GiB. In that case, it is necessary to have an unlimited dimension.

File:

• trunk/LMDZ.MARS/util/zrecast.F90 (modified) (6 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -59 +59 @@
 ! EM 11/2012 : Adapted so it can be used on "generic" model outputs; planet
 !              constants (radius, R, etc.) are now read from file
+! TN 01/2013 : Adapted for large output files with at least 2 variables > 2 GiB
+!
 implicit none
 
@@ -138 +140 @@
 logical :: auto_mcd_levels ! Flag: specific case for MCD automatic above
                            ! local surface levels
+integer,dimension(4) :: edges,corner ! needed to write variables for big files
+
 !===============================================================================
 ! 1. Input parameters
@@ -1019 +1023 @@
 
 ! time
-ierr=NF_DEF_DIM(outfid,"Time",timelength,time_dimid)
+ierr=NF_DEF_DIM(outfid,"Time",NF_UNLIMITED,time_dimid)
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Error: Could not define Time dimension"
@@ -1506 +1510 @@
 
 ! Write time
-ierr=NF_PUT_VAR_REAL(outfid,time_varid,time)
+ierr=NF_PUT_VARA_REAL(outfid,time_varid,1,timelength,time)
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Error: Could not write Time data to output file"
@@ -1518 +1522 @@
 
 ! Write surface pressure
-
-ierr=NF_PUT_VAR_REAL(outfid,ps_varid,ps)
+corner(:)=1
+edges(1)=lonlength
+edges(2)=latlength
+edges(3)=timelength
+ierr=NF_PUT_VARA_REAL(outfid,ps_varid,corner(1:3),edges(1:3),ps)
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Error: Could not write ps data to output file"
@@ -1581 +1588 @@
                       miss_val,ps,press,za_gcm,plevel,outdata)
   ! write result to output file
-  ierr=NF_PUT_VAR_REAL(outfid,za_varid,outdata)
+  corner(:)=1
+  edges(1)=lonlength
+  edges(2)=latlength
+  edges(3)=nblev
+  edges(4)=timelength
+  ierr=NF_PUT_VARA_REAL(outfid,var_id(i),corner,edges,outdata)
   if (ierr.ne.NF_NOERR) then
     write(*,*) "Error, Failed to write zareoid to output file"