Changeset 1948 for LMDZ5/trunk/libf/phy1d/1DUTILS.h

Timestamp:

Jan 27, 2014, 11:29:29 AM (11 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/1DUTILS.h (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