Changeset 1363 for LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F

Timestamp:

Apr 16, 2010, 11:50:10 AM (14 years ago)

Author:

Ehouarn Millour

Message:

Added possibility to run in "Shallow Water" mode. SW mode is automatically set if the number of atmospheric layers is 1.
To use SW mode, the model should be compiled without physics (makelmdz_fcm -p nophys) and/or without calls to the physics (i.e. set iflag_phys=0 in gcm.def).

-Updated some definitions and comments in run.def & gcm.def

-Fixed misplaced #ifdef CPP_EARTH in calfis.F + some write(lunout,*)

-Specific initialisation of fields for SW are done in sw_case_williamson91_6 (called by iniacademic, when read_start=.false.)

• Checked (on Vargas & Brodie) that these changes don't alter usual bench GCM outputs.

EM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F

@@ -1354 +1354 @@
                 call Gather_Field(ucov,ip1jmp1,llm,0)
                 call Gather_Field(teta,ip1jmp1,llm,0)
+                call Gather_Field(pk,ip1jmp1,llm,0)
                 call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
@@ -1559 +1560 @@
                 call Gather_Field(ucov,ip1jmp1,llm,0)
                 call Gather_Field(teta,ip1jmp1,llm,0)
+                call Gather_Field(pk,ip1jmp1,llm,0)
                 call Gather_Field(phi,ip1jmp1,llm,0)
                 do iq=1,nqtot
@@ -1624 +1626 @@
                 call Gather_Field(phis,ip1jmp1,1,0)
                 if (mpi_rank==0) then
-                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
+                 CALL writehist(itau,vcov,ucov,teta,phi,q,masse,ps,phis)
                 endif
 !                CALL writehist_p(histid, histvid, itau,vcov ,