Changeset 28 for trunk/mesoscale


Ignore:
Timestamp:
Jan 25, 2011, 4:18:50 PM (14 years ago)
Author:
aslmd
Message:

LMD_MM_MARS: ancienne et nouvelle physique sur le mm coeur dynamique

Location:
trunk/mesoscale/LMD_MM_MARS
Files:
2 added
17 edited

Legend:

Unmodified
Added
Removed

  • trunk/mesoscale/LMD_MM_MARS/SIMU/runmeso

    r11 r28  
    6161#echo '***********************'
    6262echo '** Simulation directory  >>> '${simu_folder}
     63   ######## NEWPHYS
     64   ## if [[ ! ( -f "run.def" ) ]]
     65   ## then
     66   ##    echo 'new physics: you got to have run.def !'
     67   ##    exit
     68   ## fi
    6369if [[ ! ( -f "namelist.input" ) ]]
    6470then

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM

    r11 r28  
    184184state  real   hi_gc          igj   dyn_em  1  Z  i1  "HI"               "water ice"               "kg/kg"
    185185state  real   tsoil_gc       igj   dyn_em  1  Z  i1  "TSOIL"            "soil temperatures"       "K"
     186state  real   isoil_gc       igj   dyn_em  1  Z  i1  "ISOIL"            "soil thermal inertia"    "tiu"
     187state  real   dsoil_gc       igj   dyn_em  1  Z  i1  "DSOIL"            "soil depths"             "m"
    186188####
    187189#### variables dans wrf
    … …  
    194196state  real   mars_cice       ij   misc  1  -  i012rd  "MARS_CICE"   "co2 ice"                 "kg"
    195197state  real   mars_tsoil     ilj   misc  1  -  i012rd  "MARS_TSOIL"  "soil temperatures"       "K"
     198state  real   mars_isoil     ilj   misc  1  -  i012rd  "MARS_ISOIL"  "soil thermal inertia"    "tiu"
     199state  real   mars_dsoil     ilj   misc  1  -  i012rd  "MARS_DSOIL"  "soil depths"             "m"
    196200state  real   mars_gw        ilj   misc  1  -  i012rd  "MARS_GW"     "gwparam"                 " "
    197201####

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis.inc

    r11 r28  
    1414               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    1515               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    16                walbedodat,winertiedat,wphisfi,                         &
     16               walbedodat,,wphisfi,                         &
    1717               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    18                wtheta,wpsi)
     18               wtheta,wpsi,wvolcapa)
    1919

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis1.inc

    r11 r28  
    1414               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    1515               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    16                walbedodat,winertiedat,wphisfi,                         &
     16#ifdef NEWPHYS
     17               walbedodat, wphisfi, wvolcapa,                          &
     18#else
     19               walbedodat, winertiedat, wphisfi,                       &
     20#endif
    1721               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    18                wtheta,wpsi)
    19 
     22               wtheta, wpsi)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis2.inc

    r11 r28  
    1717               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    1818               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    19                walbedodat,winertiedat,wphisfi,                         &
     19#ifdef NEWPHYS
     20               walbedodat, wphisfi, wvolcapa,                          &
     21#else
     22               walbedodat, winertiedat, wphisfi,                       &
     23#endif
    2024               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    21                wtheta,wpsi)
     25               wtheta, wpsi)
     26
    2227
    2328       CASE(2)
    … …  
    3136               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    3237               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    33                walbedodat,winertiedat,wphisfi,                         &
     38#ifdef NEWPHYS
     39               walbedodat, wphisfi, wvolcapa,                          &
     40#else
     41               walbedodat, winertiedat, wphisfi,                       &
     42#endif
    3443               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    35                wtheta,wpsi)
     44               wtheta, wpsi)
     45
    3646!
    3747!       CASE(3:)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis3.inc

    r11 r28  
    1717               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    1818               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    19                walbedodat,winertiedat,wphisfi,                         &
     19#ifdef NEWPHYS
     20               walbedodat, wphisfi, wvolcapa,                          &
     21#else
     22               walbedodat, winertiedat, wphisfi,                       &
     23#endif
    2024               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    21                wtheta,wpsi)
     25               wtheta, wpsi)
    2226
    2327       CASE(2)
    … …  
    3135               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    3236               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    33                walbedodat,winertiedat,wphisfi,                         &
     37#ifdef NEWPHYS
     38               walbedodat, wphisfi, wvolcapa,                          &
     39#else
     40               walbedodat, winertiedat, wphisfi,                       &
     41#endif
    3442               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    35                wtheta,wpsi)
     43               wtheta, wpsi)
    3644
    3745       CASE(3)
    … …  
    4553               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    4654               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    47                walbedodat,winertiedat,wphisfi,                         &
     55#ifdef NEWPHYS
     56               walbedodat, wphisfi, wvolcapa,                          &
     57#else
     58               walbedodat, winertiedat, wphisfi,                       &
     59#endif
    4860               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    49                wtheta,wpsi)
     61               wtheta, wpsi)
    5062
    5163

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis4.inc

    r11 r28  
    1717               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    1818               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    19                walbedodat,winertiedat,wphisfi,                         &
     19#ifdef NEWPHYS
     20               walbedodat, wphisfi, wvolcapa,                          &
     21#else
     22               walbedodat, winertiedat, wphisfi,                       &
     23#endif
    2024               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    21                wtheta,wpsi)
     25               wtheta, wpsi)
    2226
    2327       CASE(2)
    … …  
    3135               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    3236               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    33                walbedodat,winertiedat,wphisfi,                         &
     37#ifdef NEWPHYS
     38               walbedodat, wphisfi, wvolcapa,                          &
     39#else
     40               walbedodat, winertiedat, wphisfi,                       &
     41#endif
    3442               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    35                wtheta,wpsi)
     43               wtheta, wpsi)
    3644
    3745       CASE(3)
    … …  
    4553               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    4654               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    47                walbedodat,winertiedat,wphisfi,                         &
     55#ifdef NEWPHYS
     56               walbedodat, wphisfi, wvolcapa,                          &
     57#else
     58               walbedodat, winertiedat, wphisfi,                       &
     59#endif
    4860               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    49                wtheta,wpsi)
     61               wtheta, wpsi)
    5062
    5163       CASE(4)
    … …  
    5971              wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
    6072              wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
    61               walbedodat,winertiedat,wphisfi,                         &
    62               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
    63               wtheta,wpsi)
     73#ifdef NEWPHYS
     74               walbedodat, wphisfi, wvolcapa,                          &
     75#else
     76               walbedodat, winertiedat, wphisfi,                       &
     77#endif
     78               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
     79               wtheta, wpsi)
    6480
    6581!

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq.inc

    r11 r28  
    1010               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    1111               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     12               wisoil, wdsoil, &
    1213               wecri_phys_sec,                &
    1314               output_tab2d, output_tab3d, flag_LES)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq1.inc

    r11 r28  
    1010               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    1111               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     12#ifdef NEWPHYS
     13               wisoil, wdsoil, &
     14#endif
    1215               wecri_phys_sec,                &
    1316               output_tab2d, output_tab3d, flag_LES)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq2.inc

    r11 r28  
    1313               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    1414               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     15#ifdef NEWPHYS
     16               wisoil, wdsoil, &
     17#endif
    1518               wecri_phys_sec,                &
    1619               output_tab2d, output_tab3d, flag_LES)
    … …  
    2225               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    2326               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     27#ifdef NEWPHYS
     28               wisoil, wdsoil, &
     29#endif
    2430               wecri_phys_sec,                &
    2531               output_tab2d, output_tab3d, flag_LES)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq3.inc

    r11 r28  
    1313               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    1414               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     15#ifdef NEWPHYS
     16               wisoil, wdsoil, &
     17#endif
    1518               wecri_phys_sec,                &
    1619               output_tab2d, output_tab3d, flag_LES)
    … …  
    2225               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    2326               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     27#ifdef NEWPHYS
     28               wisoil, wdsoil, &
     29#endif
    2430               wecri_phys_sec,                &
    2531               output_tab2d, output_tab3d, flag_LES)
    … …  
    3137               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    3238               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     39#ifdef NEWPHYS
     40               wisoil, wdsoil, &
     41#endif
    3342               wecri_phys_sec,                &
    3443               output_tab2d, output_tab3d, flag_LES)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq4.inc

    r11 r28  
    1313               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    1414               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     15#ifdef NEWPHYS
     16               wisoil, wdsoil, &
     17#endif
    1518               wecri_phys_sec,                &
    1619               output_tab2d, output_tab3d, flag_LES)
    … …  
    2225               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    2326               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     27#ifdef NEWPHYS
     28               wisoil, wdsoil, &
     29#endif
    2430               wecri_phys_sec,                &
    2531               output_tab2d, output_tab3d, flag_LES)
    … …  
    3137               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    3238               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     39#ifdef NEWPHYS
     40               wisoil, wdsoil, &
     41#endif
    3342               wecri_phys_sec,                &
    3443               output_tab2d, output_tab3d, flag_LES)
    … …  
    4049               pdu, pdv, pdt, pdq, pdpsrf, tracerdyn,          &
    4150               wtsurf,wtsoil,wemis,wq2,wqsurf,wco2ice,                 &
     51#ifdef NEWPHYS
     52               wisoil, wdsoil, &
     53#endif
    4254               wecri_phys_sec,                &
    4355               output_tab2d, output_tab3d, flag_LES)

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F

    r11 r28  
    284284              ENDDO 
    285285
     286
     287!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     288#ifdef NEWPHYS
     289              ! here, that bit is necessary for new soil model !         
     290              IF (config_flags%init_TI .ne. 0.) THEN
     291                  grid%therm_inert = grid%therm_inert*0. + config_flags%init_TI
     292                  print *, 'constant thermal inertia ', config_flags%init_TI
     293              ENDIF
     294
     295              DO k = 1, config_flags%num_soil_layers 
     296                grid%mars_isoil(i,k,j)=grid%em_isoil_gc(i,k+1,j)
     297                grid%mars_dsoil(i,k,j)=grid%em_dsoil_gc(i,k+1,j)
     298              ENDDO
     299
     300              DO k = 1, config_flags%num_soil_layers
     301                !!!!!!!!!!!!!!!!! DONE in soil_setting.F
     302                IF (grid%mars_dsoil(i,k,j) == -999.) THEN  !! old soil depths (or) no info in files
     303                   grid%mars_dsoil(i,k,j) = sqrt(887.75/3.14)*((2.**(k-0.5))-1.) * grid%therm_inert(i,j) / wvolcapa
     304                       !!! ATTENTION il faut interpoler si le nombre de niveaux change
     305                       !!! voir soil_setting.F (olddepthdef=.true. ; interpol=.true.)
     306                       !!! mais: en meso-echelle on a juste a prendre le mm nombre de niveaux que le GCM   
     307                ENDIF
     308                IF (grid%mars_isoil(i,k,j) == -999.) THEN  !! old soil model (or) no 3D thermal inertia   
     309                   grid%mars_isoil(i,k,j) = grid%therm_inert(i,j)
     310                ELSE             
     311                   IF (grid%mars_dsoil(i,k,j) .le. sqrt(88775./3.14) * grid%therm_inert(i,j) / wvolcapa) THEN
     312                           grid%mars_isoil(i,k,j) = grid%therm_inert(i,j)  !! if depth < skin depth, we use hi-res TI
     313                   ELSE  !! if depth > skin depth, we use low-res (GCM) TI
     314                         !! except for a transition layer
     315                         !!   EM: and, well, it would be wrong to sum up TI values
     316                         !!   EM: (cf. last page of soil model technical document)
     317                           IF (grid%mars_dsoil(i,k-1,j) .le. sqrt(88775./3.14) * grid%therm_inert(i,j) / wvolcapa) THEN
     318                                   grid%mars_isoil(i,k,j) = &
     319                                   sqrt( &
     320                                         ( grid%mars_dsoil(i,k+1,j) - grid%mars_dsoil(i,k-1,j) ) &
     321                                         / &
     322                                         ( ( (grid%mars_dsoil(i,k,j) - grid%mars_dsoil(i,k-1,j)) &
     323                                              / (grid%mars_isoil(i,k-1,j)*grid%mars_isoil(i,k-1,j)) ) &
     324                                           + &
     325                                           ( (grid%mars_dsoil(i,k+1,j) - grid%mars_dsoil(i,k,j)) &
     326                                              / (grid%mars_isoil(i,k+1,j)*grid%mars_isoil(i,k+1,j)) ) &
     327                                         ) &
     328                                       )
     329                           ENDIF       !! grid%mars_isoil(i,k-1,j) was changed at previous step to value grid%therm_inert(i,j)
     330                                       !! grid%mars_isoil(i,k+1,j) is defined to large-scale value grid%em_isoil_gc   
     331                   ENDIF       
     332                ENDIF
     333                !!!!!!!!!!!!!!!!! DONE in soil_setting.F
     334
     335              ENDDO 
     336#endif
     337!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     338
    286339grid%mars_gw(i,1,j)=grid%st040100(i,j)  !!ZMEA
    287340grid%mars_gw(i,2,j)=grid%st100200(i,j)  !!ZSTD
    … …  
    362415  DO j = jts, MIN(jte,jde-1)
    363416  DO i = its, MIN(ite,ide-1)
    364         grid%em_u_gc(i,:,j)=grid%em_u_gc(config_flags%init_WX,:,config_flags%init_WY)   ! zonal wind
    365         grid%em_v_gc(i,:,j)=grid%em_v_gc(config_flags%init_WX,:,config_flags%init_WY)   ! meridional wind
     417     grid%em_u_gc(i,:,j)=grid%em_u_gc(config_flags%init_WX,:,config_flags%init_WY) ! zonal wind
     418     grid%em_v_gc(i,:,j)=grid%em_v_gc(config_flags%init_WX,:,config_flags%init_WY) ! meridional wind
    366419  ENDDO
    367420  ENDDO
    368421  !! FIX for the STAGGERED SPECIFICITY
    369422  grid%em_u_gc(MIN(ite,ide-1)+1,:,:)=grid%em_u_gc(MIN(ite,ide-1),:,:)
    370   grid%em_v_gc(:,:,MIN(jte,jde-1)+1)=grid%em_v_gc(:,:,MIN(jte,jde-1))   
     423  grid%em_v_gc(:,:,MIN(jte,jde-1)+1)=grid%em_v_gc(:,:,MIN(jte,jde-1))
    371424
    372425  !! CHECK
    … …  
    396449IF (config_flags%init_LES) THEN
    397450
    398         print *, '*** LES MODE ***'
    399         print *, 'setting uniform values and profiles'
     451        print *, '*** LES MODE ***'
     452        print *, 'setting uniform values and profiles'
    400453        print *, 'u', grid%em_u_gc(its+1,:,jts+1)
    401454        print *, 'v', grid%em_v_gc(its+1,:,jts+1)
    … …  
    430483grid%mars_tsoil(i,:,j)=grid%mars_tsoil(its+1,:,jts+1)
    431484
     485!!!
     486grid%mars_isoil(i,:,j)=grid%mars_isoil(its+1,:,jts+1)
     487grid%mars_dsoil(i,:,j)=grid%mars_dsoil(its+1,:,jts+1)
    432488
    433489!! T.Michaels trick to break symmetry

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F

    r11 r28  
    854854     &        ,MARS_CICE=grid%mars_cice,MARS_EMISS=grid%mars_emiss   & 
    855855     &        ,MARS_TSOIL=grid%mars_tsoil       &
     856#ifdef NEWPHYS
     857     &        ,MARS_ISOIL=grid%mars_isoil       &
     858     &        ,MARS_DSOIL=grid%mars_dsoil       &
     859#endif
    856860     &        ,MARS_GW=grid%mars_gw  & 
    857861     &        ,NUM_SOIL_LAYERS=config_flags%num_soil_layers    &

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

    r11 r28  
    1010!           - general cleaning of code and comments - October 2008
    1111!           - additions for idealized cases - January 2009
     12!           - additions for new soil model in physics - January 2010
    1213!*******************************************************************************
    1314MODULE module_lmd_driver
    … …  
    2930        MARS_ALB,MARS_TI,MARS_CICE,MARS_EMISS, &
    3031        MARS_TSOIL, &
     32#ifdef NEWPHYS
     33        MARS_ISOIL, &
     34        MARS_DSOIL, &
     35#endif
    3136        MARS_GW, &
    3237        NUM_SOIL_LAYERS, &
    … …  
    105110!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    106111REAL, DIMENSION( ims:ime, NUM_SOIL_LAYERS, jms:jme ), INTENT(IN   )  :: &
    107      MARS_TSOIL         
     112     MARS_TSOIL
     113#ifdef NEWPHYS
     114REAL, DIMENSION( ims:ime, NUM_SOIL_LAYERS, jms:jme ), INTENT(IN   )  :: &
     115     MARS_ISOIL, MARS_DSOIL         
     116#endif
    108117REAL, DIMENSION( ims:ime, 5, jms:jme ), INTENT(IN   )  :: &
    109118     MARS_GW
    … …  
    147156   LOGICAL :: firstcall,lastcall,tracerdyn
    148157   REAL,DIMENSION(:),ALLOCATABLE :: q2_val, qsurf_val, tsoil_val
     158#ifdef NEWPHYS
     159   REAL,DIMENSION(:),ALLOCATABLE :: isoil_val, dsoil_val
     160#endif
    149161   REAL,DIMENSION(:),ALLOCATABLE :: aire_vec,lat_vec,lon_vec
    150162   REAL,DIMENSION(:),ALLOCATABLE :: walbedodat,winertiedat,wphisfi
    … …  
    153165   ! v--- can they be modified ?
    154166   REAL,DIMENSION(:),ALLOCATABLE :: wtsurf,wco2ice,wemis
    155    REAL,DIMENSION(:,:),ALLOCATABLE :: wq2,wqsurf,wtsoil
     167   REAL,DIMENSION(:,:),ALLOCATABLE :: wq2,wqsurf,wtsoil
     168#ifdef NEWPHYS
     169   REAL,DIMENSION(:,:),ALLOCATABLE :: wisoil,wdsoil
     170#endif
    156171   ! ----------
    157172   REAL,DIMENSION(:,:),ALLOCATABLE :: pplev,pplay,pphi,pu,pv,pt,pw
    … …  
    451466ALLOCATE(qsurf_val(nq))
    452467ALLOCATE(tsoil_val(nsoil))
     468#ifdef NEWPHYS
     469ALLOCATE(isoil_val(nsoil))
     470ALLOCATE(dsoil_val(nsoil))
     471#endif
    453472ALLOCATE(wq2(ngrid,nlayer+1))
    454473ALLOCATE(wqsurf(ngrid,nq))
    455474ALLOCATE(wtsoil(ngrid,nsoil))
     475#ifdef NEWPHYS
     476ALLOCATE(wisoil(ngrid,nsoil))
     477ALLOCATE(wdsoil(ngrid,nsoil))
     478#endif
    456479ALLOCATE(pplev(ngrid,nlayer+1))
    457480ALLOCATE(pplay(ngrid,nlayer))
    … …  
    594617tsoil_val = tsoil_val*0. + tsurf_val
    595618ENDIF
     619#ifdef NEWPHYS
     620isoil_val(:)=MARS_ISOIL(i,:,j)
     621dsoil_val(:)=MARS_DSOIL(i,:,j)
     622#endif
    596623
    597624!-------------------!
    … …  
    630657 PRINT *,'v_prof ',v_prof
    631658 PRINT *,'tsoil ',tsoil_val
     659#ifdef NEWPHYS
     660 PRINT *,'isoil ',isoil_val
     661 PRINT *,'dsoil ',dsoil_val
     662#endif
    632663ENDIF
    633664
    … …  
    660691wqsurf(subs,:) = qsurf_val(:)
    661692wtsoil(subs,:) = tsoil_val(:)
     693#ifdef NEWPHYS
     694wisoil(subs,:) = isoil_val(:)
     695wdsoil(subs,:) = dsoil_val(:)
     696#endif
    662697wtheta(subs) = theta_val
    663698wpsi(subs) = psi_val
    … …  
    680715pplev(subs,1:nlayer) = p8w_prof(1:nlayer)  !! NB: last level: no data
    681716pplay(subs,:) = p_prof(:)
    682 pt(subs,:) = t_prof(:)  ! temperature
     717pt(subs,:) = t_prof(:)
    683718pu(subs,:) = u_prof(:)
    684719pv(subs,:) = v_prof(:)
    … …  
    711746DEALLOCATE(qsurf_val)
    712747DEALLOCATE(tsoil_val)
     748#ifdef NEWPHYS
     749DEALLOCATE(isoil_val)
     750DEALLOCATE(dsoil_val)
     751#endif
    713752DEALLOCATE(dz8w_prof)
    714753DEALLOCATE(z_prof)
    … …  
    738777IF (firstcall .EQV. .true.) THEN
    739778print *, '** Mars ** LMD INITIALIZATION'
    740 include "../call_meso_inifis.inc"
     779#include "../call_meso_inifis.inc"
     780!!! le # est important pour newphys
    741781DEALLOCATE(aire_vec)
    742782DEALLOCATE(lat_vec)
    … …  
    777817pdt(:,:)=0.
    778818pdq(:,:,:)=0.
    779 include "../call_meso_physiq.inc"
     819#include "../call_meso_physiq.inc"
     820!!! le # est important pour newphys
    780821DEALLOCATE(pplev)
    781822DEALLOCATE(pplay)
    … …  
    792833DEALLOCATE(wqsurf)
    793834DEALLOCATE(wtsoil)
     835#ifdef NEWPHYS
     836DEALLOCATE(wisoil)
     837DEALLOCATE(wdsoil)
     838#endif
    794839
    795840

  • trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/share/module_model_constants.F

    r11 r28  
    7272   REAL    , PARAMETER :: wdtemisiceN   = 0.4  ! time scale for snow metamorphism (north) !2 
    7373   REAL    , PARAMETER :: wdtemisiceS   = 0.4  ! time scale for snow metamorphism (south) !2
     74#ifdef NEWPHYS
     75   REAL    , PARAMETER :: wvolcapa      = 1.e6       ! volumetric capacity of soil (new soil model)
     76#endif
    7477!
    7578!-------------------------------

  • trunk/mesoscale/LMD_MM_MARS/makemeso

    r24 r28  
    299299                       sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT"+g > yeah
    300300                                #| sed s+"fastsse"+"fast"+g > yeah
    301                        mv -f yeah configure.wrf ;;
     301                       mv -f yeah configure.wrf ;;
     302                       #sed s+"-fastsse"+"-O2 -Munroll -Mcache_align"+g configure.wrf > yeah
     303                       ##sed s+"-fastsse"+"-O2 -fpic"+g configure.wrf > yeah ## marche pas, ILM problem
     304                       ##sed s+"-fastsse"+"-mcmodel=medium -Mlarge_arrays"+g configure.wrf > yeah ## marche pas, ILM problem
     305                       #mv -f yeah configure.wrf ;;
    302306                     # MPICH, 64 bits, OK with periodic BC but no positive definite
    303307     mpi_64_test)      cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd ..
    … …  
    372376  # -- default is 10
    373377  soilsize=10
     378  #soilsize=18 ## nouvelle physique
    374379
    375380  # GCM environment variables
Note: See TracChangeset for help on using the changeset viewer.