Changeset 321 for trunk/LMDZ.MARS/libf/dyn3d

Timestamp:

Oct 21, 2011, 10:20:46 AM (13 years ago)

Author:

emillour

Message:

Mars GCM:

• Corrected small bug in newstart: initracer was not always used and thus some tracer indexes (igm_co2, igcm_h2o_vap,...) were not set. This however means that we now also call inifis from newstart and that we read in flags set in 'callphys.def' (required for sanity checks in initracer). Also adapted 'inichim_readcallphys': removed some obsolescent tests on number of tracers for given combinations of options.

EM

File:

• trunk/LMDZ.MARS/libf/dyn3d/newstart.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/dyn3d/newstart.F

@@ -347 +347 @@
       CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
 
+! also initialize various physics flags/settings which might be needed
+!    (for instance initracer needs to know about some flags, and/or
+!      'datafile' path may be changed by user)
+      call inifis(ngridmx,llm,day_ini,daysec,dtphys,
+     &                latfi,lonfi,airefi,rad,g,r,cpp)
+
 c=======================================================================
 c   lecture topographie, albedo, inertie thermique, relief sous-maille
@@ -361 +367 @@
         relief="mola"
 c     enddo
-
-! before using datareadnc, "datafile" must be set (normaly done in inifis)
-      datafile="/u/forget/WWW/datagcm/datafile" ! default value
-      call getin("datadir",datafile) ! in case user specified another path
 
       CALL datareadnc(relief,phis,alb,surfith,z0S,
@@ -457 +459 @@
         tnom(1:nqmx)=noms(1:nqmx)
         write(*,*) 'Newstart: updated tracer names'
+      else
+       ! initialize tracer names and indexes (igcm_co2, igcm_h2o_vap, ...)
+        call initracer(qsurf,co2ice)
       endif
 
@@ -883 +888 @@
           DO l=1,llm
             DO j=1,jjp1
-              DO i=1,iip1
+              DO i=1,iip1-1
                 q(i,j,l,igcm_h2o_vap)=150.e-6 * (rlatu(j)+pi/2.) / pi
               ENDDO
+              ! We want to have the very same value at lon -180 and lon 180
+              q(iip1,j,l,igcm_h2o_vap) = q(1,j,l,igcm_h2o_vap)
             ENDDO
           ENDDO