Changeset 667 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

May 23, 2012, 4:32:58 PM (13 years ago)

Author:

acolaitis

Message:

ADAPTED RESTART IN MESOSCALE AND LES.
Warning: MMM always puts the initialization state as the first frame in outputs :

init,frame1,frame2,...,frameN

Hence, when using runs restarting from one another, one will get:

init,frame1,frame2,...,frameN
frameN,frameN+1,...,frame2N

Dont forget to remove the initialization frame before concatenating your netcdfs !!

HOW TO USE THE RESTART IN LES AND MESOSCALE ?
Exemple :
=========

---

INITIALIZE FIRST RUN (ideal.exe or real.exe)
FIRST RUN

---

---

FIRST RUN

---

EXTRACT FROM namelist.input:

&time_control
run_seconds = 60 !! Duration of simulation [LES: !=0]
history_interval_s = 10 !! Frequency of outputs in seconds [LES: !=0]
frames_per_outfile = 37 !! Size of time dimension in files
start_year = 9999 !! Idealized simulation [LES: 9999]
start_month = 01
start_day = 01
start_hour = 00
start_minute = 00
start_second = 00
end_year = 9999 !! Idealized simulation [LES: 9999]
restart = .false. !! (*) whether this run is a restart run
restart_interval = 1 !! (*) restart output file interval in minutes
io_form_history = 2 !! (*) Choice of NETCDF for ouputs
io_form_restart = 2 !! (*) Choice of NETCDF for ouputs
io_form_input = 2 !! (*) Choice of NETCDF for ouputs
io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs
debug_level = 0 !! (*) Verbose level
/

This will produce a restart file : wrfrst_d01_9999-01-01_00:01:00
(and a regular output file)

---

SECOND RUN, RESTART FROM FIRST
YOU DONT NEED TO RECOMPUTE START FILES (no need for ideal.exe or real.exe)

---

EXTRACT FROM namelist.input:

&time_control
run_seconds = 60 !! Duration of simulation [LES: !=0]
history_interval_s = 10 !! Frequency of outputs in seconds [LES: !=0]
frames_per_outfile = 37 !! Size of time dimension in files
start_year = 9999 !! Idealized simulation [LES: 9999]
start_month = 01
start_day = 01
start_hour = 00
start_minute = 01
start_second = 00
end_year = 9999 !! Idealized simulation [LES: 9999]
restart = .true. !! (*) whether this run is a restart run
restart_interval = 1 !! (*) restart output file interval in minutes
io_form_history = 2 !! (*) Choice of NETCDF for ouputs
io_form_restart = 2 !! (*) Choice of NETCDF for ouputs
io_form_input = 2 !! (*) Choice of NETCDF for ouputs
io_form_boundary = 2 !! (*) Choice of NETCDF for ouputs
debug_level = 0 !! (*) Verbose level
/

This requires wrfrst_d01_9999-01-01_00:01:00 to be present
This will produce a restart file : wrfrst_d01_9999-01-01_00:02:00

---

etc ad infinitum !

---

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (4 diffs)

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

@@ -233 +233 @@
    REAL :: wappel_phys
    LOGICAL :: flag_LES
-
+   LOGICAL, SAVE :: flag_first_restart
 !**************************************************
 ! IMPORTANT: pour travailler avec grid nesting
@@ -359 +359 @@
 !----------------------------------------------!
 IF (elaps .eq. 0.) THEN
+   flag_first_restart = .false.
+ELSE
+   flag_first_restart=flag_first_restart.OR.(.NOT. ALLOCATED(dp_save))
+ENDIF
+
+IF ((elaps .eq. 0.).or.flag_first_restart) THEN
 firstcall=.true.  !! for continuity with GCM, physics are always called at the first WRF timestep
   !firstcall=.false. !! just in case you'd want to get rid of the physics 
@@ -367 +373 @@
    PRINT *, '**** check **** OK I ALLOCATE one save ARRAY for all NESTS ', max_dom, id
    !! here are the arrays that would be useful to save physics tendencies
-   ALLOCATE(dp_save(ngrid,max_dom)) 
+   ALLOCATE(dp_save(ngrid,max_dom))
    ALLOCATE(du_save(ngrid,nlayer,max_dom))
-   ALLOCATE(dv_save(ngrid,nlayer,max_dom)) 
-   ALLOCATE(dt_save(ngrid,nlayer,max_dom))     
-   ALLOCATE(dq_save(ngrid,nlayer,nq,max_dom)) 
+   ALLOCATE(dv_save(ngrid,nlayer,max_dom))
+   ALLOCATE(dt_save(ngrid,nlayer,max_dom))
+   ALLOCATE(dq_save(ngrid,nlayer,nq,max_dom))
    dp_save(:,:)=0.    !! initialize these arrays ...
    du_save(:,:,:)=0.
@@ -378 +384 @@
    dq_save(:,:,:,:)=0.
 ENDIF
+IF (id .lt. max_dom) THEN
+   flag_first_restart=.true.
+ELSE
+   flag_first_restart=.false.
+ENDIF
 #else
 ALLOCATE(dp_save(ngrid)) !! here are the arrays that would be useful to save physics tendencies
 ALLOCATE(du_save(ngrid,nlayer))
-ALLOCATE(dv_save(ngrid,nlayer)) 
-ALLOCATE(dt_save(ngrid,nlayer))     
-ALLOCATE(dq_save(ngrid,nlayer,nq)) 
+ALLOCATE(dv_save(ngrid,nlayer))
+ALLOCATE(dt_save(ngrid,nlayer))
+ALLOCATE(dq_save(ngrid,nlayer,nq))
 dp_save(:)=0.    !! initialize these arrays ...
 du_save(:,:)=0.
 dv_save(:,:)=0.
-dt_save(:,:)=0.    
+dt_save(:,:)=0.
 dq_save(:,:,:)=0.
+flag_first_restart=.false.
 #endif
 !! put here some general information you'd like to print just once
     print *, 'TILES: ', i_start,i_end, j_start, j_end  ! numbers for simple runs, arrays for parallel runs
-    print *, 'DOMAIN: ', ids, ide, jds, jde 
+    print *, 'DOMAIN: ', ids, ide, jds, jde
     print *, 'MEMORY: ', ims, ime, jms, jme
     print *, 'ADVECTED TRACERS: ', num_3d_s-1