Changeset 1120 for trunk/LMDZ.TITAN/libf

Timestamp:

Dec 3, 2013, 3:05:40 PM (11 years ago)

Author:

slebonnois

Message:

SL: Titan and Venus modifications following a modif in dyn3d[par].

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys.F90 (modified) (2 diffs)
• physiq.F (modified) (1 diff)
• rcm1d.F (modified) (16 diffs)

trunk/LMDZ.TITAN/libf/phytitan/clesphys.h

@@ -10 +10 @@
        LOGICAL ok_orodr,ok_orolf,ok_gw_nonoro
        INTEGER nbapp_rad, nbapp_chim, iflag_con, iflag_ajs
-       REAL    ecritphy
+       REAL    ecriphy
        INTEGER lev_histmth, lev_histday
        REAL    solaire
@@ -39 +39 @@
 
        COMMON/clesphys_r/                                               &
-     &     ecritphy, solaire, z0, lmixmin, ksta, inertie, emis,         &
+     &     ecriphy, solaire, z0, lmixmin, ksta, inertie, emis,          &
      &     tx,tcorrect,p_prodaer,xnuf,xvis,xir
 

trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90

@@ -139 +139 @@
 !Config         en jours.
 !
-       ecritphy = 1.
-       CALL getin('ecritphy', ecritphy)
+       ecriphy = 1.
+       CALL getin('ecritphy', ecriphy)
 !
 !Config Key  = if_ebil
@@ -479 +479 @@
   write(numout,*)' Sortie journaliere = ', ok_journe
   write(numout,*)' Sortie instantanee = ', ok_instan
-  write(numout,*)' frequence sorties = ', ecritphy  
+  write(numout,*)' frequence sorties = ', ecriphy  
   write(numout,*)' Sortie bilan d energie, if_ebil =', if_ebil
   write(numout,*)' Duree de l annee = ',year_day

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -591 +591 @@
          ENDIF
 
-         ecrit_ins = NINT(RDAY/dtime*ecritphy)  ! Fraction de jour reglable
+         ecrit_ins = NINT(RDAY/dtime*ecriphy)  ! Fraction de jour reglable
          IF (ok_instan) THEN
          WRITE(lunout,*)'La frequence de sortie instant. est de ', 

trunk/LMDZ.TITAN/libf/phytitan/rcm1d.F

@@ -4 +4 @@
       use control_mod
       use comgeomphy
+      USE phys_state_var_mod
       use cpdet_mod, only: ini_cpdet
       IMPLICIT NONE
@@ -22 +23 @@
 c   
 c=======================================================================
-
-c Version TITAN a tester et verifier  
-c  - verifier pour Ls... 
-c  - faire un profile.F ...
 
 #include "dimensions.h"
@@ -52 +49 @@
       REAL play(llm)   ! Pressure at the middle of the layers (Pa)
       REAL plev(llm+1) ! intermediate pressure levels (pa)
-      REAL psurf,tsurf(1)      
+      REAL psurf      
       REAL u(llm),v(llm)  ! zonal, meridional wind
       REAL gru,grv   ! prescribed "geostrophic" background wind
       REAL temp(llm)   ! temperature at the middle of the layers
       REAL,allocatable :: q(:,:) ! tracer mixing ratio (e.g. kg/kg)
-      REAL tsoil(nsoilmx)   ! subsurface soik temperature (K)
       REAL zlay(llm)   ! altitude estimee dans les couches (km)
       REAL long(1),lati(1),area(1)
       REAL cufi(1),cvfi(1)
-      REAL phisfi(1),albedo(1)
-      REAL solsw(1),sollwdown(1),dlw(1),radsol(1)
-      REAL zmea(1), zstd(1)
-      REAL zsig(1), zgam(1), zthe(1)
-      REAL zpic(1), zval(1)
+      REAL phisfi(1)
 
 c    Physical and dynamical tandencies (e.g.  m.s-2, K/s, Pa/s)
