Changeset 1953

Timestamp:

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

Author:

fhourdin

Message:

Correction in for readprofiles_fire

Location:

LMDZ5/trunk/libf/phy1d

Files:

• 1DUTILS.h (modified) (1 diff)
• 1D_read_forc_cases.h (modified) (4 diffs)

LMDZ5/trunk/libf/phy1d/1DUTILS.h

@@ -3017 +3017 @@
 
        open(ilesfile,file='trac.inp.001',status='old',iostat=ierr)
-       nt1=ntrac
-       nt2=ntrac
         if (ierr /= 0) then
             print*,'WARNING : trac.inp does not exist'

LMDZ5/trunk/libf/phy1d/1D_read_forc_cases.h

@@ -3 +3 @@
 ! forcing_radconv = .T. : Pure radiative-convective equilibrium:
 !----------------------------------------------------------------------
+
+
+      nq1=0
+      nq2=0
 
       if (forcing_les .or. forcing_radconv
@@ -62 +66 @@
         ug(l)   = ugprof(kmax)-frac*( ugprof(kmax)- ugprof(kmax-1))
         vg(l)   = vgprof(kmax)-frac*( vgprof(kmax)- vgprof(kmax-1))
-        q(l,nq1:nq2)=qprof(kmax,nq1:nq2)
+        IF (nq2>0) q(l,nq1:nq2)=qprof(kmax,nq1:nq2)
      s               -frac*(qprof(kmax,nq1:nq2)-qprof(kmax-1,nq1:nq2))
         omega(l)=   wfls(kmax)-frac*(   wfls(kmax)-   wfls(kmax-1))
@@ -87 +91 @@
             ug(l)   = ugprof(k)-frac*( ugprof(k)- ugprof(k-1))
             vg(l)   = vgprof(k)-frac*( vgprof(k)- vgprof(k-1))
-            q(l,nq1:nq2)=qprof(k,nq1:nq2)
+            IF (nq2>0) q(l,nq1:nq2)=qprof(k,nq1:nq2)
      s                   -frac*(qprof(k,nq1:nq2)-qprof(k-1,nq1:nq2))
             omega(l)=   wfls(k)-frac*(   wfls(k)-   wfls(k-1))
@@ -108 +112 @@
             vg(l)   = vgprof(1)
             omega(l)=   wfls(1)
-            q(l,nq1:nq2)=qprof(1,nq1:nq2)
+            IF (nq2>0) q(l,nq1:nq2)=qprof(1,nq1:nq2)
             dq_dyn(l,1)  =dqtdtls(1)
             dt_cooling(l)=thlpcar(1)