Changeset 860 for trunk/LMDZ.MARS/util

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.

Location:

trunk/LMDZ.MARS/util

Files:

• hrecast.F90 (modified) (6 diffs)
• zrecast.F90 (modified) (6 diffs)

trunk/LMDZ.MARS/util/hrecast.F90

@@ -14 +14 @@
 !     from 90 to -90 (again, you may omit the endpoints).
 ! EM 09/2009
+! TN 01/2013 : Adapted for large output files with at least 2 variables > 2 GiB
 
 implicit none
@@ -91 +92 @@
 logical :: out_has_poles ! Flag: true if output grid includes North and South
                          ! poles 
+                         
+integer, dimension(4) :: edges,corner ! needed to write variables for big files
+
 !===============================================================================
 ! 1. Input parameters
@@ -537 +541 @@
 
 ! 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
   stop "Error: Could not define latitude dimension"
@@ -858 +862 @@
 
 ! 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
   stop "Error: Could not write Time data to output file"
@@ -944 +948 @@
     enddo
     ! write interpolated data to output
-    ierr=NF_PUT_VAR_REAL(outfid,var_id(i),out_3d_data)
+    corner(:)=1
+    edges(1)=lonlength
+    edges(2)=latlength
+    edges(3)=timelength
+    ierr=NF_PUT_VARA_REAL(outfid,var_id(i),corner(1:3),edges(1:3),out_3d_data)
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error: Could not write ",trim(var(i))," data to output file"
@@ -995 +1003 @@
     enddo
     ! write interpolated data to output
-    ierr=NF_PUT_VAR_REAL(outfid,var_id(i),out_4d_data)
+    corner(:)=1
+    edges(1)=lonlength
+    edges(2)=latlength
+    edges(3)=altlength
+    edges(4)=timelength
+    ierr=NF_PUT_VARA_REAL(outfid,var_id(i),corner,edges,out_4d_data)
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Error: Could not write ",trim(var(i))," data to output file"

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"