source: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis2.inc @ 1634

Last change on this file since 1634 was 1212, checked in by aslmd, 12 years ago

LMDZ.MARS + MESOSCALE

A quite major commit, at least for MESOSCALE.
In a word: any ngrid deserves to be free.

  • no need to recompile when changing number of horizontal grid points or number of processors
  • latest version of LMDZ.MARS physics can be used
  • WARNING! Nesting is still yet to be fixed (since r1027)

Also some small bug fixes to LMDZ.MARS.

Changes in LMDZ.MARS

--> fixed a potential bug in thermal plume model because zlmax was computed both in thermcell_main_mars and calltherm_interface... so made it an OUT argument of calltherm_interface. also: changed the name to limz. and added precompiling flags to avoid the use of planetwide in MESOSCALE. in MESOSCALE we just go high enough (nlayer-5) and do not care about computational cost (although we certainly gain from not using MAXVAL).
--> moved allocations upward in inifis. does not change anything for GCM, but make MESOSCALE modifications simpler, and overall make inifis better organized: first allocations, then reading callphys.def file.
--> added precompiling flags around lines that are both useless for MESOSCALE (notably I/O) and recently adapted to parallel computations in the GCM
--> tidied up what is MESOSCALE vs. GCM in surfini

Changes in MESOSCALE

--> changed makemeso to allow dynamically set nx ny nprocs
--> changed makemeso to remove links to Fortran code adapted to parallel GCM and useless for mesoscale
--> changed ngridmx to ngrid in inifis includes

File size: 2.3 KB
Line 
1!
2! part of the module_lmd_driver program
3! that is (number of domains)-dependent
4!
5
6
7SELECT CASE (id)
8
9       CASE(1)
10       CALL inifis(ngrid,nlayer, &
11#ifdef NEWPHYS
12               nq,  &
13#endif
14               wday_ini,wdaysec,                &
15               wappel_phys,                     &
16               lat_vec,lon_vec,aire_vec,1/reradius,g,r_d,cp,           &
17#ifdef NEWPHYS
18               dt,                                                     &
19#else
20               nq,dt,                                                &
21#endif
22               womeg,wmugaz,                                           &
23               wyear_day,wperiheli,waphelie,wperi_day,wobliquit,       &
24               wz0,wemin_turb,wlmixmin,                                &
25               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
26               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
27#ifdef NEWPHYS
28               walbedodat, wphisfi, wvolcapa, wz0tab,                  &
29#else
30               walbedodat, winertiedat, wphisfi,                       &
31#endif
32               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
33               wtheta, wpsi)
34
35
36       CASE(2)
37       CALL d2_inifis(ngrid,nlayer, &
38#ifdef NEWPHYS
39               nq,  &
40#endif
41               wday_ini,wdaysec,                   &
42               wappel_phys,                        &
43               lat_vec,lon_vec,aire_vec,1/reradius,g,r_d,cp,           &
44#ifdef NEWPHYS
45               dt,                                                     &
46#else
47               nq,dt,                                                &
48#endif
49               womeg,wmugaz,                                           &
50               wyear_day,wperiheli,waphelie,wperi_day,wobliquit,       &
51               wz0,wemin_turb,wlmixmin,                                &
52               wemissiv,wemissiceN,wemissiceS,walbediceN,walbediceS,   &
53               wiceradiusN,wiceradiusS,wdtemisiceN,wdtemisiceS,        &
54#ifdef NEWPHYS
55               walbedodat, wphisfi, wvolcapa, wz0tab,                  &
56#else
57               walbedodat, winertiedat, wphisfi,                       &
58#endif
59               wzmea,wzstd,wzsig,wzgam,wzthe,                          &
60               wtheta, wpsi)
61
62!
63!       CASE(3:)
64!       print *, '** Mars ** Not supported ...'
65!       stop
66!
67END SELECT
Note: See TracBrowser for help on using the repository browser.