Changeset 397 for trunk

Timestamp:

Nov 18, 2011, 12:49:09 PM (13 years ago)

Author:

emillour

Message:

Mars GCM utility concatnc : minor upgrade. Handle cases when ap and bp are not available but aps and bps are.
FF

File:

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

trunk/LMDZ.MARS/util/concatnc.F90

@@ -16 +16 @@
 !   of data (which is the case for outputs from 'zrecast')
 !   (EM, April 2010)
+! + handle absence of ap() and bp() if aps and bps are available 
+!    (case of stats file) FF, November 2011
 ! ********************************************************
 
@@ -748 +750 @@
 logical :: phis ! is "phisinit" available ?
 logical :: hybrid ! are "aps" and "bps" available ?
+logical :: apbp ! are "ap" and "bp" available ?
 
 !==============================================================================
@@ -797 +800 @@
   if (ierr.ne.NF_NOERR) then
     write(*,*) "Ooops. Failed to get ap ID. OK."
-    hybrid=.false.
+    apbp=.false.
   else
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,ap)
-    hybrid=.true.
+    apbp=.true.
     if (ierr.ne.NF_NOERR) then
       stop "Error: Failed reading ap"
@@ -816 +819 @@
   if (ierr.ne.NF_NOERR) then
     write(*,*) "Ooops. Failed to get bp ID. OK."
-    hybrid=.false.
+    apbp=.false.
   else
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bp)
-    hybrid=.true.
+    apbp=.true.
     if (ierr.ne.NF_NOERR) then
       stop "Error: Failed reading bp"
@@ -900 +903 @@
   endif
 
-! define ap
-  call def_var(nout,"ap","hybrid sigma at interlayers"," ",1,&
+  if (apbp) then
+!   define ap
+    call def_var(nout,"ap","hybrid sigma at interlayers"," ",1,&
              (/interlayerdimout/),apid,ierr)
-  if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var ap"
-  endif
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to def_var ap"
+    endif
 
 ! write ap
 #ifdef NC_DOUBLE
-  ierr=NF_PUT_VAR_DOUBLE(outfid,apid,ap)
-#else
-  ierr=NF_PUT_VAR_REAL(outfid,apid,ap)
-#endif
-  if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write ap"
-  endif
+    ierr=NF_PUT_VAR_DOUBLE(outfid,apid,ap)
+#else
+    ierr=NF_PUT_VAR_REAL(outfid,apid,ap)
+#endif
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to write ap"
+    endif
 
 ! define bp
-  call def_var(nout,"bp","hybrid sigma at interlayers"," ",1,&
+    call def_var(nout,"bp","hybrid sigma at interlayers"," ",1,&
              (/interlayerdimout/),bpid,ierr)
-  if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var bp"
-  endif
-
-! write bp
-#ifdef NC_DOUBLE
-  ierr=NF_PUT_VAR_DOUBLE(outfid,bpid,bp)
-#else
-  ierr=NF_PUT_VAR_REAL(outfid,bpid,bp)
-#endif
-  if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write bp"
-  endif
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to def_var bp"
+    endif
+
+!   write bp
+#ifdef NC_DOUBLE
+    ierr=NF_PUT_VAR_DOUBLE(outfid,bpid,bp)
+#else
+    ierr=NF_PUT_VAR_REAL(outfid,bpid,bp)
+#endif
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to write bp"
+    endif
+  endif ! of if (apbp)
 
 else