ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c MESOSCALE MODEL
c 
c       No need to use startfi.nc
c               > part of the job of phyetat0 is done in inifis
c               > remaining initializations are passed here from the WRF variables
c               > beware, some operations were done by phyetat0 (ex: tracers)
c                       > if any problems, look in phyetat0
c
c    LMD_MM_MARS
c       day_ini,tsurf,tsoil,emis,q2,qsurf,co2ice are inputs
c               and locally saved variables
c                       (no need to call phyetat0)
c
      tsurf(:)=wtsurf(:)
      PRINT*,'check: tsurf ',tsurf(1),tsurf(ngridmx)
      tsoil(:,:)=wtsoil(:,:)
      PRINT*,'check: tsoil ',tsoil(1,1),tsoil(ngridmx,nsoilmx)
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     !!!new physics
      PRINT*,'no variable z0 field for the moment. use z0_default.'
      z0(:) = z0_default
      PRINT*, 'check z0 ', z0(1)
      inertiedat(:,:)=wisoil(:,:)
      PRINT*,'check: inert ',inertiedat(1,1),inertiedat(ngridmx,nsoilmx)
      mlayer(0:nsoilmx-1)=wdsoil(1,:)
      PRINT*,'check: midlayer ', mlayer(:)
            !!!!!!!!!!!!!!!!! DONE in soil_setting.F 
            ! 1.5 Build layer(); following the same law as mlayer()
            ! Assuming layer distribution follows mid-layer law:
            ! layer(k)=lay1*alpha**(k-1)
            lay1=sqrt(mlayer(0)*mlayer(1))
            alpha=mlayer(1)/mlayer(0)
            do iloop=1,nsoilmx
              layer(iloop)=lay1*(alpha**(iloop-1))
            enddo

      PRINT*,'check: layer ', layer(:)

            !!!!!!!!!!!!!!!!! DONE in soil_setting.F
      tnom(:)=wtnom(:)   !! est rempli dans advtrac.h
      PRINT*,'check: tracernames ', tnom
     !!!new physics
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      emis(:)=wemis(:)
      PRINT*,'check: emis ',emis(1),emis(ngridmx)
      q2(:,:)=wq2(:,:)
      PRINT*,'check: q2 ',q2(1,1),q2(ngridmx,nlayermx+1)
      qsurf(:,:)=wqsurf(:,:)
      PRINT*,'check: qsurf ',qsurf(1,1),qsurf(ngridmx,nqmx)
      co2ice(:)=wco2ice(:)
      PRINT*,'check: co2 ',co2ice(1),co2ice(ngridmx)
      day_ini=wday_ini

c       artificially filling dyn3d/control.h is also required
c       > iphysiq is put in WRF to be set easily (cf ptimestep)
c       > day_step is simply deduced:
c
      day_step=daysec/ptimestep
      PRINT*,'Call to LMD physics:',day_step,' per Martian day'
c
      iphysiq=ptimestep
c
      !DO iq=1, nq 
      !  PRINT*, tnom(iq), pq(:,:,iq)
      !ENDDO

c
ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc