source: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_ini.F @ 268

Last change on this file since 268 was 226, checked in by aslmd, 13 years ago

MESOSCALE/LMDZ.MARS.new
--> modified to impact last changes

MESOSCALE/LMD_MM_MARS/makemeso
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq?.inc
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis?.inc
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
--> modified to get rid of ecri_phys

and make changes related to meso_physiq and meso_inifis

LMDZ.MARS/libf/phymars
--> see LMDZ.MARS/README

15/07/2011 == AS

  • Modified the mesoscale part so that the previous change by EM does not imply an error in the mesoscale case. More development is needed though to get the "varying z0" capability in the mesoscale model.
  • Worked on versions of meso_physiq and meso_inifis as close as possible to physiq and inifis for more continuity in the process of impacting changes (and even possibly to reach a common version of physiq and inifis).

    The main point is to make the mesoscale significant specific parts

    coded into include files in meso_inc so that meso_physiq and meso_inifis looks very close to physiq and inifis.

    This is completely transparent for GCM users who does not need the

    contents of meso_inc.
  • Slight cosmetic changes to physiq.f and inifis.F --- some of them e.g. to prepare convergence between meso_physiq and physiq
File size: 2.6 KB
Line 
1ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
2c MESOSCALE MODEL
3c
4c       No need to use startfi.nc
5c               > part of the job of phyetat0 is done in inifis
6c               > remaining initializations are passed here from the WRF variables
7c               > beware, some operations were done by phyetat0 (ex: tracers)
8c                       > if any problems, look in phyetat0
9c
10c    LMD_MM_MARS
11c       day_ini,tsurf,tsoil,emis,q2,qsurf,co2ice are inputs
12c               and locally saved variables
13c                       (no need to call phyetat0)
14c
15      tsurf(:)=wtsurf(:)
16      PRINT*,'check: tsurf ',tsurf(1),tsurf(ngridmx)
17      tsoil(:,:)=wtsoil(:,:)
18      PRINT*,'check: tsoil ',tsoil(1,1),tsoil(ngridmx,nsoilmx)
19     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
20     !!!new physics
21      PRINT*,'no variable z0 field for the moment. use z0_default.'
22      z0(:) = z0_default
23      PRINT*, 'check z0 ', z0(1)
24      inertiedat(:,:)=wisoil(:,:)
25      PRINT*,'check: inert ',inertiedat(1,1),inertiedat(ngridmx,nsoilmx)
26      mlayer(0:nsoilmx-1)=wdsoil(1,:)
27      PRINT*,'check: midlayer ', mlayer(:)
28            !!!!!!!!!!!!!!!!! DONE in soil_setting.F
29            ! 1.5 Build layer(); following the same law as mlayer()
30            ! Assuming layer distribution follows mid-layer law:
31            ! layer(k)=lay1*alpha**(k-1)
32            lay1=sqrt(mlayer(0)*mlayer(1))
33            alpha=mlayer(1)/mlayer(0)
34            do iloop=1,nsoilmx
35              layer(iloop)=lay1*(alpha**(iloop-1))
36            enddo
37
38      PRINT*,'check: layer ', layer(:)
39
40            !!!!!!!!!!!!!!!!! DONE in soil_setting.F
41      tnom(:)=wtnom(:)   !! est rempli dans advtrac.h
42      PRINT*,'check: tracernames ', tnom
43     !!!new physics
44     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
45      emis(:)=wemis(:)
46      PRINT*,'check: emis ',emis(1),emis(ngridmx)
47      q2(:,:)=wq2(:,:)
48      PRINT*,'check: q2 ',q2(1,1),q2(ngridmx,nlayermx+1)
49      qsurf(:,:)=wqsurf(:,:)
50      PRINT*,'check: qsurf ',qsurf(1,1),qsurf(ngridmx,nqmx)
51      co2ice(:)=wco2ice(:)
52      PRINT*,'check: co2 ',co2ice(1),co2ice(ngridmx)
53      day_ini=wday_ini
54
55c       artificially filling dyn3d/control.h is also required
56c       > iphysiq is put in WRF to be set easily (cf ptimestep)
57c       > day_step is simply deduced:
58c
59      day_step=daysec/ptimestep
60      PRINT*,'Call to LMD physics:',day_step,' per Martian day'
61c
62      iphysiq=ptimestep
63c
64      !DO iq=1, nq
65      !  PRINT*, tnom(iq), pq(:,:,iq)
66      !ENDDO
67
68c
69ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Note: See TracBrowser for help on using the repository browser.