Ignore:
Timestamp:
May 25, 2012, 8:06:39 PM (13 years ago)
Author:
acolaitis
Message:

LES RESTART
===========

Several corrections for LES restart. Added a number of save variables from physiq to restart netcdf, so that continuity between runs can truly be achieved.
Added some lines in makemeso so that debug option "-g" also works for WRF (in addtion to the GCM physics) when compiling with ifort

===============================
WARNING WARNING WARNING WARNING
===============================

  • FROM THIS REVISION, YOU MUST MODIFY MANUALLY SOME FILES AS FOLLOW, BEFORE YOU CAN RE COMPILE THE LES :

simply copy the call_meso_physiq*.inc files in $MMM/SRC/WRFV2/ into $MMM/SRC/LES/WRFV2/.
=>> this is usually done during LES installation by $MMM/SRC/LES/LMD_LES_MARS_install

  • IF YOU WANT TO MAKE A RUN WITH RESTART FILES, YOU MUST MAKE A CLEAN RECOMPILATION OF THE MODEL WITH THE UPDATED REGISTRY AND UPDATED GCM FILES

=>> it is advised to remove your working directory (ex: lesnewphys_mpifort64) and start again with makemeso

  • IF YOU WANT TO MAKE A RUN AND GENERATE RESTART FILES, YOU MUST RECOMPILE THE IDEAL.EXE STEP AND PERFORM IT

