Changeset 887 for trunk/LMDZ.VENUS/libf

Timestamp:

Feb 18, 2013, 2:46:14 PM (12 years ago)

Author:

slebonnois

Message:

SL: add rcm1d tool in Venus and Titan physics to run the code in 1D, and associated small modifications

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• physiq.F (modified) (1 diff)
• profile.F (modified) (3 diffs)
• rcm1d.F (added)

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

@@ -1302 +1302 @@
 #endif
 
-c ---------- Sorties testphys1d -------------
-      
-      if (klon.eq.1) then
-        call writeg1d(klon,klev,t_seri,"Temp","Temperature")
-        call writeg1d(klon,1,ftsol,"tsurf","Surface Temp")
-      DO k = 1, klev
-      DO i = 1, klon
-c        tmpout(i,k) = real(heat(i,k))/RDAY
-         tmpout(i,k) = heat(i,k)/RDAY
-      ENDDO
-      ENDDO
-        call writeg1d(klon,klev,tmpout,
-     .                 "heat","Solar heating")
-      DO k = 1, klev
-      DO i = 1, klon
-c        tmpout(i,k) = real(dtrad(i,k))
-         tmpout(i,k) = dtrad(i,k)
-      ENDDO
-      ENDDO
-        call writeg1d(klon,klev,tmpout,
-     .                 "dtrad","IR cooling")
-        call writeg1d(klon,klev,lwnet,"lwnet","Net LW flux")
-        call writeg1d(klon,klev,swnet,"swnet","Net SW flux")
-        call writeg1d(klon,klev,fluxt,"flux_vdf","Turbulent flux")
-        call writeg1d(klon,klev,flux_ajs,"flux_ajs","Dry adjust. flux")
-        call writeg1d(klon,klev,flux_ec,"flux_ec","Ec flux")
-c       call writeg1d(klon,1,solsw,"surfsw","Net SW flux at surface")
-c       call writeg1d(klon,1,sollw,"surflw","Net LW flux at surface")
-        call writeg1d(klon,1,radsol,"surfnet","Net flux at surface")
-        call writeg1d(klon,klev,d_t_vdf,"dt_vdf","DT from clmain")
-        call writeg1d(klon,klev,d_t_ajs,"dt_ajs","DT from ajsec")
-        call writeg1d(klon,klev,d_t_ec,"dt_ec","DT from Ec")
-      endif
-      
 c====================================================================
 c Si c'est la fin, il faut conserver l'etat de redemarrage

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

@@ -2 +2 @@
       IMPLICIT NONE
 c=======================================================================
-c     Subroutine utilisee dans le modele 1-D  "testphys1d"
+c     Subroutine utilisee dans le modele 1-D  "rcm1d"
 c     pour l'initialisation du profil atmospherique
 c=======================================================================
+c
+c   VERSION VENUS
 c
 c   differents profils d'atmospheres. T=f(z)
@@ -56 +58 @@
 c   ----------------
 
-c la lecture se fait dans le testphys1d.def, ouvert par testphys1d.F
+c la lecture se fait dans le rcm1d.def, ouvert par rcm1d.F
       READ(unit,*) 
       READ(unit,*) 
@@ -66 +68 @@
       READ(unit,*) largeur
       READ(unit,*) hauteur
-      CLOSE(unit)
 
 c-----------------------------------------------------------------------