Changeset 632 for trunk/LMDZ.MARS/util/concatnc.F90

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

File:

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

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