Changeset 1948

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)

Location:

LMDZ5/trunk/libf

Files:

• phy1d/1DUTILS.h (modified) (3 diffs)
• phy1d/1D_read_forc_cases.h (modified) (4 diffs)
• phy1d/lmdz1d.F (modified) (3 diffs)
• phylmd/physiq.F90 (modified) (3 diffs)

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

@@ -2968 +2968 @@
 
 !=====================================================================
-      subroutine readprofiles(nlev_max,kmax,height,
+      subroutine readprofiles(nlev_max,kmax,ntrac,height,
      .           thlprof,qtprof,uprof,
      .           vprof,e12prof,ugprof,vgprof,
      .           wfls,dqtdxls,dqtdyls,dqtdtls,
-     .           thlpcar)
+     .           thlpcar,tracer,nt1,nt2)
       implicit none
 
-        integer nlev_max,kmax,kmax2
+        integer nlev_max,kmax,kmax2,ntrac
         logical :: llesread = .true.
 
@@ -2982 +2982 @@
      .       ugprof(nlev_max),vgprof(nlev_max),wfls(nlev_max),
      .       dqtdxls(nlev_max),dqtdyls(nlev_max),dqtdtls(nlev_max),
-     .           thlpcar(nlev_max)
+     .           thlpcar(nlev_max),tracer(nlev_max,ntrac)
 
         integer, parameter :: ilesfile=1
-        integer :: ierr,irad,imax,jtot,k
+        integer :: ierr,irad,imax,jtot,k,itrac,nt1,nt2
         logical :: lmoist,lcoriol,ltimedep
         real :: xsize,ysize
@@ -3015 +3015 @@
         end do
         close(ilesfile)
+
+       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'
+        else
+        read (ilesfile,*) kmax2,nt1,nt2
+        if (nt2>ntrac) then
+          stop'Augmenter le nombre de traceurs dans traceur.def'
+        endif
+        if (kmax .ne. kmax2) then
+          print *, 'fichiers prof.inp et lscale.inp incompatibles :'
+          print *, 'nbre de niveaux : ',kmax,' et ',kmax2
+          stop 'lecture profiles'
+        endif
+        do k=1,kmax
+          read (ilesfile,*) height(k),(tracer(k,itrac),itrac=nt1,nt2)
+        end do
+        close(ilesfile)
+        endif
 
         return

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

@@ -24 +24 @@
 !----------------------------------------------------------------------
 
-      call readprofiles(nlev_max,kmax,height,
+      call readprofiles(nlev_max,kmax,nqtot,height,
      .           tttprof,qtprof,uprof,vprof,
      .           e12prof,ugprof,vgprof,
      .           wfls,dqtdxls,dqtdyls,dqtdtls,
-     .           thlpcar) 
+     .           thlpcar,qprof,nq1,nq2)
       endif
 
@@ -62 +62 @@
         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)
+     s               -frac*(qprof(kmax,nq1:nq2)-qprof(kmax-1,nq1:nq2))
         omega(l)=   wfls(kmax)-frac*(   wfls(kmax)-   wfls(kmax-1))
 
@@ -85 +87 @@
             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)
+     s                   -frac*(qprof(k,nq1:nq2)-qprof(k-1,nq1:nq2))
             omega(l)=   wfls(k)-frac*(   wfls(k)-   wfls(k-1))
             dq_dyn(l,1)=dqtdtls(k)-frac*(dqtdtls(k)-dqtdtls(k-1))
@@ -104 +108 @@
             vg(l)   = vgprof(1)
             omega(l)=   wfls(1)
+            q(l,nq1:nq2)=qprof(1,nq1:nq2)
             dq_dyn(l,1)  =dqtdtls(1)
             dt_cooling(l)=thlpcar(1)

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)) 

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -2296 +2296 @@
        'AVANT LA CONVECTION SECHE , iflag_thermals=' &
        ,iflag_thermals,'   nsplit_thermals=',nsplit_thermals
-  if(iflag_thermals.lt.0) then
+  if(iflag_thermals<=0) then
      !  Rien
      !  ====
-     IF(prt_level>9)WRITE(lunout,*)'pas de convection'
+     IF(prt_level>9)WRITE(lunout,*)'pas de modele de panache thermique'
+     IF(prt_level>9.and.iflag_thermals<0)WRITE(lunout,*)'pas d ajustement convectif'
 
 
@@ -2446 +2447 @@
      ! Dans le cas contraire, on demarre au niveau 1.
 
-     if (iflag_thermals.ge.13.or.iflag_thermals.eq.0) then
+     if (iflag_thermals>=13.or.iflag_thermals<=0) then
 
         if(iflag_thermals.eq.0) then
@@ -2461 +2462 @@
         ! non nulles numeriquement pour des mailles non concernees.
 
-        if (iflag_thermals.eq.0) then
+        if (iflag_thermals==0) then
            CALL ajsec_convV2(paprs, pplay, t_seri,q_seri &
                 , d_t_ajsb, d_q_ajsb)
-        else
+        else if (iflag_thermals>0) then
            CALL ajsec(paprs, pplay, t_seri,q_seri,limbas &
                 , d_t_ajsb, d_q_ajsb)