Changeset 1948 for LMDZ5/trunk/libf/phy1d/lmdz1d.F

Timestamp:

Jan 27, 2014, 11:29:29 AM (10 years ago)

Author:

idelkadi

Message:

• Option pour tourner sans thermiques ni ajustement sec si ifalg_thermals<0 (phylmd/physiq.F90)

• Modifications permettant d'imposer des profils initiaux de traceurs

lus dans un tracer.inp.001 (phylmd/1DUTILS.h, 1D_read_forc_cases.h et lmdz1d.F)

File:

• LMDZ5/trunk/libf/phy1d/lmdz1d.F (modified) (3 diffs)

LMDZ5/trunk/libf/phy1d/lmdz1d.F

@@ -90 +90 @@
 
         integer :: kmax = llm
-        integer nlev_max,llm700
-        parameter (nlev_max = 1000)
+        integer llm700,nq1,nq2
+        INTEGER, PARAMETER :: nlev_max=1000, nqmx=1000
         real timestep, frac, timeit
         real height(nlev_max),tttprof(nlev_max),qtprof(nlev_max),
@@ -97 +97 @@
      .              ugprof(nlev_max),vgprof(nlev_max),wfls(nlev_max),
      .              dqtdxls(nlev_max),dqtdyls(nlev_max),
-     .              dqtdtls(nlev_max),thlpcar(nlev_max)
+     .              dqtdtls(nlev_max),thlpcar(nlev_max),
+     .              qprof(nlev_max,nqmx)
 
         real    :: fff
@@ -427 +428 @@
       call infotrac_init 
 
+      if (nqtot>nqmx) STOP'Augmenter nqmx dans lmdz1d.F'
       allocate(q(llm,nqtot)) ; q(:,:)=0.
       allocate(dq(llm,nqtot))