Changeset 1286


Ignore:
Timestamp:
Dec 17, 2009, 2:47:10 PM (15 years ago)
Author:
Laurent Fairhead
Message:

Les parametres definissant la transition eau glace dans les nuages sont
maintenant lus dans physiq.def. Les valeurs par defaut donnent les memes
resultats que precedemment. JLD

Location:
LMDZ4/trunk/libf
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r1279 r1286  
    185185      PARAMETER (testita = 9)
    186186
    187       logical , parameter :: flag_verif = .true.
     187      logical , parameter :: flag_verif = .false.
    188188     
    189189
  • LMDZ4/trunk/libf/dyn3dpar/leapfrog_p.F

    r1279 r1286  
    172172      INTEGER testita
    173173      PARAMETER (testita = 9)
     174
     175      logical , parameter :: flag_verif = .false.
    174176     
    175177c declaration liees au parallelisme
     
    12861288
    12871289
    1288             IF( itau. EQ. itaufinp1 ) then 
     1290            IF( itau. EQ. itaufinp1 ) then
     1291
     1292              if (flag_verif) then
     1293                write(79,*) 'ucov',ucov
     1294                write(80,*) 'vcov',vcov
     1295                write(81,*) 'teta',teta
     1296                write(82,*) 'ps',ps
     1297                write(83,*) 'q',q
     1298                WRITE(85,*) 'q1 = ',q(:,:,1)
     1299                WRITE(86,*) 'q3 = ',q(:,:,3)
     1300              endif
     1301 
    12891302
    12901303c$OMP MASTER
  • LMDZ4/trunk/libf/phylmd/conf_phys.F90

    r1279 r1286  
    126126  INTEGER,SAVE :: iflag_pdf_omp
    127127  REAL,SAVE :: rad_froid_omp, rad_chau1_omp, rad_chau2_omp
     128  REAL,SAVE :: t_glace_min_omp, t_glace_max_omp
    128129  REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
    129130  REAL,SAVE :: qsol0_omp
     
    834835
    835836!
     837!Config Key  = t_glace_min
     838!Config Desc = 
     839!Config Def  = 258.
     840!Config Help =
     841!
     842  t_glace_min_omp = 258.
     843  call getin('t_glace_min',t_glace_min_omp)
     844
     845!
     846!Config Key  = t_glace_max
     847!Config Desc = 
     848!Config Def  = 273.13
     849!Config Help =
     850!
     851  t_glace_max_omp = 273.13
     852  call getin('t_glace_max',t_glace_max_omp)
     853
     854!
    836855!Config Key  = top_height
    837856!Config Desc =
     
    13741393    rad_chau1 = rad_chau1_omp
    13751394    rad_chau2 = rad_chau2_omp
     1395    t_glace_min = t_glace_min_omp
     1396    t_glace_max = t_glace_max_omp
    13761397    top_height = top_height_omp
    13771398    overlap = overlap_omp
     
    15441565  write(numout,*)' tau_ratqs = ',tau_ratqs
    15451566  write(numout,*)' top_height = ',top_height
     1567  write(numout,*)' rad_froid = ',rad_froid
     1568  write(numout,*)' rad_chau1 = ',rad_chau1
     1569  write(numout,*)' rad_chau2 = ',rad_chau2
     1570  write(numout,*)' t_glace_min = ',t_glace_min
     1571  write(numout,*)' t_glace_max = ',t_glace_max
    15461572  write(numout,*)' overlap = ',overlap
    15471573  write(numout,*)' cdmmax = ',cdmmax
  • LMDZ4/trunk/libf/phylmd/newmicro.F

    r1279 r1286  
    7373      REAL coef, coef_froi, coef_chau
    7474      PARAMETER (coef_chau=0.13, coef_froi=0.09)
    75       REAL seuil_neb, t_glace
    76       PARAMETER (seuil_neb=0.001, t_glace=273.0-15.0)
     75      REAL seuil_neb
     76      PARAMETER (seuil_neb=0.001)
    7777      INTEGER nexpo ! exponentiel pour glace/eau
    7878      PARAMETER (nexpo=6)
     
    143143         DO k = 1, klev
    144144            DO i = 1, klon
    145                zfice2(i,k) = 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
     145c               zfice2(i,k) = 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
     146               zfice2(i,k) = 1.0 - (t(i,k)-t_glace_min) /
     147     &                             (t_glace_max-t_glace_min)
    146148               zfice2(i,k) = MIN(MAX(zfice2(i,k),0.0),1.0)
    147149c     IM Total Liquid/Ice water content                                   
     
    339341                 
    340342                  zfice1 = MIN(
    341      &                 MAX( 1.0 - (t(i,k)-t_glace) / (273.13-t_glace)
    342      &                 ,0.0),1.0)**nexpo
     343     &                 MAX( 1.0 - (t(i,k)-t_glace_min) /
     344     &                    (t_glace_max-t_glace_min),0.0),1.0)**nexpo
    343345                 
    344346                  radius = rad_chaud * (1.-zfice1) + rad_froid * zfice1
  • LMDZ4/trunk/libf/phylmd/nuage.h

    r766 r1286  
    11!
    2 ! $Header$
     2! $Id$
    33!
    4       REAL rad_froid, rad_chau1, rad_chau2
     4      REAL rad_froid, rad_chau1, rad_chau2, t_glace_max, t_glace_min
    55
    6       common /nuagecom/ rad_froid,rad_chau1, rad_chau2
     6      common /nuagecom/ rad_froid,rad_chau1, rad_chau2,t_glace_max,     &
     7     &                  t_glace_min
    78!$OMP THREADPRIVATE(/nuagecom/)
Note: See TracChangeset for help on using the changeset viewer.