Changeset 1675 for trunk/LMDZ.VENUS/libf

Timestamp:

Mar 16, 2017, 9:51:39 AM (8 years ago)

Author:

slebonnois

Message:

SL: update of xml file for xios output + a few modifications for run with uppr atmosphere

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• flott_gwd_ran.F90 (modified) (1 diff)
• hgardfou.F (modified) (3 diffs)
• load_ksi.F (modified) (2 diffs)
• moldiff_red.F90 (modified) (4 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90

@@ -138 +138 @@
 !(similar to observations)
     KMIN = 1.E-5        ! Min horizontal wavenumber
+!   KMIN = 6.3E-6       ! Min horizontal wavenumber
     KMAX = 1.E-4        ! Max horizontal wavenumber
     !Online output: one value only

trunk/LMDZ.VENUS/libf/phyvenus/hgardfou.F

@@ -22 +22 @@
       DATA firstcall /.TRUE./
       IF (firstcall) THEN
-         PRINT*, 'hgardfou garantit la temperature dans [20,1200] K'
+         PRINT*, 'hgardfou garantit la temperature dans [15,1200] K'
          firstcall = .FALSE.
       ENDIF
@@ -29 +29 @@
       DO k = 1, klev
          DO i = 1, klon
+!!!! MODIF GG to avoid crash after 78--> 95 extension!!
+!!!  WARNING: it has to be review/removed when the extension to the
+!!!  thermosphere is completed (physical processes and ionosphere added)
+          IF (k.LT.85) THEN
             zt(i) = t(i,k)
+          ENDIF
          ENDDO
 #ifdef CRAY
@@ -49 +54 @@
          ENDIF
 #ifdef CRAY
-         CALL WHENFLT(klon, zt, 1, 20.0, jadrs, jbad)
+         CALL WHENFLT(klon, zt, 1, 15.0, jadrs, jbad)
 #else
          jbad = 0
          DO i = 1, klon
-         IF (zt(i).LT.20.0) THEN
+         IF (zt(i).LT.15.0) THEN
             jbad = jbad + 1
             jadrs(jbad) = i

trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F

@@ -45 +45 @@
       nlve = klev
 
+cc      GG modif below 
+c----------------------------------
+c   Initialisation of values to 0 
+c     (for all vertical levels)
+c----------------------------------
+
+      ksive(0:klev+1,0:klev+1,nnuve,nbmat) = 0.0
+
 c ------------------------
 c Loading the ksi file
@@ -77 +85 @@
         read(10,*)
         read(10,*) m,Nb
-        if (m.ne.nlve) then
+cc      GG modif below  
+        if (nlve.le.78.and.m.ne.nlve) then
          write(*,*) 'This is subroutine load_ksi'
          print*,'Dimension problem between ksi.txt and nlve'

trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90

@@ -1462 +1462 @@
 
         compteur=0
-        do while (pnew2 .ge. pp(il))    
+        do while (Pnew2 .ge. pp(il))    
         compteur=compteur+1
         do nn=1,nq
@@ -1474 +1474 @@
         if (compteur .ge. 100000) then
         print*,'error moldiff_red infinite loop'
-        print*,ig,il,pp(il),tt(nl),pnew2,qnew(il,:),Znew2
+        print*,ig,il,pp(il),tt(nl),Pnew2,qnew(il,:),Znew2
         stop
         endif
@@ -1564 +1564 @@
 
         compteur=0
-        do while (pnew2 .ge. pp(il))    
+        do while (Pnew2 .ge. pp(il))    
         compteur=compteur+1
         do nn=1,nq
@@ -1576 +1576 @@
         if (compteur .ge. 100000) then
         print*,'pb moldiff_red infinite loop'
-        print*,ig,nl,T(nl),pnew2,qnew(il,:),Znew2
+        print*,ig,nl,T(nl),Pnew2,qnew(il,:),Znew2
         stop
         endif