Changeset 927 for trunk/LMDZ.VENUS/libf

Timestamp:

Apr 3, 2013, 2:20:46 PM (12 years ago)

Author:

slebonnois

Message:

SL: mise a jour de quelques docs + correction bug dans newstart (Venus)

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• comg1d.h (deleted)
• newstart.F (modified) (4 diffs)
• writeg1d.F (deleted)

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

@@ -238 +238 @@
       ip1jmp1old = (imold+1-1/imold)*(jmold+1-1/jmold)
       
+c dans run.def
+      CALL getin('planet_type',planet_type)
+      call ini_cpdet
+
 c=======================================================================
 c   CHANGEMENT DE CONSTANTES CONTENUES DANS tab_cntrl
@@ -1024 +1028 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       write (*,*) 'T ', T(1,jjp1,1)  ! INFO
+! pour info:
+! Si extension verticale, la T est extrapolee constante au-dessus de lmold
 
 c passage grille physique pour restartphy.nc 
@@ -1055 +1061 @@
           do i=1,iip1
             ucov( i,j,l ) = unat( i,j,l ) * cu(i,j)
-c           ucov( i,j,l ) = 0
+! pour info:
+! Si extension verticale, on impose u=0 au-dessus de lmold
+            if (l.gt.lmold) ucov( i,j,l ) = 0
           end do
         end do
@@ -1065 +1073 @@
           do i=1,iim
             vcov( i,j,l ) = vnat( i,j,l ) * cv(i,j)
-c           vcov( i,j,l ) = 0
+! pour info:
+! Si extension verticale, on impose v=0 au-dessus de lmold
+            if (l.gt.lmold) vcov( i,j,l ) = 0
           end do
           vcov( iip1,j,l ) = vcov( 1,j,l )