Changeset 632 for trunk/LMDZ.MARS/util

Timestamp:

Apr 25, 2012, 10:23:52 AM (13 years ago)

Author:

emillour

Message:

Mars GCM utilities: minor improvements:

• zrecast also looks for a 'diagfi1.nc' file when looking for 'phisinit' variable
• concatnc now also includes 'aires' in its output file

EM

Location:

trunk/LMDZ.MARS/util

Files:

• concatnc.F90 (modified) (6 diffs)
• zrecast.F90 (modified) (1 diff)

trunk/LMDZ.MARS/util/concatnc.F90

@@ -713 +713 @@
 !==============================================================================
 ! Purpose:
-! Copy ap() , bp(), aps(), bps() and phisinit() from input file to outpout file
+! Copy ap() , bp(), aps(), bps(), aire() and phisinit()
+! from input file to outpout file
 !==============================================================================
 ! Remarks:
@@ -743 +744 @@
 real,dimension(:),allocatable :: ap,bp ! hybrid vertical coordinates
 real,dimension(:),allocatable :: sigma ! sigma levels
+real,dimension(:,:),allocatable :: aire ! mesh areas
 real,dimension(:,:),allocatable :: phisinit ! Ground geopotential
 integer :: apsid,bpsid,sigmaid,phisinitid
@@ -748 +750 @@
 integer :: ierr
 integer :: tmpvarid ! temporary variable ID
+logical :: area ! is "aire" available ?
 logical :: phis ! is "phisinit" available ?
 logical :: hybrid ! are "aps" and "bps" available ?
@@ -843 +846 @@
 endif ! of if (.not.hybrid)
 
+! mesh area
+allocate(aire(lonlen,latlen),stat=ierr)
+if (ierr.ne.0) then
+  write(*,*) "init2: failed to allocate aire!"
+  stop
+endif
+ierr=NF_INQ_VARID(infid,"aire",tmpvarid)
+if (ierr.ne.NF_NOERR) then
+  write(*,*)"init2 warning: Failed to get aire ID."
+  area = .false.
+else
+  ierr=NF_GET_VAR_REAL(infid,tmpvarid,aire)
+  if (ierr.ne.NF_NOERR) then
+    stop "init2 Error: Failed reading aire"
+  endif
+  area = .true.
+endif
+
 ! ground geopotential phisinit
 allocate(phisinit(lonlen,latlen),stat=ierr)
@@ -958 +979 @@
 
 !==============================================================================
-! 2.2. phisinit()
-!==============================================================================
+! 2.2. aire() and phisinit()
+!==============================================================================
+
+if (area) then
+  ! define aire
+  call def_var(nout,"aire","Mesh area","m2",2,&
+           (/londimout,latdimout/),tmpvarid,ierr)
+  if (ierr.ne.NF_NOERR) then
+     stop "init2 Error: Failed to def_var aire"
+  endif
+  
+  ! write aire
+#ifdef NC_DOUBLE
+  ierr=NF_PUT_VAR_DOUBLE(outfid,tmpvarid,aire)
+#else
+  ierr=NF_PUT_VAR_REAL(outfid,tmpvarid,aire)
+#endif
+  if (ierr.ne.NF_NOERR) then
+    stop "init2 Error: Failed to write aire"
+  endif
+endif ! of if (area)
 
 IF (phis) THEN
@@ -990 +1030 @@
 if (allocated(sigma)) deallocate(sigma)
 if (allocated(phisinit)) deallocate(phisinit)
+if (allocated(aire)) deallocate(aire)
 
 end subroutine init2

trunk/LMDZ.MARS/util/zrecast.F90

@@ -499 +499 @@
   if (ierr.ne.NF_NOERR) then
     write(*,*) "Problem: Could not find/open that file"
-    infile2="phisinit.nc"
+    infile2="diagfi1.nc"
     write(*,*) "         Trying file ",trim(infile2)
     ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
     if (ierr.ne.NF_NOERR) then
-      write(*,*) "Error: Could not open that file either"
-      stop "Might as well stop here"
+      write(*,*) "Problem: Could not find/open that file"
+      infile2="phisinit.nc"
+      write(*,*) "         Trying file ",trim(infile2)
+      ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
+      if (ierr.ne.NF_NOERR) then
+        write(*,*) "Error: Could not open that file either"
+        stop "Might as well stop here"
+      endif
     endif
   endif