@@ -85 +77 @@
       COMMON/cpdetvenus/cppdyn,nu_venus,t0_venus
       REAL cppdyn,nu_venus,t0_venus
-      real pi
 
 c=======================================================================
@@ -168 +159 @@
 c     print*,nbapp_rad
 c LU DANS PHYSIQ.DEF...
-      nbapp_rad = 10.
+      nbapp_rad = 100.
 
       PRINT *,'nombre de jours simules ?'
@@ -176 +167 @@
       ndt=ndt*day_step     
       dtphys=daysec/day_step  
-      dtime=dtphys
 
 c Pression de surface sur la planete
@@ -196 +186 @@
       long(1)=0.E+0
 
-c  Initialisation albedo 
-c  ----------------------
-c ne sert pas ici...
-      albedo(1)=0.3
-c      PRINT *,'Albedo du sol nu ?'
-c      READ(unit,*) albedo(1)
-c      PRINT *,albedo(1)
-
 c   Initialisation speciales "physiq"
 c   ---------------------------------
@@ -209 +191 @@
       CALL init_phys_lmdz(iim,jjm,llm,1,(/1/))
       call initcomgeomphy
-      call infotrac_init
       call ini_cpdet
 
@@ -271 +252 @@
       pks=cpp*(psurf/preff)**rcp
 
+c  init des variables pour phyredem
+c  --------------------------------
+      call phys_state_var_init
+
 c  profil de temperature et altitude au premier appel
 c  --------------------------------------------------
@@ -289 +274 @@
       print*,"               Pression        Altitude     Temperature"
       ilayer=1
-      tsurf(1)=tmp3(0)
+      ftsol(1)=tmp3(0)
        temp(1)=tmp3(1)
        zlay(1)=tmp3(1)*tmp1(1)
-      print*,"           0",tsurf(1)
+      print*,"           0",ftsol(1)
       print*,ilayer,play(ilayer),zlay(ilayer),temp(ilayer)
       DO ilayer=2,nlayer
@@ -303 +288 @@
 c     ~~~~~~~~~~~~~~~~~~~~~~~
       DO isoil=1,nsoil
-         tsoil(isoil)=93.
+         ftsoil(1,isoil)=ftsol(1)
       ENDDO
 
@@ -331 +316 @@
          print*,'Ls=',zls*180./pi
 
+c  Initialisation albedo 
+c  ----------------------
+
+      falbe(1)=0.3
+
 c  Ecriture de "startphy.nc"
 c  -------------------------
@@ -338 +328 @@
 
       solsw(1)    = 0.
-      sollwdown(1)= 0.
-      dlw(1)      = 0.
+      sollw(1)    = 0.
+      fder(1)     = 0.
       radsol(1)   = 0.
      
@@ -345 +335 @@
       soil_model  = .true.
 
-      call phyredem("startphy.nc  ",
-     .              lati,long,
-     .              tsurf,tsoil,albedo,
-     .              solsw,sollwdown,dlw,radsol,
-     .    zmea, zstd, zsig, zgam, zthe, zpic, zval, 
-     .              temp)
+      call phyredem("startphy.nc")
+
+c  deallocation des variables phyredem
+c  -----------------------------------
+      call phys_state_var_end
 
 c=======================================================================
@@ -463 +452 @@
 
       OPEN(11,file='profile.new')
-      write (11,*) tsurf
       DO ilayer=1,nlayer
         write (11,*) zlay(ilayer),temp(ilayer),tmp1(ilayer)
@@ -476 +464 @@
 #include "../dyn3d/disvert_noterre.F"
 #include "../dyn3d/abort_gcm.F"
-!#include "../dyn3d/dump2d.F"
-
-c***********************************************************************
-      function ssum(n,sx,incx)
-c
-      IMPLICIT NONE
-c
-      integer n,incx,i,ix
-      real ssum,sx((n-1)*incx+1)
-c
-      ssum=0.
-      ix=1
-      do 10 i=1,n
-         ssum=ssum+sx(ix)
-         ix=ix+incx
-10    continue
-c
-      return
-      end