=>> current wrfinput will not work, as it does not contain the new variables
===============================
WARNING WARNING WARNING WARNING
===============================

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r667 r674  
    3232        MARS_WICE, &
    3333        MARS_TSOIL, &
     34        MARS_Q2, &
     35        MARS_FLUXRAD, &
     36        MARS_WSTAR, &
     37        MARS_TSURF, &
    3438#ifdef NEWPHYS
    3539        MARS_ISOIL, &
    … …  
    5357#include "module_lmd_driver_output1.inc"
    5458#endif
    55         SLPX,SLPY)
     59        SLPX,SLPY,RESTART)
    5660! NB: module_lmd_driver_output1.inc : output arguments generated from Registry
    5761
    … …  
    115119     MSFT,MSFU,MSFV, &
    116120     XLAT,XLONG,HT,  &
    117      MARS_ALB,MARS_TI,MARS_EMISS,MARS_CICE, &
    118      MARS_WICE, &
     121     MARS_ALB,MARS_TI,MARS_EMISS, &
    119122     SLPX,SLPY
     123REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT   )  :: &
     124     MARS_CICE,MARS_WICE, &
     125     MARS_FLUXRAD,MARS_WSTAR,MARS_TSURF
    120126! 3D arrays
    121127REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ) :: &
    122128     dz8w,p8w,p,exner,t,t8w,rho,u,v,w,z,th
     129REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
     130     MARS_Q2
    123131!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    124132REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
    … …  
    126134INTEGER, INTENT(IN   ) :: HISTORY_INTERVAL
    127135!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    128 REAL, DIMENSION( ims:ime, NUM_SOIL_LAYERS, jms:jme ), INTENT(IN  )  :: &
     136REAL, DIMENSION( ims:ime, NUM_SOIL_LAYERS, jms:jme ), INTENT(INOUT )  :: &
    129137     MARS_TSOIL
    130138#ifdef NEWPHYS
    … …  
    141149REAL, DIMENSION( ims:ime, kms:kme, jms:jme, 1:num_3d_s ), INTENT(INOUT ) :: &
    142150     scalar
     151! Logical
     152LOGICAL, INTENT(IN ) :: restart
    143153
    144154!-------------------------------------------
    … …  
    187197   REAL :: zmea_val,zstd_val,zsig_val,zgam_val,zthe_val
    188198   REAL :: theta_val, psi_val
     199   REAL :: wstar_val,fluxrad_val
    189200   LOGICAL :: firstcall,lastcall,tracerdyn
    190201   REAL,DIMENSION(:),ALLOCATABLE :: q2_val, qsurf_val, tsoil_val
    … …  
    198209   REAL,DIMENSION(:),ALLOCATABLE :: wtheta, wpsi
    199210   ! v--- can they be modified ?
    200    REAL,DIMENSION(:),ALLOCATABLE :: wtsurf,wco2ice,wemis
     211   REAL,DIMENSION(:),ALLOCATABLE :: wtsurf,wco2ice,wemis,wwstar,wfluxrad
    201212   REAL,DIMENSION(:,:),ALLOCATABLE :: wq2,wqsurf,wtsoil
    202213#ifdef NEWPHYS
    … …  
    257268#endif
    258269
     270!! FOR RESTART
     271      REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: &
     272             save_tsoil_restart
     273      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: &
     274             save_tsurf_restart
     275      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: &
     276             save_fluxrad_restart
     277      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: &
     278             save_co2ice_restart
     279      REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: &
     280             save_q2_restart
     281      REAL, DIMENSION(:,:), ALLOCATABLE, SAVE :: &
     282             save_qsurf_restart
     283      REAL, DIMENSION(:), ALLOCATABLE, SAVE :: &
     284             save_wstar_restart
     285
    259286!!!IDEALIZED IDEALIZED
    260287      REAL :: lat_input, lon_input, ls_input, lct_input
    … …  
    402429flag_first_restart=.false.
    403430#endif
     431
     432! Restart save arrays
     433ALLOCATE(save_tsoil_restart(ngrid,nsoil))
     434ALLOCATE(save_fluxrad_restart(ngrid))
     435ALLOCATE(save_co2ice_restart(ngrid))
     436ALLOCATE(save_q2_restart(ngrid,nlayer+1))
     437ALLOCATE(save_qsurf_restart(ngrid,nq))
     438ALLOCATE(save_wstar_restart(ngrid))
     439ALLOCATE(save_tsurf_restart(ngrid))
     440save_tsoil_restart(:,:)=0.
     441save_fluxrad_restart(:)=0.
     442save_co2ice_restart(:)=0.
     443save_q2_restart(:,:)=0.
     444save_qsurf_restart(:,:)=0.
     445save_wstar_restart(:)=0.
     446save_tsurf_restart(:)=0.
     447
    404448!! put here some general information you'd like to print just once
    405449    print *, 'TILES: ', i_start,i_end, j_start, j_end  ! numbers for simple runs, arrays for parallel runs
    … …  
    622666ALLOCATE(wqsurf(ngrid,nq))       !!!!!
    623667ALLOCATE(wtsoil(ngrid,nsoil))    !!!!!
     668ALLOCATE(wfluxrad(ngrid))
     669ALLOCATE(wwstar(ngrid))
    624670#ifdef NEWPHYS
    625671ALLOCATE(wisoil(ngrid,nsoil))    !!!!!
    … …  
    769815  !!! CAS DU CO2
    770816  DO iii=1,nq
    771    IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
     817   IF ( wtnom(iii) .eq. 'co2' .and. (.not. restart)) q_prof(:,iii) = 0.95
    772818  ENDDO
    773819  IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
    774    IF (firstcall .EQV. .true.) THEN
     820   IF (firstcall .EQV. .true. .and. (.not. restart)) THEN
    775821      q_prof(:,:) = 0.95
    776822   ENDIF
    … …  
    892938! Ground temperature, emissivity, CO2 ice cover !
    893939!-----------------------------------------------!
    894 tsurf_val=tsk(i,j)
     940IF (.not. restart) THEN
     941   tsurf_val=tsk(i,j) ! because tsk is not updated in lmd_driver, the stored value in restart is at t, not t+dt
     942                      ! note to aslmd: i dont know why we dont update tsk, so I use MARS_TSURF to be conservative
     943ELSE
     944   tsurf_val=MARS_TSURF(i,j)
     945ENDIF
    895946emis_val=MARS_EMISS(i,j)
    896947co2ice_val=MARS_CICE(i,j)
    … …  
    912963  IF ( (i == ips) .AND. (j == jps) ) PRINT *,'** Mars ** IDEALIZED SIMULATION tsoil is set to tsurf'
    913964  do k=1,nsoil
    914    tsoil_val(k) = tsurf_val
     965   IF (.not.restart) THEN
     966     tsoil_val(k) = tsurf_val
     967   ELSE
     968     !this is a restart run. We must not set tsoil to tsurf in the init.
     969     !tsoil was saved in physiq.F under the name MARS_TSOIL in the restart file
     970     !(see Registry)
     971     tsoil_val(k)=MARS_TSOIL(i,k,j)
     972   ENDIF
     973
    915974#ifdef NEWPHYS
    916975   IF ( nsoil .lt. 18 ) THEN
    … …  
    10281087!-------------------------!
    10291088!-------------------------!
    1030 q2_val(:)=0.      !PBL wind variance
     1089IF (.not. restart) THEN
     1090   q2_val(:)=0.      !PBL wind variance
     1091   fluxrad_val=0.
     1092   wstar_val=0.
     1093ELSE
     1094   q2_val(:)=MARS_Q2(i,:,j)
     1095   fluxrad_val=MARS_FLUXRAD(i,j)
     1096   wstar_val=MARS_WSTAR(i,j)
     1097ENDIF
    10311098
    10321099!-----------------!
    … …  
    10481115wemis(subs) = emis_val
    10491116wq2(subs,:) = q2_val(:)
     1117wfluxrad(subs) = fluxrad_val
     1118wwstar(subs) = wstar_val
    10501119wqsurf(subs,:) = qsurf_val(:)
    10511120wtsoil(subs,:) = tsoil_val(:)
    … …  
    12061275DEALLOCATE(pw)
    12071276DEALLOCATE(pq)
    1208 DEALLOCATE(wtsurf)
    1209 DEALLOCATE(wco2ice)
     1277!DEALLOCATE(wtsurf)
     1278!DEALLOCATE(wco2ice)
    12101279DEALLOCATE(wemis)
    1211 DEALLOCATE(wq2)
    1212 DEALLOCATE(wqsurf)
    1213 DEALLOCATE(wtsoil)
     1280!DEALLOCATE(wq2)
     1281!DEALLOCATE(wqsurf)
     1282!DEALLOCATE(wtsoil)
     1283!DEALLOCATE(wwstar)
     1284!DEALLOCATE(wfluxrad)
    12141285#ifdef NEWPHYS
    12151286DEALLOCATE(wisoil)
    … …  
    12711342dq_save(:,:,:)=pdq(:,:,:)
    12721343#endif
    1273 
     1344save_tsoil_restart(:,:)=wtsoil(:,:)
     1345save_fluxrad_restart(:)=wfluxrad(:)
     1346save_co2ice_restart(:)=wco2ice(:)
     1347save_q2_restart(:,:)=wq2(:,:)
     1348save_qsurf_restart(:,:)=wqsurf(:,:)
     1349save_wstar_restart(:)=wwstar(:)
     1350save_tsurf_restart(:)=wtsurf(:)
     1351DEALLOCATE(wtsoil)
     1352DEALLOCATE(wfluxrad)
     1353DEALLOCATE(wco2ice)
     1354DEALLOCATE(wq2)
     1355DEALLOCATE(wqsurf)
     1356DEALLOCATE(wwstar)
     1357DEALLOCATE(wtsurf)
    12741358ENDIF call_physics
    12751359
    … …  
    13511435!---------------------------!
    13521436PSFC(i,j)=PSFC(i,j)+pdpsrf(subs)*dt    !!! here dt is needed
     1437
     1438!------------------------------------!
     1439! Save key variables for restart ! 
     1440!------------------------------------!
     1441
     1442MARS_TSOIL(i,:,j)=save_tsoil_restart(subs,:)
     1443MARS_FLUXRAD(i,j)=save_fluxrad_restart(subs)
     1444MARS_CICE(i,j)=save_co2ice_restart(subs)
     1445MARS_Q2(i,:,j)=save_q2_restart(subs,:)
     1446SELECT CASE (MARS_MODE)
     1447   CASE (1,11,12)
     1448     MARS_WICE(i,j)=save_qsurf_restart(subs,2)
     1449END SELECT
     1450MARS_WSTAR(i,j)=save_wstar_restart(subs)
     1451MARS_TSURF(i,j)=save_tsurf_restart(subs)
    13531452
    13541453!---------!
Note: See TracChangeset for help on using the changeset viewer.