<HTML> <BODY BGCOLOR=#bbeeee LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE><a name='1'>
<font color=#447700>!#define NO_RESTRICT_ACCEL<a name='2'></font>
<font color=#447700>!#define NO_GFDLETAINIT<a name='3'></font>
<font color=#447700>!#define NO_UPSTREAM_ADVECTION<a name='4'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='5'></font>
<font color=#447700>!<a name='6'></font>
<A NAME='START_DOMAIN_NMM'><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A><a name='7'>
<font color=#993300>SUBROUTINE </font><font color=#cc0000>START_DOMAIN_NMM</font>(GRID, allowed_to_read , & <A href='../../call_to/START_DOMAIN_NMM.html' TARGET='index'>1</A>,<A href='../../call_from/START_DOMAIN_NMM.html' TARGET='index'>36</A><a name='8'>
<font color=#447700>!<a name='9'></font>
#include <<A href='../../html_code/include/nmm_dummy_args.inc.html'>nmm_dummy_args.inc</A>><A NAME="nmm_dummy_args.inc_1"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='10'>
<font color=#447700>!<a name='11'></font>
& )<a name='12'>
<font color=#447700>!----------------------------------------------------------------------<a name='13'></font>
<font color=#447700>!<a name='14'></font>
USE <A href='../../html_code/frame/module_domain.F.html#MODULE_DOMAIN'>MODULE_DOMAIN</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_DOMAIN_24"><a name='15'>
USE <A href='../../html_code/frame/module_driver_constants.F.html#MODULE_DRIVER_CONSTANTS'>MODULE_DRIVER_CONSTANTS</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_DRIVER_CONSTANTS_4"><a name='16'>
USE <A href='../../html_code/share/module_model_constants.F.html#MODULE_MODEL_CONSTANTS'>module_model_constants</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_MODEL_CONSTANTS_23"><a name='17'>
USE <A href='../../html_code/frame/module_configure.F.html#MODULE_CONFIGURE'>MODULE_CONFIGURE</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_CONFIGURE_25"><a name='18'>
USE <A href='../../html_code/frame/module_wrf_error.F.html#MODULE_WRF_ERROR'>MODULE_WRF_ERROR</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_WRF_ERROR_5"><a name='19'>
USE <A href='../../html_code/share/module_MPP.F.html#MODULE_MPP'>MODULE_MPP</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_MPP_4"><a name='20'>
USE <A href='../../html_code/dyn_nmm/module_CTLBLK.F.html#MODULE_CTLBLK'>MODULE_CTLBLK</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_CTLBLK_2"><a name='21'>
USE <A href='../../html_code/frame/module_dm_stubs.F.html#MODULE_DM'>MODULE_DM</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_DM_19"><a name='22'>
<font color=#447700>!<a name='23'></font>
USE <A href='../../html_code/dyn_nmm/module_IGWAVE_ADJUST.F.html#MODULE_IGWAVE_ADJUST'>MODULE_IGWAVE_ADJUST</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_IGWAVE_ADJUST_2">,ONLY: PDTE, PFDHT, DDAMP<a name='24'>
USE <A href='../../html_code/dyn_nmm/module_ADVECTION.F.html#MODULE_ADVECTION'>MODULE_ADVECTION</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_ADVECTION_2">, ONLY: ADVE, VAD2, HAD2<a name='25'>
USE <A href='../../html_code/dyn_nmm/module_NONHY_DYNAM.F.html#MODULE_NONHY_DYNAM'>MODULE_NONHY_DYNAM</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_NONHY_DYNAM_2">, ONLY: VADZ, HADZ<a name='26'>
USE <A href='../../html_code/dyn_nmm/module_DIFFUSION_NMM.F.html#MODULE_DIFFUSION_NMM'>MODULE_DIFFUSION_NMM</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_DIFFUSION_NMM_2">,ONLY: HDIFF<a name='27'>
USE <A href='../../html_code/dyn_nmm/module_BNDRY_COND.F.html#MODULE_BNDRY_COND'>MODULE_BNDRY_COND</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_BNDRY_COND_2">, ONLY: BOCOH, BOCOV<a name='28'>
USE <A href='../../html_code/phys/module_physics_init.F.html#MODULE_PHYSICS_INIT'>MODULE_PHYSICS_INIT</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_PHYSICS_INIT_2"><a name='29'>
USE <A href='../../html_code/phys/module_ra_gfdleta.F.html#MODULE_RA_GFDLETA'>MODULE_RA_GFDLETA</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODULE_RA_GFDLETA_3"><a name='30'>
<font color=#447700>!<a name='31'></font>
USE MODULE_EXT_INTERNAL<a name='32'>
<font color=#447700>!<a name='33'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='34'></font>
<font color=#447700>!<a name='35'></font>
IMPLICIT NONE<a name='36'>
<font color=#447700>!<a name='37'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='38'></font>
<font color=#447700>!***<a name='39'></font>
<font color=#447700>!*** Arguments<a name='40'></font>
<font color=#447700>!***<a name='41'></font>
TYPE(DOMAIN),INTENT(INOUT) :: GRID<a name='42'>
LOGICAL , INTENT(IN) :: allowed_to_read<a name='43'>
<font color=#447700>!<a name='44'></font>
#include <<A href='../../html_code/include/nmm_dummy_decl.inc.html'>nmm_dummy_decl.inc</A>><A NAME="nmm_dummy_decl.inc_2"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='45'>
<font color=#447700>!<a name='46'></font>
TYPE(GRID_CONFIG_REC_TYPE) :: CONFIG_FLAGS<a name='47'>
<font color=#447700>!<a name='48'></font>
<font color=#447700>!<a name='49'></font>
<font color=#447700>!***<a name='50'></font>
<font color=#447700>!*** LOCAL DATA<a name='51'></font>
<font color=#447700>!***<a name='52'></font>
INTEGER :: IDS,IDE,JDS,JDE,KDS,KDE &<a name='53'>
& ,IMS,IME,JMS,JME,KMS,KME &<a name='54'>
& ,IPS,IPE,JPS,JPE,KPS,KPE<a name='55'>
<font color=#447700>!<a name='56'></font>
INTEGER :: ERROR,LOOP<a name='57'>
<a name='58'>
REAL,ALLOCATABLE,DIMENSION(:) :: PHALF<a name='59'>
<font color=#447700>!<a name='60'></font>
REAL :: EPSB=0.1,EPSIN=9.8<a name='61'>
<font color=#447700>!<a name='62'></font>
INTEGER :: JHL=7<a name='63'>
<font color=#447700>!<a name='64'></font>
INTEGER :: I,IEND,IER,IERR,IFE,IFS,IHH,IHL,IHRSTB,II,IRTN &<a name='65'>
& ,ISIZ1,ISIZ2,ISTART,IX,J,J00,JFE,JFS,JHH,JJ &<a name='66'>
& ,JM1,JM2,JM3,JP1,JP2,JP3,JX &<a name='67'>
& ,K,K400,KBI,KBI2,KCCO2,KNT,KNTI,KOFF,KOFV &<a name='68'>
& ,LB,LLMH,LMHK,LMVK,LRECBC &<a name='69'>
& ,N,NMAP,NRADLH,NRADSH,NREC,NS,RECL,STAT &<a name='70'>
& ,STEPBL,STEPCU,STEPRA<a name='71'>
<font color=#447700>!<a name='72'></font>
INTEGER :: ILPAD2,IRPAD2,JBPAD2,JTPAD2<a name='73'>
INTEGER :: ITS,ITE,JTS,JTE,KTS,KTE<a name='74'>
<font color=#447700>!<a name='75'></font>
INTEGER,DIMENSION(3) :: LPTOP<a name='76'>
<font color=#447700>!<a name='77'></font>
REAL :: ADDL,APELM,APELMNW,APEM1,CAPA,CLOGES,DPLM,DZLM,EPS,ESE &<a name='78'>
& ,FAC1,FAC2,PDIF,PLM,PM1,PSFCK,PSS,PSUM,QLM,RANG &<a name='79'>
& ,SLPM,TERM1,THLM,TIME,TLM,TSFCK,ULM,VLM<a name='80'>
<font color=#447700>!<a name='81'></font>
<font color=#447700>!!! REAL :: BLDT,CWML,EXNSFC,G_INV,PLYR,PSFC,ROG,SFCZ,THSIJ,TL<a name='82'></font>
REAL :: CWML,EXNSFC,G_INV,PLYR,PSURF,ROG,SFCZ,THSIJ,TL<a name='83'>
REAL :: TSTART,TEND,TPREC,THEAT,TCLOD,TRDSW,TRDLW,TSRFC<a name='84'>
<a name='85'>
<font color=#447700>!<a name='86'></font>
<font color=#447700>!!! REAL,ALLOCATABLE,DIMENSION(:,:) :: RAINBL,RAINNC,RAINNC &<a name='87'></font>
INTEGER,ALLOCATABLE,DIMENSION(:,:) :: LOWLYR<a name='88'>
REAL,ALLOCATABLE,DIMENSION(:) :: SFULL,SMID<a name='89'>
<font color=#447700>!state real DZS l dyn_em - Z ir <a name='90'></font>
<font color=#447700>!state real CLDFRA ikj dyn_em 1 - r<a name='91'></font>
<font color=#447700>!state real RQCBLTEN ikj dyn_em 1 - r<a name='92'></font>
<font color=#447700>!state real RQIBLTEN ikj dyn_em 1 - r<a name='93'></font>
<font color=#447700>!state real RQVBLTEN ikj dyn_em 1 - r<a name='94'></font>
<font color=#447700>!state real RTHBLTEN ikj dyn_em 1 - r<a name='95'></font>
<font color=#447700>!state real RUBLTEN ikj dyn_em 1 - r<a name='96'></font>
<font color=#447700>!state real RVBLTEN ikj dyn_em 1 - r<a name='97'></font>
<font color=#447700>!state real RQCCUTEN ikj dyn_em 1 - r<a name='98'></font>
<font color=#447700>!state real RQICUTEN ikj dyn_em 1 - r<a name='99'></font>
<font color=#447700>!state real RQRCUTEN ikj dyn_em 1 - r<a name='100'></font>
<font color=#447700>!state real RQSCUTEN ikj dyn_em 1 - r<a name='101'></font>
<font color=#447700>!state real RQVCUTEN ikj dyn_em 1 - r<a name='102'></font>
<font color=#447700>!state real RTHCUTEN ikj dyn_em 1 - r<a name='103'></font>
<font color=#447700>!state real RTHRATEN ikj dyn_em 1 - r<a name='104'></font>
<font color=#447700>!state real RTHRATENLW ikj dyn_em 1 - r<a name='105'></font>
<font color=#447700>!state real RTHRATENSW ikj dyn_em 1 - r<a name='106'></font>
<font color=#447700>!state real TSLB ilj dyn_em 1 Z irh<a name='107'></font>
<font color=#447700>!state real ZS l dyn_em - Z ir <a name='108'></font>
REAL,ALLOCATABLE,DIMENSION(:) :: DZS,ZS<a name='109'>
REAL,ALLOCATABLE,DIMENSION(:,:,:) :: RQCBLTEN,RQIBLTEN &<a name='110'>
& ,RQVBLTEN,RTHBLTEN &<a name='111'>
& ,RUBLTEN,RVBLTEN &<a name='112'>
& ,RQCCUTEN,RQICUTEN,RQRCUTEN &<a name='113'>
& ,RQSCUTEN,RQVCUTEN,RTHCUTEN &<a name='114'>
& ,RTHRATEN &<a name='115'>
& ,RTHRATENLW,RTHRATENSW &<a name='116'>
& ,TSLB<a name='117'>
REAL,ALLOCATABLE,DIMENSION(:,:) :: EMISS,GLW,GSW,HFX &<a name='118'>
& ,MAVAIL,NCA &<a name='119'>
& ,QFX,RAINBL,RAINC,RAINNC &<a name='120'>
& ,RAINCV,RAINNCV &<a name='121'>
& ,SNOWC,THC,TMN,TSFC &<a name='122'>
& ,XLAND,XLAT,XLONG<a name='123'>
<a name='124'>
REAL,ALLOCATABLE,DIMENSION(:,:) :: Z0_DUM<a name='125'>
<font color=#447700>!<a name='126'></font>
<font color=#447700>!!! REAL,ALLOCATABLE,DIMENSION(:,:,:) :: W0AVG,ZINT,ZMID<a name='127'></font>
#if 0<a name='128'>
REAL,ALLOCATABLE,DIMENSION(:,:,:) :: W0AVG<a name='129'>
#endif<a name='130'>
LOGICAL :: E_BDY,N_BDY,S_BDY,W_BDY,WARM_RAIN<a name='131'>
integer :: jam,retval<a name='132'>
character(20) :: seeout="hi08.t00z.nhbmeso"<a name='133'>
real :: dummyx(791)<a name='134'>
integer myproc<a name='135'>
real :: dsig,dsigsum,pdbot,pdtot,rpdtot<a name='136'>
real :: fisx,ht,prodx,rg<a name='137'>
integer :: i_t=096,j_t=195,n_t=11<a name='138'>
integer :: i_u=49,j_u=475,n_u=07<a name='139'>
integer :: i_v=49,j_v=475,n_v=07<a name='140'>
<a name='141'>
#ifdef DEREF_KLUDGE<a name='142'>
<font color=#447700>! see http://www.mmm.ucar.edu/wrf/WG2/topics/deref_kludge.htm<a name='143'></font>
INTEGER :: sm31 , em31 , sm32 , em32 , sm33 , em33<a name='144'>
INTEGER :: sm31x, em31x, sm32x, em32x, sm33x, em33x<a name='145'>
INTEGER :: sm31y, em31y, sm32y, em32y, sm33y, em33y<a name='146'>
#endif<a name='147'>
<a name='148'>
#include "<A href='../../html_code/include/deref_kludge.h.html'>deref_kludge.h</A>"<A NAME="deref_kludge.h_3"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='149'>
<a name='150'>
<font color=#447700>!<a name='151'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='152'></font>
#define COPY_IN<a name='153'>
#include <<A href='../../html_code/include/nmm_scalar_derefs.inc.html'>nmm_scalar_derefs.inc</A>><A NAME="nmm_scalar_derefs.inc_4"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='154'>
#ifdef DM_PARALLEL<a name='155'>
# include <<A href='../../html_code/include/nmm_data_calls.inc.html'>nmm_data_calls.inc</A>><A NAME="nmm_data_calls.inc_5"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='156'>
#endif<a name='157'>
<font color=#447700>!----------------------------------------------------------------------<a name='158'></font>
<font color=#447700>!**********************************************************************<a name='159'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='160'></font>
<font color=#447700>!<a name='161'></font>
CALL <A href='../../html_code/frame/module_domain.F.html#GET_IJK_FROM_GRID'>GET_IJK_FROM_GRID</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="GET_IJK_FROM_GRID_11">(GRID, &<a name='162'>
& IDS,IDE,JDS,JDE,KDS,KDE, &<a name='163'>
& IMS,IME,JMS,JME,KMS,KME, &<a name='164'>
& IPS,IPE,JPS,JPE,KPS,KPE)<a name='165'>
<font color=#447700>!<a name='166'></font>
ITS=IPS<a name='167'>
ITE=IPE<a name='168'>
JTS=JPS<a name='169'>
JTE=JPE<a name='170'>
KTS=KPS<a name='171'>
KTE=KPE<a name='172'>
<a name='173'>
CALL <A href='../../html_code/frame/module_configure.F.html#MODEL_TO_GRID_CONFIG_REC'>model_to_grid_config_rec</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="MODEL_TO_GRID_CONFIG_REC_11">(grid%id,model_config_rec &<a name='174'>
& ,config_flags)<a name='175'>
<font color=#447700>!<a name='176'></font>
RESTRT=config_flags%restart<a name='177'>
write(0,*) 'set RESTRT to: ', RESTRT<a name='178'>
<a name='179'>
#if 1<a name='180'>
IF(IME.GT. NMM_MAX_DIM )THEN<a name='181'>
WRITE(wrf_err_message,*) &<a name='182'>
'start_domain_nmm ime (',ime,') > ',NMM_MAX_DIM, &<a name='183'>
'. Increase NMM_MAX_DIM in configure.wrf, clean, and recompile.'<a name='184'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_ERROR_FATAL'>WRF_ERROR_FATAL</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_ERROR_FATAL_54">(wrf_err_message)<a name='185'>
ENDIF<a name='186'>
<font color=#447700>!<a name='187'></font>
IF(JME.GT. NMM_MAX_DIM )THEN<a name='188'>
WRITE(wrf_err_message,*) &<a name='189'>
'start_domain_nmm jme (',jme,') > ',NMM_MAX_DIM, &<a name='190'>
'. Increase NMM_MAX_DIM in configure.wrf, clean, and recompile.'<a name='191'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_ERROR_FATAL'>WRF_ERROR_FATAL</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_ERROR_FATAL_55">(wrf_err_message)<a name='192'>
ENDIF<a name='193'>
#else<a name='194'>
IF(IMS.GT.-2.OR.IME.GT. NMM_MAX_DIM )THEN<a name='195'>
WRITE(wrf_err_message,*) &<a name='196'>
'start_domain_nmm ims(',ims,' > -2 or ime (',ime,') > ',NMM_MAX_DIM, &<a name='197'>
'. Increase NMM_MAX_DIM in configure.wrf, clean, and recompile.'<a name='198'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_ERROR_FATAL'>WRF_ERROR_FATAL</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_ERROR_FATAL_56">(wrf_err_message)<a name='199'>
ENDIF<a name='200'>
<font color=#447700>!<a name='201'></font>
IF(JMS.GT.-2.OR.JME.GT. NMM_MAX_DIM )THEN<a name='202'>
WRITE(wrf_err_message,*) &<a name='203'>
'start_domain_nmm jms(',jms,' > -2 or jme (',jme,') > ',NMM_MAX_DIM, &<a name='204'>
'. Increase NMM_MAX_DIM in configure.wrf, clean, and recompile.'<a name='205'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_ERROR_FATAL'>WRF_ERROR_FATAL</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_ERROR_FATAL_57">(wrf_err_message)<a name='206'>
ENDIF<a name='207'>
#endif<a name='208'>
<font color=#447700>!<a name='209'></font>
<font color=#447700>!---------------------------------------------------------------------- <a name='210'></font>
<font color=#447700>!<a name='211'></font>
WRITE(0,196)IHRST,IDAT<a name='212'>
WRITE(LIST,196)IHRST,IDAT<a name='213'>
196 FORMAT(' FORECAST BEGINS ',I2,' GMT ',2(I2,'/'),I4)<a name='214'>
<font color=#447700>!!!!!!tlb<a name='215'></font>
<font color=#447700>!!!! For now, set NPES to 1<a name='216'></font>
NPES=1<a name='217'>
<font color=#447700>!!!!!!tlb<a name='218'></font>
MY_IS_GLB=IPS<a name='219'>
MY_IE_GLB=IPE-1<a name='220'>
MY_JS_GLB=JPS<a name='221'>
MY_JE_GLB=JPE-1<a name='222'>
<font color=#447700>!<a name='223'></font>
IM=IPE-1<a name='224'>
JM=JPE-1<a name='225'>
<font color=#447700>!!!!!!!!!<a name='226'></font>
<font color=#447700>!! All "my" variables defined below have had the IDE or JDE specification<a name='227'></font>
<font color=#447700>!! reduced by 1<a name='228'></font>
<font color=#447700>!!!!!!!!!!!<a name='229'></font>
<a name='230'>
MYIS=MAX(IDS,IPS)<a name='231'>
MYIE=MIN(IDE-1,IPE)<a name='232'>
MYJS=MAX(JDS,JPS)<a name='233'>
MYJE=MIN(JDE-1,JPE)<a name='234'>
<a name='235'>
MYIS1 =MAX(IDS+1,IPS)<a name='236'>
MYIE1 =MIN(IDE-2,IPE)<a name='237'>
MYJS2 =MAX(JDS+2,JPS)<a name='238'>
MYJE2 =MIN(JDE-3,JPE)<a name='239'>
<font color=#447700>!<a name='240'></font>
MYIS_P1=MAX(IDS,IPS-1)<a name='241'>
MYIE_P1=MIN(IDE-1,IPE+1)<a name='242'>
MYIS_P2=MAX(IDS,IPS-2)<a name='243'>
MYIE_P2=MIN(IDE-1,IPE+2)<a name='244'>
MYIS_P3=MAX(IDS,IPS-3)<a name='245'>
MYIE_P3=MIN(IDE-1,IPE+3)<a name='246'>
MYJS_P3=MAX(JDS,JPS-3)<a name='247'>
MYJE_P3=MIN(JDE-1,JPE+3)<a name='248'>
MYIS_P4=MAX(IDS,IPS-4)<a name='249'>
MYIE_P4=MIN(IDE-1,IPE+4)<a name='250'>
MYJS_P4=MAX(JDS,JPS-4)<a name='251'>
MYJE_P4=MIN(JDE-1,JPE+4)<a name='252'>
MYIS_P5=MAX(IDS,IPS-5)<a name='253'>
MYIE_P5=MIN(IDE-1,IPE+5)<a name='254'>
MYJS_P5=MAX(JDS,JPS-5)<a name='255'>
MYJE_P5=MIN(JDE-1,JPE+5)<a name='256'>
<font color=#447700>!<a name='257'></font>
MYIS1_P1=MAX(IDS+1,IPS-1)<a name='258'>
MYIE1_P1=MIN(IDE-2,IPE+1)<a name='259'>
MYIS1_P2=MAX(IDS+1,IPS-2)<a name='260'>
MYIE1_P2=MIN(IDE-2,IPE+2)<a name='261'>
<font color=#447700>!<a name='262'></font>
MYJS1_P1=MAX(JDS+1,JPS-1)<a name='263'>
MYJS2_P1=MAX(JDS+2,JPS-1)<a name='264'>
MYJE1_P1=MIN(JDE-2,JPE+1)<a name='265'>
MYJE2_P1=MIN(JDE-3,JPE+1)<a name='266'>
MYJS1_P2=MAX(JDS+1,JPS-2)<a name='267'>
MYJE1_P2=MIN(JDE-2,JPE+2)<a name='268'>
MYJS2_P2=MAX(JDS+2,JPS-2)<a name='269'>
MYJE2_P2=MIN(JDE-3,JPE+2)<a name='270'>
MYJS1_P3=MAX(JDS+1,JPS-3)<a name='271'>
MYJE1_P3=MIN(JDE-2,JPE+3)<a name='272'>
MYJS2_P3=MAX(JDS+2,JPS-3)<a name='273'>
MYJE2_P3=MIN(JDE-3,JPE+3)<a name='274'>
<font color=#447700>!!!!!!!!!!!<a name='275'></font>
<font color=#447700>!<a name='276'></font>
#ifdef DM_PARALLEL<a name='277'>
<a name='278'>
call <A href='../../html_code/frame/module_dm_stubs.F.html#WRF_GET_MYPROC'>wrf_get_myproc</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_GET_MYPROC_9">(myproc)<a name='279'>
<a name='280'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_1.inc.html'>HALO_NMM_INIT_1.inc</A>><A NAME="HALO_NMM_INIT_1.inc_6"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='281'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_2.inc.html'>HALO_NMM_INIT_2.inc</A>><A NAME="HALO_NMM_INIT_2.inc_7"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='282'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_3.inc.html'>HALO_NMM_INIT_3.inc</A>><A NAME="HALO_NMM_INIT_3.inc_8"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='283'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_4.inc.html'>HALO_NMM_INIT_4.inc</A>><A NAME="HALO_NMM_INIT_4.inc_9"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='284'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_5.inc.html'>HALO_NMM_INIT_5.inc</A>><A NAME="HALO_NMM_INIT_5.inc_10"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='285'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_6.inc.html'>HALO_NMM_INIT_6.inc</A>><A NAME="HALO_NMM_INIT_6.inc_11"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='286'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_7.inc.html'>HALO_NMM_INIT_7.inc</A>><A NAME="HALO_NMM_INIT_7.inc_12"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='287'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_8.inc.html'>HALO_NMM_INIT_8.inc</A>><A NAME="HALO_NMM_INIT_8.inc_13"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='288'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_9.inc.html'>HALO_NMM_INIT_9.inc</A>><A NAME="HALO_NMM_INIT_9.inc_14"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='289'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_10.inc.html'>HALO_NMM_INIT_10.inc</A>><A NAME="HALO_NMM_INIT_10.inc_15"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='290'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_11.inc.html'>HALO_NMM_INIT_11.inc</A>><A NAME="HALO_NMM_INIT_11.inc_16"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='291'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_12.inc.html'>HALO_NMM_INIT_12.inc</A>><A NAME="HALO_NMM_INIT_12.inc_17"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='292'>
<a name='293'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_13.inc.html'>HALO_NMM_INIT_13.inc</A>><A NAME="HALO_NMM_INIT_13.inc_18"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='294'>
<a name='295'>
<font color=#447700>! CALL wrf_shutdown<a name='296'></font>
<font color=#447700>! stop<a name='297'></font>
<a name='298'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_14.inc.html'>HALO_NMM_INIT_14.inc</A>><A NAME="HALO_NMM_INIT_14.inc_19"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='299'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_15.inc.html'>HALO_NMM_INIT_15.inc</A>><A NAME="HALO_NMM_INIT_15.inc_20"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='300'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_16.inc.html'>HALO_NMM_INIT_16.inc</A>><A NAME="HALO_NMM_INIT_16.inc_21"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='301'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_17.inc.html'>HALO_NMM_INIT_17.inc</A>><A NAME="HALO_NMM_INIT_17.inc_22"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='302'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_18.inc.html'>HALO_NMM_INIT_18.inc</A>><A NAME="HALO_NMM_INIT_18.inc_23"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='303'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_19.inc.html'>HALO_NMM_INIT_19.inc</A>><A NAME="HALO_NMM_INIT_19.inc_24"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='304'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_20.inc.html'>HALO_NMM_INIT_20.inc</A>><A NAME="HALO_NMM_INIT_20.inc_25"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='305'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_21.inc.html'>HALO_NMM_INIT_21.inc</A>><A NAME="HALO_NMM_INIT_21.inc_26"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='306'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_22.inc.html'>HALO_NMM_INIT_22.inc</A>><A NAME="HALO_NMM_INIT_22.inc_27"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='307'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_23.inc.html'>HALO_NMM_INIT_23.inc</A>><A NAME="HALO_NMM_INIT_23.inc_28"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='308'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_24.inc.html'>HALO_NMM_INIT_24.inc</A>><A NAME="HALO_NMM_INIT_24.inc_29"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='309'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_25.inc.html'>HALO_NMM_INIT_25.inc</A>><A NAME="HALO_NMM_INIT_25.inc_30"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='310'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_26.inc.html'>HALO_NMM_INIT_26.inc</A>><A NAME="HALO_NMM_INIT_26.inc_31"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='311'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_27.inc.html'>HALO_NMM_INIT_27.inc</A>><A NAME="HALO_NMM_INIT_27.inc_32"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='312'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_28.inc.html'>HALO_NMM_INIT_28.inc</A>><A NAME="HALO_NMM_INIT_28.inc_33"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='313'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_29.inc.html'>HALO_NMM_INIT_29.inc</A>><A NAME="HALO_NMM_INIT_29.inc_34"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='314'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_30.inc.html'>HALO_NMM_INIT_30.inc</A>><A NAME="HALO_NMM_INIT_30.inc_35"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='315'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_31.inc.html'>HALO_NMM_INIT_31.inc</A>><A NAME="HALO_NMM_INIT_31.inc_36"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='316'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_32.inc.html'>HALO_NMM_INIT_32.inc</A>><A NAME="HALO_NMM_INIT_32.inc_37"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='317'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_33.inc.html'>HALO_NMM_INIT_33.inc</A>><A NAME="HALO_NMM_INIT_33.inc_38"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='318'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_34.inc.html'>HALO_NMM_INIT_34.inc</A>><A NAME="HALO_NMM_INIT_34.inc_39"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='319'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_35.inc.html'>HALO_NMM_INIT_35.inc</A>><A NAME="HALO_NMM_INIT_35.inc_40"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='320'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_36.inc.html'>HALO_NMM_INIT_36.inc</A>><A NAME="HALO_NMM_INIT_36.inc_41"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='321'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_37.inc.html'>HALO_NMM_INIT_37.inc</A>><A NAME="HALO_NMM_INIT_37.inc_42"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='322'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_38.inc.html'>HALO_NMM_INIT_38.inc</A>><A NAME="HALO_NMM_INIT_38.inc_43"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='323'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_39.inc.html'>HALO_NMM_INIT_39.inc</A>><A NAME="HALO_NMM_INIT_39.inc_44"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='324'>
#endif<a name='325'>
<a name='326'>
DO J=MYJS_P4,MYJE_P4<a name='327'>
IHEG(J)=MOD(J+1,2)<a name='328'>
IHWG(J)=IHEG(J)-1<a name='329'>
IVEG(J)=MOD(J,2)<a name='330'>
IVWG(J)=IVEG(J)-1<a name='331'>
ENDDO<a name='332'>
<font color=#447700>!<a name='333'></font>
DO J=MYJS_P4,MYJE_P4<a name='334'>
IVW(J)=IVWG(J)<a name='335'>
IVE(J)=IVEG(J)<a name='336'>
IHE(J)=IHEG(J)<a name='337'>
IHW(J)=IHWG(J)<a name='338'>
ENDDO<a name='339'>
<font color=#447700>!<a name='340'></font>
CAPA=R_D/CP<a name='341'>
LM=KPE-KPS+1<a name='342'>
<font color=#447700>!<a name='343'></font>
IFS=IPS<a name='344'>
JFS=JPS<a name='345'>
JFE=MIN(JPE,JDE-1)<a name='346'>
IFE=MIN(IPE,IDE-1)<a name='347'>
<font color=#447700>!<a name='348'></font>
IF(.NOT.RESTRT)THEN<a name='349'>
call <A href='../../html_code/frame/module_dm_stubs.F.html#WRF_GET_MYPROC'>wrf_get_myproc</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_GET_MYPROC_10">(mype)<a name='350'>
DO J=JFS,JFE<a name='351'>
DO I=IFS,IFE<a name='352'>
LLMH=LMH(I,J)<a name='353'>
KOFF=KPE-1-LLMH<a name='354'>
PDSL(I,J) =PD(I,J)*RES(I,J)<a name='355'>
PREC(I,J) =0.<a name='356'>
ACPREC(I,J)=0.<a name='357'>
CUPREC(I,J)=0.<a name='358'>
rg=1./g<a name='359'>
ht=fis(i,j)*rg<a name='360'>
<font color=#447700>!!! fisx=ht*g<a name='361'></font>
<font color=#447700>! fisx=max(fis(i,j),0.)<a name='362'></font>
<font color=#447700>! prodx=Z0(I,J)*Z0MAX<a name='363'></font>
<font color=#447700>! Z0(I,J) =SM(I,J)*Z0SEA+(1.-SM(I,J))* &<a name='364'></font>
<font color=#447700>! & (Z0(I,J)*Z0MAX+FISx *FCM+Z0LAND)<a name='365'></font>
<font color=#447700>!!! & (prodx +FISx *FCM+Z0LAND)<a name='366'></font>
QS(I,J) =0.<a name='367'>
AKMS(I,J) =0.<a name='368'>
AKHS(I,J) =0.<a name='369'>
TWBS(I,J) =0.<a name='370'>
QWBS(I,J) =0.<a name='371'>
CLDEFI(I,J)=1.<a name='372'>
<font color=#447700>!!!! HTOP(I,J) =REAL(LLMH)<a name='373'></font>
<font color=#447700>!!!! HBOT(I,J) =REAL(LLMH)<a name='374'></font>
HTOP(I,J) =REAL(KTS)<a name='375'>
HTOPD(I,J) =REAL(KTS)<a name='376'>
HTOPS(I,J) =REAL(KTS)<a name='377'>
HBOT(I,J) =REAL(KTE)<a name='378'>
HBOTD(I,J) =REAL(KTE)<a name='379'>
HBOTS(I,J) =REAL(KTE)<a name='380'>
<font color=#447700>!***<a name='381'></font>
<font color=#447700>!*** AT THIS POINT, WE MUST CALCULATE THE INITIAL POTENTIAL TEMPERATURE<a name='382'></font>
<font color=#447700>!*** OF THE SURFACE AND OF THE SUBGROUND.<a name='383'></font>
<font color=#447700>!*** EXTRAPOLATE DOWN FOR INITIAL SURFACE POTENTIAL TEMPERATURE.<a name='384'></font>
<font color=#447700>!*** ALSO DO THE SHELTER PRESSURE.<a name='385'></font>
<font color=#447700>!***<a name='386'></font>
PM1=AETA1(KOFF+1)*PDTOP+AETA2(KOFF+1)*PDSL(I,J)+PT<a name='387'>
APEM1=(1.E5/PM1)**CAPA<a name='388'>
<a name='389'>
IF (NMM_TSK(I,J) .ge. 200.) THEN <font color=#447700>! have a specific skin temp, use it<a name='390'></font>
THS(I,J)=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='391'>
TSFCK=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))<a name='392'>
ELSE <font color=#447700>! use lowest layer as a proxy<a name='393'></font>
THS(I,J)=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='394'>
TSFCK=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))<a name='395'>
ENDIF<a name='396'>
<a name='397'>
if (I .eq. IFE/2 .and. J .eq. JFE/2) then<a name='398'>
write(6,*) 'I,J,T(I,KOFF+1,J),NMM_TSK(I,J):: ', I,J,T(I,KOFF+1,J),NMM_TSK(I,J)<a name='399'>
write(6,*) 'THS(I,J): ', THS(I,J)<a name='400'>
endif<a name='401'>
<a name='402'>
PSFCK=PD(I,J)+PDTOP+PT<a name='403'>
<font color=#447700>!<a name='404'></font>
IF(SM(I,J).LT.0.5) THEN<a name='405'>
QS(I,J)=PQ0/PSFCK*EXP(A2*(TSFCK-A3)/(TSFCK-A4))<a name='406'>
ELSEIF(SM(I,J).GT.0.5) THEN<a name='407'>
THS(I,J)=SST(I,J)*(1.E5/(PD(I,J)+PDTOP+PT))**CAPA<a name='408'>
ENDIF<a name='409'>
<font color=#447700>!<a name='410'></font>
TERM1=-0.068283/T(I,KOFF+1,J)<a name='411'>
PSHLTR(I,J)=(PD(I,J)+PDTOP+PT)*EXP(TERM1)<a name='412'>
<font color=#447700>!<a name='413'></font>
USTAR(I,J)=0.1<a name='414'>
THZ0(I,J)=THS(I,J)<a name='415'>
QZ0(I,J)=QS(I,J)<a name='416'>
UZ0(I,J)=0.<a name='417'>
VZ0(I,J)=0.<a name='418'>
<font color=#447700>! <a name='419'></font>
ENDDO<a name='420'>
ENDDO<a name='421'>
<a name='422'>
<font color=#447700>!***<a name='423'></font>
<font color=#447700>!*** INITIALIZE 3D MASKS<a name='424'></font>
<font color=#447700>!***<a name='425'></font>
DO J=JFS,JFE<a name='426'>
DO K=KPS,KPE<a name='427'>
DO I=IFS,IFE<a name='428'>
HTM(I,K,J)=1.<a name='429'>
VTM(I,K,J)=1.<a name='430'>
ENDDO<a name='431'>
ENDDO<a name='432'>
ENDDO<a name='433'>
<font color=#447700>!***<a name='434'></font>
<font color=#447700>!*** INITIALIZE CLOUD FIELDS<a name='435'></font>
<font color=#447700>!***<a name='436'></font>
DO J=JFS,JFE<a name='437'>
DO K=KPS,KPE<a name='438'>
DO I=IFS,IFE<a name='439'>
CWM(I,K,J)=0.<a name='440'>
ENDDO<a name='441'>
ENDDO<a name='442'>
ENDDO<a name='443'>
<font color=#447700>!***<a name='444'></font>
<font color=#447700>!*** INITIALIZE ACCUMULATOR ARRAYS TO ZERO.<a name='445'></font>
<font color=#447700>!***<a name='446'></font>
ARDSW=0.0<a name='447'>
ARDLW=0.0<a name='448'>
ASRFC=0.0<a name='449'>
AVRAIN=0.0<a name='450'>
AVCNVC=0.0<a name='451'>
<font color=#447700>!<a name='452'></font>
DO J=JFS,JFE<a name='453'>
DO I=IFS,IFE<a name='454'>
ACFRCV(I,J)=0.<a name='455'>
NCFRCV(I,J)=0<a name='456'>
ACFRST(I,J)=0.<a name='457'>
NCFRST(I,J)=0<a name='458'>
ACSNOW(I,J)=0.<a name='459'>
ACSNOM(I,J)=0.<a name='460'>
SSROFF(I,J)=0.<a name='461'>
BGROFF(I,J)=0.<a name='462'>
ALWIN(I,J) =0.<a name='463'>
ALWOUT(I,J)=0.<a name='464'>
ALWTOA(I,J)=0.<a name='465'>
ASWIN(I,J) =0.<a name='466'>
ASWOUT(I,J)=0.<a name='467'>
ASWTOA(I,J)=0.<a name='468'>
SFCSHX(I,J)=0.<a name='469'>
SFCLHX(I,J)=0.<a name='470'>
SUBSHX(I,J)=0.<a name='471'>
SNOPCX(I,J)=0.<a name='472'>
SFCUVX(I,J)=0.<a name='473'>
SFCEVP(I,J)=0.<a name='474'>
POTEVP(I,J)=0.<a name='475'>
POTFLX(I,J)=0.<a name='476'>
ENDDO<a name='477'>
ENDDO<a name='478'>
<font color=#447700>!***<a name='479'></font>
<font color=#447700>!*** INITIALIZE SATURATION SPECIFIC HUMIDITY OVER THE WATER.<a name='480'></font>
<font color=#447700>!***<a name='481'></font>
EPS=R_D/R_V<a name='482'>
<font color=#447700>!<a name='483'></font>
DO J=JFS,JFE<a name='484'>
DO I=IFS,IFE<a name='485'>
IF(SM(I,J).GT.0.5)THEN<a name='486'>
CLOGES =-CM1/SST(I,J)-CM2*ALOG10(SST(I,J))+CM3<a name='487'>
ESE = 10.**(CLOGES+2.)<a name='488'>
QS(I,J)= SM(I,J)*EPS*ESE/(PD(I,J)+PDTOP+PT-ESE*(1.-EPS))<a name='489'>
ENDIF<a name='490'>
ENDDO<a name='491'>
ENDDO<a name='492'>
<font color=#447700>!*** <a name='493'></font>
<font color=#447700>!*** INITIALIZE TURBULENT KINETIC ENERGY (TKE) TO A SMALL<a name='494'></font>
<font color=#447700>!*** VALUE (EPSQ2) ABOVE GROUND. SET TKE TO ZERO IN THE<a name='495'></font>
<font color=#447700>!*** THE LOWEST MODEL LAYER. IN THE LOWEST TWO ATMOSPHERIC<a name='496'></font>
<font color=#447700>!*** ETA LAYERS SET TKE TO A SMALL VALUE (Q2INI).<a name='497'></font>
<font color=#447700>!***<a name='498'></font>
DO J=JFS,JFE<a name='499'>
DO K=KPS,KPE-1<a name='500'>
DO I=IFS,IFE<a name='501'>
Q2(I,K,J)=HTM(I,K+1,J)*HBM2(I,J)*EPSQ2<a name='502'>
ENDDO<a name='503'>
ENDDO<a name='504'>
ENDDO<a name='505'>
<font color=#447700>!<a name='506'></font>
DO J=JFS,JFE<a name='507'>
DO I=IFS,IFE<a name='508'>
Q2(I,LM,J) = 0.<a name='509'>
LLMH = LMH(I,J)<a name='510'>
Q2(I,LLMH-2,J)= HBM2(I,J)*Q2INI<a name='511'>
Q2(I,LLMH-1,J)= HBM2(I,J)*Q2INI<a name='512'>
ENDDO<a name='513'>
ENDDO<a name='514'>
<font color=#447700>!*** <a name='515'></font>
<font color=#447700>!*** PAD ABOVE GROUND SPECIFIC HUMIDITY IF IT IS TOO SMALL.<a name='516'></font>
<font color=#447700>!*** INITIALIZE LATENT HEATING ACCUMULATION ARRAYS.<a name='517'></font>
<font color=#447700>!***<a name='518'></font>
DO J=JFS,JFE<a name='519'>
DO K=KPS,KPE<a name='520'>
DO I=IFS,IFE<a name='521'>
IF(Q(I,K,J).LT.EPSQ)Q(I,K,J)=EPSQ*HTM(I,K,J)<a name='522'>
TRAIN(I,K,J)=0.<a name='523'>
TCUCN(I,K,J)=0.<a name='524'>
ENDDO<a name='525'>
ENDDO<a name='526'>
ENDDO<a name='527'>
<font color=#447700>!<a name='528'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='529'></font>
<font color=#447700>!*** END OF SCRATCH START INITIALIZATION BLOCK.<a name='530'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='531'></font>
<font color=#447700>!<a name='532'></font>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_14">('INIT: INITIALIZED ARRAYS FOR CLEAN START')<a name='533'>
ENDIF <font color=#447700>! <--- (not restart)<a name='534'></font>
<a name='535'>
<a name='536'>
<font color=#447700>!<a name='537'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='538'></font>
<font color=#447700>!*** INITIALIZE PHYSICS VARIABLES IF STARTING THIS RUN FROM SCRATCH.<a name='539'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='540'></font>
<font color=#447700>!<a name='541'></font>
IF(NEST)THEN<a name='542'>
DO J=JFS,JFE<a name='543'>
DO I=IFS,IFE<a name='544'>
<font color=#447700>!<a name='545'></font>
LLMH=LMH(I,J)<a name='546'>
KOFF=KPE-1-LLMH<a name='547'>
<font color=#447700>!<a name='548'></font>
IF(T(I,KOFF+1,J).EQ.0.)THEN<a name='549'>
T(I,KOFF+1,J)=T(I,KOFF+2,J)<a name='550'>
ENDIF<a name='551'>
<font color=#447700>!<a name='552'></font>
TERM1=-0.068283/T(I,KOFF+1,J)<a name='553'>
PSHLTR(I,J)=(PD(I,J)+PDTOP+PT)*EXP(TERM1)<a name='554'>
ENDDO<a name='555'>
ENDDO<a name='556'>
ENDIF<a name='557'>
<font color=#447700>!<a name='558'></font>
IF(.NOT.RESTRT)THEN<a name='559'>
DO J=JFS,JFE<a name='560'>
DO I=IFS,IFE<a name='561'>
LLMH=LMH(I,J)<a name='562'>
KOFF=KPE-1-LLMH<a name='563'>
PDSL(I,J) =PD(I,J)*RES(I,J)<a name='564'>
PREC(I,J) =0.<a name='565'>
ACPREC(I,J)=0.<a name='566'>
CUPREC(I,J)=0.<a name='567'>
<font color=#447700>! Z0(I,J) =SM(I,J)*Z0SEA+(1.-SM(I,J))* &<a name='568'></font>
<font color=#447700>! (FIS(I,J)*FCM+Z0LAND+Z0(I,J))<a name='569'></font>
QS(I,J) =0.<a name='570'>
AKMS(I,J) =0.<a name='571'>
AKHS(I,J) =0.<a name='572'>
TWBS(I,J) =0.<a name='573'>
QWBS(I,J) =0.<a name='574'>
CLDEFI(I,J)=1.<a name='575'>
<font color=#447700>!!!! HTOP(I,J) =REAL(LLMH)<a name='576'></font>
<font color=#447700>!!!! HBOT(I,J) =REAL(LLMH)<a name='577'></font>
HTOP(I,J) =REAL(KTS)<a name='578'>
HBOT(I,J) =REAL(KTE)<a name='579'>
<font color=#447700>!***<a name='580'></font>
<font color=#447700>!*** AT THIS POINT, WE MUST CALCULATE THE INITIAL POTENTIAL TEMPERATURE<a name='581'></font>
<font color=#447700>!*** OF THE SURFACE AND OF THE SUBGROUND.<a name='582'></font>
<font color=#447700>!*** EXTRAPOLATE DOWN FOR INITIAL SURFACE POTENTIAL TEMPERATURE.<a name='583'></font>
<font color=#447700>!*** ALSO DO THE SHELTER PRESSURE.<a name='584'></font>
<font color=#447700>!***<a name='585'></font>
PM1=AETA1(KOFF+1)*PDTOP+AETA2(KOFF+1)*PDSL(I,J)+PT<a name='586'>
APEM1=(1.E5/PM1)**CAPA<a name='587'>
IF (NMM_TSK(I,J) .ge. 200.) THEN <font color=#447700>! have a specific skin temp, use it<a name='588'></font>
THS(I,J)=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='589'>
TSFCK=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))<a name='590'>
ELSE <font color=#447700>! use lowest layer as a proxy<a name='591'></font>
THS(I,J)=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='592'>
TSFCK=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))<a name='593'>
ENDIF<a name='594'>
<font color=#447700>!!!!! THS(I,J)=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='595'></font>
<font color=#447700>! THS(I,J)=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))*APEM1<a name='596'></font>
<font color=#447700>!! TSFCK=T(I,KOFF+1,J)*(1.+P608*Q(I,KOFF+1,J))<a name='597'></font>
<font color=#447700>! TSFCK=NMM_TSK(I,J)*(1.+P608*Q(I,KOFF+1,J))<a name='598'></font>
<a name='599'>
PSFCK=PD(I,J)+PDTOP+PT<a name='600'>
<font color=#447700>!<a name='601'></font>
IF(SM(I,J).LT.0.5) THEN<a name='602'>
QS(I,J)=PQ0/PSFCK*EXP(A2*(TSFCK-A3)/(TSFCK-A4))<a name='603'>
ELSEIF(SM(I,J).GT.0.5) THEN<a name='604'>
<font color=#447700>!reinstated below 1020<a name='605'></font>
THS(I,J)=SST(I,J)*(1.E5/(PD(I,J)+PDTOP+PT))**CAPA<a name='606'>
ENDIF<a name='607'>
<a name='608'>
IF (THS(I,J) .lt. 200) then<a name='609'>
write(6,*) 'bad THS in start_domain_nmm: ', I,J,THS(I,J)<a name='610'>
endif<a name='611'>
<font color=#447700>!<a name='612'></font>
TERM1=-0.068283/T(I,KOFF+1,J)<a name='613'>
PSHLTR(I,J)=(PD(I,J)+PDTOP+PT)*EXP(TERM1)<a name='614'>
<font color=#447700>!<a name='615'></font>
USTAR(I,J)=0.1<a name='616'>
THZ0(I,J)=THS(I,J)<a name='617'>
QZ0(I,J)=QS(I,J)<a name='618'>
UZ0(I,J)=0.<a name='619'>
VZ0(I,J)=0.<a name='620'>
<font color=#447700>! <a name='621'></font>
ENDDO<a name='622'>
ENDDO<a name='623'>
<font color=#447700>!***<a name='624'></font>
<font color=#447700>!*** INITIALIZE CLOUD FIELDS<a name='625'></font>
<font color=#447700>!***<a name='626'></font>
DO J=JFS,JFE<a name='627'>
DO K=KPS,KPE<a name='628'>
DO I=IFS,IFE<a name='629'>
CWM(I,K,J)=0.<a name='630'>
ENDDO<a name='631'>
ENDDO<a name='632'>
ENDDO<a name='633'>
<font color=#447700>!***<a name='634'></font>
<font color=#447700>!*** INITIALIZE ACCUMULATOR ARRAYS TO ZERO.<a name='635'></font>
<font color=#447700>!***<a name='636'></font>
ARDSW=0.0<a name='637'>
ARDLW=0.0<a name='638'>
ASRFC=0.0<a name='639'>
AVRAIN=0.0<a name='640'>
AVCNVC=0.0<a name='641'>
<font color=#447700>!<a name='642'></font>
DO J=JFS,JFE<a name='643'>
DO I=IFS,IFE<a name='644'>
ACFRCV(I,J)=0.<a name='645'>
NCFRCV(I,J)=0<a name='646'>
ACFRST(I,J)=0.<a name='647'>
NCFRST(I,J)=0<a name='648'>
ACSNOW(I,J)=0.<a name='649'>
ACSNOM(I,J)=0.<a name='650'>
SSROFF(I,J)=0.<a name='651'>
BGROFF(I,J)=0.<a name='652'>
ALWIN(I,J) =0.<a name='653'>
ALWOUT(I,J)=0.<a name='654'>
ALWTOA(I,J)=0.<a name='655'>
ASWIN(I,J) =0.<a name='656'>
ASWOUT(I,J)=0.<a name='657'>
ASWTOA(I,J)=0.<a name='658'>
SFCSHX(I,J)=0.<a name='659'>
SFCLHX(I,J)=0.<a name='660'>
SUBSHX(I,J)=0.<a name='661'>
SNOPCX(I,J)=0.<a name='662'>
SFCUVX(I,J)=0.<a name='663'>
SFCEVP(I,J)=0.<a name='664'>
POTEVP(I,J)=0.<a name='665'>
POTFLX(I,J)=0.<a name='666'>
ENDDO<a name='667'>
ENDDO<a name='668'>
<font color=#447700>!***<a name='669'></font>
<font color=#447700>!*** INITIALIZE SATURATION SPECIFIC HUMIDITY OVER THE WATER.<a name='670'></font>
<font color=#447700>!***<a name='671'></font>
EPS=R_D/R_V<a name='672'>
<font color=#447700>!<a name='673'></font>
DO J=JFS,JFE<a name='674'>
DO I=IFS,IFE<a name='675'>
IF(SM(I,J).GT.0.5)THEN<a name='676'>
CLOGES =-CM1/SST(I,J)-CM2*ALOG10(SST(I,J))+CM3<a name='677'>
ESE = 10.**(CLOGES+2.)<a name='678'>
QS(I,J)= SM(I,J)*EPS*ESE/(PD(I,J)+PDTOP+PT-ESE*(1.-EPS))<a name='679'>
ENDIF<a name='680'>
ENDDO<a name='681'>
ENDDO<a name='682'>
<font color=#447700>!*** <a name='683'></font>
<font color=#447700>!*** INITIALIZE TURBULENT KINETIC ENERGY (TKE) TO A SMALL<a name='684'></font>
<font color=#447700>!*** VALUE (EPSQ2) ABOVE GROUND. SET TKE TO ZERO IN THE<a name='685'></font>
<font color=#447700>!*** THE LOWEST MODEL LAYER. IN THE LOWEST TWO ATMOSPHERIC<a name='686'></font>
<font color=#447700>!*** ETA LAYERS SET TKE TO A SMALL VALUE (Q2INI).<a name='687'></font>
<font color=#447700>!***<a name='688'></font>
DO J=JFS,JFE<a name='689'>
DO K=KPS,KPE-1<a name='690'>
DO I=IFS,IFE<a name='691'>
<font color=#447700>!!!tlb Q2(I,K,J)=HTM(I,K+1,J)*HBM2(I,J)*EPSQ2<a name='692'></font>
Q2(I,K+1,J)=HTM(I,K,J)*HBM2(I,J)*EPSQ2<a name='693'>
ENDDO<a name='694'>
ENDDO<a name='695'>
ENDDO<a name='696'>
<font color=#447700>!<a name='697'></font>
DO J=JFS,JFE<a name='698'>
DO I=IFS,IFE<a name='699'>
Q2(I,KPS,J) = 0.<a name='700'>
LLMH = LMH(I,J)<a name='701'>
Q2(I,KOFF+2,J)= HBM2(I,J)*Q2INI<a name='702'>
Q2(I,KOFF+3,J)= HBM2(I,J)*Q2INI<a name='703'>
ENDDO<a name='704'>
ENDDO<a name='705'>
<font color=#447700>!*** <a name='706'></font>
<font color=#447700>!*** PAD ABOVE GROUND SPECIFIC HUMIDITY IF IT IS TOO SMALL.<a name='707'></font>
<font color=#447700>!*** INITIALIZE LATENT HEATING ACCUMULATION ARRAYS.<a name='708'></font>
<font color=#447700>!***<a name='709'></font>
DO J=JFS,JFE<a name='710'>
DO K=KPS,KPE<a name='711'>
DO I=IFS,IFE<a name='712'>
IF(Q(I,K,J).LT.EPSQ)Q(I,K,J)=EPSQ*HTM(I,K,J)<a name='713'>
TRAIN(I,K,J)=0.<a name='714'>
TCUCN(I,K,J)=0.<a name='715'>
ENDDO<a name='716'>
ENDDO<a name='717'>
ENDDO<a name='718'>
<font color=#447700>!<a name='719'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='720'></font>
<font color=#447700>!*** END OF SCRATCH START INITIALIZATION BLOCK.<a name='721'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='722'></font>
<font color=#447700>!<a name='723'></font>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_15">('INIT: INITIALIZED ARRAYS FOR CLEAN START')<a name='724'>
ENDIF<a name='725'>
<font color=#447700>!<a name='726'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='727'></font>
<font color=#447700>!*** RESTART INITIALIZING. CHECK TO SEE IF WE NEED TO ZERO<a name='728'></font>
<font color=#447700>!*** ACCUMULATION ARRAYS.<a name='729'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='730'></font>
<a name='731'>
TSPH=3600./GRID%DT <font color=#447700>! needed?<a name='732'></font>
NPHS0=GRID%NPHS<a name='733'>
<a name='734'>
<font color=#447700>!!!! How do we pass in this TSTART information in reality???<a name='735'></font>
<a name='736'>
TSTART=00.0<a name='737'>
TPREC=06.0<a name='738'>
THEAT=06.0<a name='739'>
TCLOD=06.0<a name='740'>
TRDSW=03.0<a name='741'>
TRDLW=03.0<a name='742'>
TSRFC=03.0<a name='743'>
<a name='744'>
NSTART = INT(TSTART*TSPH+0.5)<a name='745'>
<a name='746'>
NTSD = NSTART <font color=#447700>! (this NTSD value not honored by integrate)<a name='747'></font>
<a name='748'>
<a name='749'>
<font color=#447700>!! want non-zero values for NPREC, NHEAT type vars to avoid problems<a name='750'></font>
<font color=#447700>!! with mod statements below.<a name='751'></font>
<a name='752'>
NPREC = INT(TPREC *TSPH+0.5)<a name='753'>
NHEAT = INT(THEAT *TSPH+0.5)<a name='754'>
NCLOD = INT(TCLOD *TSPH+0.5)<a name='755'>
NRDSW = INT(TRDSW *TSPH+0.5)<a name='756'>
NRDLW = INT(TRDLW *TSPH+0.5)<a name='757'>
NSRFC = INT(TSRFC *TSPH+0.5)<a name='758'>
<a name='759'>
IF(RESTRT)THEN<a name='760'>
<font color=#447700>!<a name='761'></font>
<font color=#447700>!***<a name='762'></font>
<font color=#447700>!*** AVERAGE CLOUD AMOUNT ARRAY<a name='763'></font>
<font color=#447700>!***<a name='764'></font>
IF(MOD(NTSD,NCLOD).LT.GRID%NPHS)THEN<a name='765'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_16">(' ZERO AVG CLD AMT ARRAY')<a name='766'>
DO J=JFS,JFE<a name='767'>
DO I=IFS,IFE<a name='768'>
ACFRCV(I,J)=0.<a name='769'>
NCFRCV(I,J)=0<a name='770'>
ACFRST(I,J)=0.<a name='771'>
NCFRST(I,J)=0<a name='772'>
ENDDO<a name='773'>
ENDDO<a name='774'>
ENDIF<a name='775'>
<font color=#447700>!*** <a name='776'></font>
<font color=#447700>!*** GRID-SCALE AND CONVECTIVE LATENT HEATING ARRAYS.<a name='777'></font>
<font color=#447700>!*** <a name='778'></font>
IF(MOD(NTSD,NHEAT).LT.GRID%NCNVC)THEN<a name='779'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_17">(' ZERO ACCUM LATENT HEATING ARRAYS')<a name='780'>
<font color=#447700>!<a name='781'></font>
AVRAIN=0.<a name='782'>
AVCNVC=0.<a name='783'>
DO J=JFS,JFE<a name='784'>
DO K=KPS,KPE<a name='785'>
DO I=IFS,IFE<a name='786'>
TRAIN(I,K,J)=0.<a name='787'>
TCUCN(I,K,J)=0.<a name='788'>
ENDDO<a name='789'>
ENDDO<a name='790'>
ENDDO<a name='791'>
ENDIF<a name='792'>
<font color=#447700>!***<a name='793'></font>
<font color=#447700>!*** IF THIS IS NOT A NESTED RUN, INITIALIZE TKE<a name='794'></font>
<font color=#447700>!***<a name='795'></font>
<font color=#447700>! IF(.NOT.NEST)THEN<a name='796'></font>
<font color=#447700>! DO K=1,LM<a name='797'></font>
<font color=#447700>! DO J=JFS,JFE<a name='798'></font>
<font color=#447700>! DO I=IFS,IFE<a name='799'></font>
<font color=#447700>! Q2(I,K,J)=AMAX1(Q2(I,K,J)*HBM2(I,J),EPSQ2)<a name='800'></font>
<font color=#447700>! ENDDO<a name='801'></font>
<font color=#447700>! ENDDO<a name='802'></font>
<font color=#447700>! ENDDO<a name='803'></font>
<font color=#447700>! ENDIF<a name='804'></font>
<font color=#447700>!***<a name='805'></font>
<font color=#447700>!*** CLOUD EFFICIENCY<a name='806'></font>
<font color=#447700>!***<a name='807'></font>
<font color=#447700>! DO J=JFS,JFE<a name='808'></font>
<font color=#447700>! DO I=IFS,IFE<a name='809'></font>
<font color=#447700>!!! CLDEFI(I,J)=AVGEFI*SM(I,J)+STEFI*(1.-SM(I,J))<a name='810'></font>
<font color=#447700>! CLDEFI(I,J)=1.<a name='811'></font>
<font color=#447700>! ENDDO<a name='812'></font>
<font color=#447700>! ENDDO<a name='813'></font>
<font color=#447700>!***<a name='814'></font>
<font color=#447700>!*** TOTAL AND CONVECTIVE PRECIPITATION ARRAYS.<a name='815'></font>
<font color=#447700>!*** TOTAL SNOW AND SNOW MELT ARRAYS.<a name='816'></font>
<font color=#447700>!*** STORM SURFACE AND BASE GROUND RUN OFF ARRAYS.<a name='817'></font>
<font color=#447700>! <a name='818'></font>
IF(MOD(NTSD,NPREC).LT.GRID%NPHS)THEN<a name='819'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_18">(' ZERO ACCUM PRECIP ARRAYS')<a name='820'>
DO J=JFS,JFE<a name='821'>
DO I=IFS,IFE<a name='822'>
ACPREC(I,J)=0.<a name='823'>
CUPREC(I,J)=0.<a name='824'>
ACSNOW(I,J)=0.<a name='825'>
ACSNOM(I,J)=0.<a name='826'>
SSROFF(I,J)=0.<a name='827'>
BGROFF(I,J)=0.<a name='828'>
ENDDO<a name='829'>
ENDDO<a name='830'>
ENDIF<a name='831'>
<font color=#447700>!*** <a name='832'></font>
<font color=#447700>!*** LONG WAVE RADIATION ARRAYS.<a name='833'></font>
<font color=#447700>!*** <a name='834'></font>
IF(MOD(NTSD,NRDLW).LT.GRID%NPHS)THEN<a name='835'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_19">(' ZERO ACCUM LW RADTN ARRAYS')<a name='836'>
ARDLW=0.<a name='837'>
DO J=JFS,JFE<a name='838'>
DO I=IFS,IFE<a name='839'>
ALWIN(I,J) =0.<a name='840'>
ALWOUT(I,J)=0.<a name='841'>
ALWTOA(I,J)=0.<a name='842'>
ENDDO<a name='843'>
ENDDO<a name='844'>
ENDIF<a name='845'>
<font color=#447700>!*** <a name='846'></font>
<font color=#447700>!*** SHORT WAVE RADIATION ARRAYS.<a name='847'></font>
<font color=#447700>!*** <a name='848'></font>
IF(MOD(NTSD,NRDSW).LT.GRID%NPHS)THEN<a name='849'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_20">(' ZERO ACCUM SW RADTN ARRAYS')<a name='850'>
ARDSW=0.<a name='851'>
DO J=JFS,JFE<a name='852'>
DO I=IFS,IFE<a name='853'>
ASWIN(I,J) =0.<a name='854'>
ASWOUT(I,J)=0.<a name='855'>
ASWTOA(I,J)=0.<a name='856'>
ENDDO<a name='857'>
ENDDO<a name='858'>
ENDIF<a name='859'>
<font color=#447700>!*** <a name='860'></font>
<font color=#447700>!*** SURFACE SENSIBLE AND LATENT HEAT FLUX ARRAYS.<a name='861'></font>
<font color=#447700>!*** <a name='862'></font>
IF(MOD(NTSD,NSRFC).LT.GRID%NPHS)THEN<a name='863'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_21">(' ZERO ACCUM SFC FLUX ARRAYS')<a name='864'>
ASRFC=0.<a name='865'>
DO J=JFS,JFE<a name='866'>
DO I=IFS,IFE<a name='867'>
SFCSHX(I,J)=0.<a name='868'>
SFCLHX(I,J)=0.<a name='869'>
SUBSHX(I,J)=0.<a name='870'>
SNOPCX(I,J)=0.<a name='871'>
SFCUVX(I,J)=0.<a name='872'>
SFCEVP(I,J)=0.<a name='873'>
POTEVP(I,J)=0.<a name='874'>
POTFLX(I,J)=0.<a name='875'>
ENDDO<a name='876'>
ENDDO<a name='877'>
ENDIF<a name='878'>
<font color=#447700>!***<a name='879'></font>
<font color=#447700>!*** ENDIF FOR RESTART FILE ACCUMULATION ZERO BLOCK.<a name='880'></font>
<font color=#447700>!***<a name='881'></font>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_22">('INIT: INITIALIZED ARRAYS FOR RESTART START')<a name='882'>
ENDIF<a name='883'>
<font color=#447700>!<a name='884'></font>
DO J=JFS,JFE<a name='885'>
DO K=KPS,KPE<a name='886'>
DO I=IFS,IFE<a name='887'>
ZERO_3D(I,K,J)=0.<a name='888'>
ENDDO<a name='889'>
ENDDO<a name='890'>
ENDDO<a name='891'>
<font color=#447700>!----------------------------------------------------------------------<a name='892'></font>
<font color=#447700>!<a name='893'></font>
<font color=#447700>!*** FLAG FOR INITIALIZING ARRAYS, LOOKUP TABLES, & CONSTANTS USED IN<a name='894'></font>
<font color=#447700>!*** MICROPHYSICS AND RADIATION<a name='895'></font>
<font color=#447700>!<a name='896'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='897'></font>
<font color=#447700>!<a name='898'></font>
MICRO_START=.TRUE.<a name='899'>
<font color=#447700>!<a name='900'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='901'></font>
<font color=#447700>!***<a name='902'></font>
<font color=#447700>!*** INITIALIZE ADVECTION TENDENCIES TO ZERO SO THAT<a name='903'></font>
<font color=#447700>!*** BOUNDARY POINTS WILL ALWAYS BE ZERO<a name='904'></font>
<font color=#447700>!***<a name='905'></font>
DO J=JFS,JFE<a name='906'>
DO K=KPS,KPE<a name='907'>
DO I=IFS,IFE<a name='908'>
ADT(I,K,J)=0.<a name='909'>
ADU(I,K,J)=0.<a name='910'>
ADV(I,K,J)=0.<a name='911'>
ENDDO<a name='912'>
ENDDO<a name='913'>
ENDDO<a name='914'>
<font color=#447700>!----------------------------------------------------------------------<a name='915'></font>
<font color=#447700>!***<a name='916'></font>
<font color=#447700>!*** SET INDEX ARRAYS FOR UPSTREAM ADVECTION<a name='917'></font>
<font color=#447700>!***<a name='918'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='919'></font>
DO J=JFS,JFE<a name='920'>
N_IUP_H(J)=0<a name='921'>
N_IUP_V(J)=0<a name='922'>
N_IUP_ADH(J)=0<a name='923'>
N_IUP_ADV(J)=0<a name='924'>
<font color=#447700>!<a name='925'></font>
DO I=IFS,IFE<a name='926'>
IUP_H(I,J)=-999<a name='927'>
IUP_V(I,J)=-999<a name='928'>
IUP_ADH(I,J)=-999<a name='929'>
IUP_ADV(I,J)=-999<a name='930'>
ENDDO<a name='931'>
<font color=#447700>!<a name='932'></font>
ENDDO<a name='933'>
<a name='934'>
#ifndef NO_UPSTREAM_ADVECTION<a name='935'>
<font color=#447700>!<a name='936'></font>
<font color=#447700>!*** N_IUP_H HOLDS THE NUMBER OF MASS POINTS NEEDED IN EACH ROW<a name='937'></font>
<font color=#447700>!*** FOR UPSTREAM ADVECTION (FULL ROWS IN THE 3RD THROUGH 7TH<a name='938'></font>
<font color=#447700>!*** ROWS FROM THE SOUTH AND NORTH GLOBAL BOUNDARIES AND <a name='939'></font>
<font color=#447700>!*** FOUR POINTS ADJACENT TO THE WEST AND EAST GLOBAL BOUNDARIES<a name='940'></font>
<font color=#447700>!*** ON ALL OTHER INTERNAL ROWS). SIMILARLY FOR N_IUP_V.<a name='941'></font>
<font color=#447700>!*** BECAUSE OF HORIZONTAL OPERATIONS, THESE POINTS EXTEND OUTSIDE <a name='942'></font>
<font color=#447700>!*** OF THE UPSTREAM REGION SOMEWHAT.<a name='943'></font>
<font color=#447700>!*** N_IUP_ADH HOLDS THE NUMBER OF MASS POINTS NEEDED IN EACH ROW<a name='944'></font>
<font color=#447700>!*** FOR THE COMPUTATION OF THE TENDENCIES THEMSELVES (ADT, ADQ2M<a name='945'></font>
<font color=#447700>!*** AND ADQ2L); SPECIFICALLY THESE TENDENCIES ARE ONLY DONE IN<a name='946'></font>
<font color=#447700>!*** THE UPSTREAM REGION.<a name='947'></font>
<font color=#447700>!*** N_IUP_ADV HOLDS THE NUMBER OF MASS POINTS NEEDED IN EACH ROW<a name='948'></font>
<font color=#447700>!*** FOR THE VELOCITY POINT TENDENCIES.<a name='949'></font>
<font color=#447700>!*** IUP_H AND IUP_V HOLD THE ACTUAL I VALUES USED IN EACH ROW.<a name='950'></font>
<font color=#447700>!*** LIKEWISE FOR IUP_ADH AND IUP_ADV. <a name='951'></font>
<font color=#447700>!*** ALSO, SET UPSTRM FOR THOSE TASKS AROUND THE GLOBAL EDGE.<a name='952'></font>
<font color=#447700>!<a name='953'></font>
UPSTRM=.FALSE.<a name='954'>
<font color=#447700>!<a name='955'></font>
S_BDY=(JPS==JDS)<a name='956'>
N_BDY=(JPE==JDE)<a name='957'>
W_BDY=(IPS==IDS)<a name='958'>
E_BDY=(IPE==IDE)<a name='959'>
<font color=#447700>!<a name='960'></font>
JTPAD2=2<a name='961'>
JBPAD2=2<a name='962'>
IRPAD2=2<a name='963'>
ILPAD2=2<a name='964'>
<font color=#447700>!<a name='965'></font>
IF(S_BDY)THEN<a name='966'>
UPSTRM=.TRUE.<a name='967'>
JBPAD2=0<a name='968'>
<font color=#447700>!<a name='969'></font>
DO JJ=1,7<a name='970'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='971'></font>
KNTI=0<a name='972'>
DO I=MYIS_P2,MYIE_P2<a name='973'>
IUP_H(IMS+KNTI,J)=I<a name='974'>
IUP_V(IMS+KNTI,J)=I<a name='975'>
KNTI=KNTI+1<a name='976'>
ENDDO<a name='977'>
N_IUP_H(J)=KNTI<a name='978'>
N_IUP_V(J)=KNTI<a name='979'>
ENDDO<a name='980'>
<font color=#447700>!<a name='981'></font>
DO JJ=3,5<a name='982'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='983'></font>
KNTI=0<a name='984'>
ISTART=MYIS1_P2<a name='985'>
IEND=MYIE1_P2<a name='986'>
IF(E_BDY)IEND=IEND-MOD(JJ+1,2)<a name='987'>
DO I=ISTART,IEND<a name='988'>
IUP_ADH(IMS+KNTI,J)=I<a name='989'>
KNTI=KNTI+1<a name='990'>
ENDDO<a name='991'>
N_IUP_ADH(J)=KNTI<a name='992'>
<font color=#447700>!<a name='993'></font>
KNTI=0<a name='994'>
ISTART=MYIS1_P2<a name='995'>
IEND=MYIE1_P2<a name='996'>
IF(E_BDY)IEND=IEND-MOD(JJ,2)<a name='997'>
DO I=ISTART,IEND<a name='998'>
IUP_ADV(IMS+KNTI,J)=I<a name='999'>
KNTI=KNTI+1<a name='1000'>
ENDDO<a name='1001'>
N_IUP_ADV(J)=KNTI<a name='1002'>
ENDDO<a name='1003'>
ENDIF<a name='1004'>
<font color=#447700>!<a name='1005'></font>
IF(N_BDY)THEN<a name='1006'>
UPSTRM=.TRUE.<a name='1007'>
JTPAD2=0<a name='1008'>
<font color=#447700>!<a name='1009'></font>
DO JJ=JDE-7, JDE-1 <font color=#447700>! JM-6,JM<a name='1010'></font>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1011'></font>
KNTI=0<a name='1012'>
DO I=MYIS_P2,MYIE_P2<a name='1013'>
IUP_H(IMS+KNTI,J)=I<a name='1014'>
IUP_V(IMS+KNTI,J)=I<a name='1015'>
KNTI=KNTI+1<a name='1016'>
ENDDO<a name='1017'>
N_IUP_H(J)=KNTI<a name='1018'>
N_IUP_V(J)=KNTI<a name='1019'>
ENDDO<a name='1020'>
<font color=#447700>!<a name='1021'></font>
DO JJ=JDE-5, JDE-3 <font color=#447700>! JM-4,JM-2<a name='1022'></font>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1023'></font>
KNTI=0<a name='1024'>
ISTART=MYIS1_P2<a name='1025'>
IEND=MYIE1_P2<a name='1026'>
IF(E_BDY)IEND=IEND-MOD(JJ+1,2)<a name='1027'>
DO I=ISTART,IEND<a name='1028'>
IUP_ADH(IMS+KNTI,J)=I<a name='1029'>
KNTI=KNTI+1<a name='1030'>
ENDDO<a name='1031'>
N_IUP_ADH(J)=KNTI<a name='1032'>
<font color=#447700>!<a name='1033'></font>
KNTI=0<a name='1034'>
ISTART=MYIS1_P2<a name='1035'>
IEND=MYIE1_P2<a name='1036'>
IF(E_BDY)IEND=IEND-MOD(JJ,2)<a name='1037'>
DO I=ISTART,IEND<a name='1038'>
IUP_ADV(IMS+KNTI,J)=I<a name='1039'>
KNTI=KNTI+1<a name='1040'>
ENDDO<a name='1041'>
N_IUP_ADV(J)=KNTI<a name='1042'>
ENDDO<a name='1043'>
ENDIF<a name='1044'>
<font color=#447700>!<a name='1045'></font>
IF(W_BDY)THEN<a name='1046'>
UPSTRM=.TRUE.<a name='1047'>
ILPAD2=0<a name='1048'>
DO JJ=8,JDE-8 <font color=#447700>! JM-7<a name='1049'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1050'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1051'></font>
<font color=#447700>!<a name='1052'></font>
DO I=1,4<a name='1053'>
IUP_H(IMS+I-1,J)=I<a name='1054'>
IUP_V(IMS+I-1,J)=I<a name='1055'>
ENDDO<a name='1056'>
N_IUP_H(J)=4<a name='1057'>
N_IUP_V(J)=4<a name='1058'>
ENDIF<a name='1059'>
ENDDO<a name='1060'>
<font color=#447700>!<a name='1061'></font>
DO JJ=6,JDE-6 <font color=#447700>! JM-5<a name='1062'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1063'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1064'></font>
KNTI=0<a name='1065'>
IEND=2+MOD(JJ,2)<a name='1066'>
DO I=2,IEND<a name='1067'>
IUP_ADH(IMS+KNTI,J)=I<a name='1068'>
KNTI=KNTI+1<a name='1069'>
ENDDO<a name='1070'>
N_IUP_ADH(J)=KNTI<a name='1071'>
<font color=#447700>!<a name='1072'></font>
KNTI=0<a name='1073'>
IEND=2+MOD(JJ+1,2)<a name='1074'>
DO I=2,IEND<a name='1075'>
IUP_ADV(IMS+KNTI,J)=I<a name='1076'>
KNTI=KNTI+1<a name='1077'>
ENDDO<a name='1078'>
N_IUP_ADV(J)=KNTI<a name='1079'>
<font color=#447700>!<a name='1080'></font>
ENDIF<a name='1081'>
ENDDO<a name='1082'>
ENDIF<a name='1083'>
<font color=#447700>!<a name='1084'></font>
CALL <A href='../../html_code/frame/module_dm_stubs.F.html#WRF_GET_NPROCX'>WRF_GET_NPROCX</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_GET_NPROCX_1">(INPES)<a name='1085'>
<font color=#447700>!<a name='1086'></font>
IF(E_BDY)THEN<a name='1087'>
UPSTRM=.TRUE.<a name='1088'>
IRPAD2=0<a name='1089'>
DO JJ=8,JDE-8 <font color=#447700>! JM-7<a name='1090'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1091'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1092'></font>
IEND=IM-MOD(JJ+1,2)<a name='1093'>
ISTART=IEND-3<a name='1094'>
<font color=#447700>!<a name='1095'></font>
<font color=#447700>!*** IN CASE THERE IS ONLY A SINGLE GLOBAL TASK IN THE<a name='1096'></font>
<font color=#447700>!*** I DIRECTION THEN WE MUST ADD THE WESTSIDE UPSTREAM<a name='1097'></font>
<font color=#447700>!*** POINTS TO THE EASTSIDE POINTS IN EACH ROW.<a name='1098'></font>
<font color=#447700>!<a name='1099'></font>
KNTI=0<a name='1100'>
IF(INPES.EQ.1)KNTI=N_IUP_H(J)<a name='1101'>
<font color=#447700>!<a name='1102'></font>
DO II=ISTART,IEND<a name='1103'>
I=II <font color=#447700>! -MY_IS_GLB+1<a name='1104'></font>
IUP_H(IMS+KNTI,J)=I<a name='1105'>
KNTI=KNTI+1<a name='1106'>
ENDDO<a name='1107'>
N_IUP_H(J)=KNTI<a name='1108'>
ENDIF<a name='1109'>
ENDDO<a name='1110'>
<font color=#447700>!<a name='1111'></font>
DO JJ=6,JDE-6 <font color=#447700>! JM-5<a name='1112'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1113'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1114'></font>
IEND=IM-1-MOD(JJ+1,2)<a name='1115'>
ISTART=IEND-MOD(JJ,2)<a name='1116'>
KNTI=0<a name='1117'>
IF(INPES.EQ.1)KNTI=N_IUP_ADH(J)<a name='1118'>
DO II=ISTART,IEND<a name='1119'>
I=II <font color=#447700>! -MY_IS_GLB+1<a name='1120'></font>
IUP_ADH(IMS+KNTI,J)=I<a name='1121'>
KNTI=KNTI+1<a name='1122'>
ENDDO<a name='1123'>
N_IUP_ADH(J)=KNTI<a name='1124'>
ENDIF<a name='1125'>
ENDDO<a name='1126'>
<font color=#447700>!***<a name='1127'></font>
DO JJ=8,JDE-8 <font color=#447700>! JM-7<a name='1128'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1129'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1130'></font>
IEND=IM-MOD(JJ,2)<a name='1131'>
ISTART=IEND-3<a name='1132'>
KNTI=0<a name='1133'>
IF(INPES.EQ.1)KNTI=N_IUP_V(J)<a name='1134'>
<font color=#447700>!<a name='1135'></font>
DO II=ISTART,IEND<a name='1136'>
I=II <font color=#447700>! -MY_IS_GLB+1<a name='1137'></font>
IUP_V(IMS+KNTI,J)=I<a name='1138'>
KNTI=KNTI+1<a name='1139'>
ENDDO<a name='1140'>
N_IUP_V(J)=KNTI<a name='1141'>
ENDIF<a name='1142'>
ENDDO<a name='1143'>
<font color=#447700>!<a name='1144'></font>
DO JJ=6,JDE-6 <font color=#447700>! JM-5<a name='1145'></font>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1146'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1147'></font>
IEND=IM-1-MOD(JJ,2)<a name='1148'>
ISTART=IEND-MOD(JJ+1,2)<a name='1149'>
KNTI=0<a name='1150'>
IF(INPES.EQ.1)KNTI=N_IUP_ADV(J)<a name='1151'>
DO II=ISTART,IEND<a name='1152'>
I=II <font color=#447700>! -MY_IS_GLB+1<a name='1153'></font>
IUP_ADV(IMS+KNTI,J)=I<a name='1154'>
KNTI=KNTI+1<a name='1155'>
ENDDO<a name='1156'>
N_IUP_ADV(J)=KNTI<a name='1157'>
ENDIF<a name='1158'>
ENDDO<a name='1159'>
ENDIF<a name='1160'>
<font color=#447700>!----------------------------------------------------------------------<a name='1161'></font>
<font color=#447700>!!!!!!!!!!!!!!!!!!!!tlb<a name='1162'></font>
<font color=#447700>!!!Read in EM and EMT from the original NMM nhb file<a name='1163'></font>
<font color=#447700>!!! call int_get_fresh_handle( retval )<a name='1164'></font>
<font color=#447700>!!! close(retval)<a name='1165'></font>
<font color=#447700>!!! open(unit=retval,file=seeout,form='UNFORMATTED',iostat=ier)<a name='1166'></font>
<font color=#447700>!!!!!!do j=1,128<a name='1167'></font>
<font color=#447700>!!! read(seeout)<a name='1168'></font>
<font color=#447700>!!!!!! read(55)<a name='1169'></font>
<font color=#447700>!!!!!!enddo<a name='1170'></font>
<font color=#447700>!!! read(seeout)dummyx,em,emt<a name='1171'></font>
<font color=#447700>!!!!!!read(55)dummyx,em,emt<a name='1172'></font>
<font color=#447700>!!! close(retval)<a name='1173'></font>
jam=6+2*(JDE-JDS-1-9)<a name='1174'>
<font color=#447700>! read(55)(em(j),j=1,jam),(emt(j),j=1,jam)<a name='1175'></font>
<font color=#447700>!!!!!!!!!!!!!!!!!!!!tlb<a name='1176'></font>
<font color=#447700>!<a name='1177'></font>
<font color=#447700>!*** EXTRACT EM AND EMT FOR THE LOCAL SUBDOMAINS<a name='1178'></font>
<font color=#447700>!<a name='1179'></font>
DO J=MYJS_P5,MYJE_P5<a name='1180'>
EM_LOC(J)=-9.E9<a name='1181'>
EMT_LOC(J)=-9.E9<a name='1182'>
ENDDO<a name='1183'>
<font color=#447700>!!! IF(IBROW==1)THEN<a name='1184'></font>
IF(S_BDY)THEN<a name='1185'>
DO J=3,5<a name='1186'>
EM_LOC(J)=EM(J-2)<a name='1187'>
EMT_LOC(J)=EMT(J-2)<a name='1188'>
ENDDO<a name='1189'>
ENDIF<a name='1190'>
<font color=#447700>!!! IF(ITROW==1)THEN<a name='1191'></font>
IF(N_BDY)THEN<a name='1192'>
KNT=3<a name='1193'>
DO JJ=JDE-5,JDE-3 <font color=#447700>! JM-4,JM-2<a name='1194'></font>
KNT=KNT+1<a name='1195'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1196'></font>
EM_LOC(J)=EM(KNT)<a name='1197'>
EMT_LOC(J)=EMT(KNT)<a name='1198'>
ENDDO<a name='1199'>
ENDIF<a name='1200'>
<font color=#447700>!!! IF(ILCOL==1)THEN<a name='1201'></font>
IF(W_BDY)THEN<a name='1202'>
KNT=6<a name='1203'>
DO JJ=6,JDE-6 <font color=#447700>! JM-5<a name='1204'></font>
KNT=KNT+1<a name='1205'>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1206'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1207'></font>
EM_LOC(J)=EM(KNT)<a name='1208'>
EMT_LOC(J)=EMT(KNT)<a name='1209'>
ENDIF<a name='1210'>
ENDDO<a name='1211'>
ENDIF<a name='1212'>
<font color=#447700>!!! IF(IRCOL==1)THEN<a name='1213'></font>
IF(E_BDY)THEN<a name='1214'>
KNT=6+JDE-11 <font color=#447700>! JM-10<a name='1215'></font>
DO JJ=6,JDE-6 <font color=#447700>! JM-5<a name='1216'></font>
KNT=KNT+1<a name='1217'>
IF(JJ.GE.MY_JS_GLB-2.AND.JJ.LE.MY_JE_GLB+2)THEN<a name='1218'>
J=JJ <font color=#447700>! -MY_JS_GLB+1<a name='1219'></font>
EM_LOC(J)=EM(KNT)<a name='1220'>
EMT_LOC(J)=EMT(KNT)<a name='1221'>
ENDIF<a name='1222'>
ENDDO<a name='1223'>
ENDIF<a name='1224'>
#else<a name='1225'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_23">( 'start_domain_nmm: upstream advection commented out')<a name='1226'>
#endif<a name='1227'>
<font color=#447700>!<a name='1228'></font>
<font color=#447700>!***<a name='1229'></font>
<font color=#447700>!*** SET ZERO-VALUE FOR SOME OUTPUT DIAGNOSTIC ARRAYS<a name='1230'></font>
<font color=#447700>!***<a name='1231'></font>
IF(NSTART.EQ.0)THEN<a name='1232'>
<font color=#447700>!<a name='1233'></font>
DO J=JFS,JFE<a name='1234'>
DO I=IFS,IFE<a name='1235'>
PCTSNO(I,J)=-999.0<a name='1236'>
IF(SM(I,J).LT.0.5)THEN<a name='1237'>
IF(SICE(I,J).GT.0.5)THEN<a name='1238'>
<font color=#447700>!***<a name='1239'></font>
<font color=#447700>!*** SEA-ICE CASE<a name='1240'></font>
<font color=#447700>!***<a name='1241'></font>
SMSTAV(I,J)=1.0<a name='1242'>
SMSTOT(I,J)=1.0<a name='1243'>
SSROFF(I,J)=0.0<a name='1244'>
BGROFF(I,J)=0.0<a name='1245'>
CMC(I,J)=0.0<a name='1246'>
DO NS=1,GRID%NSOIL<a name='1247'>
SMC(I,NS,J)=1.0<a name='1248'>
<font color=#447700>! SH2O(I,NS,J)=0.05<a name='1249'></font>
SH2O(I,NS,J)=1.0<a name='1250'>
ENDDO<a name='1251'>
ENDIF<a name='1252'>
ELSE<a name='1253'>
<font color=#447700>!***<a name='1254'></font>
<font color=#447700>!*** WATER CASE<a name='1255'></font>
<font color=#447700>!***<a name='1256'></font>
SMSTAV(I,J)=1.0<a name='1257'>
SMSTOT(I,J)=1.0<a name='1258'>
SSROFF(I,J)=0.0<a name='1259'>
BGROFF(I,J)=0.0<a name='1260'>
SOILTB(I,J)=273.16<a name='1261'>
GRNFLX(I,J)=0.<a name='1262'>
SUBSHX(I,J)=0.0<a name='1263'>
ACSNOW(I,J)=0.0<a name='1264'>
ACSNOM(I,J)=0.0<a name='1265'>
SNOPCX(I,J)=0.0<a name='1266'>
CMC(I,J)=0.0<a name='1267'>
SNO(I,J)=0.0<a name='1268'>
DO NS=1,GRID%NSOIL<a name='1269'>
SMC(I,NS,J)=1.0<a name='1270'>
STC(I,NS,J)=273.16<a name='1271'>
<font color=#447700>! SH2O(I,NS,J)=0.05<a name='1272'></font>
SH2O(I,NS,J)=1.0<a name='1273'>
ENDDO<a name='1274'>
ENDIF<a name='1275'>
<font color=#447700>!<a name='1276'></font>
ENDDO<a name='1277'>
ENDDO<a name='1278'>
<font color=#447700>!<a name='1279'></font>
APHTIM=0.0<a name='1280'>
ARATIM=0.0<a name='1281'>
ACUTIM=0.0<a name='1282'>
<font color=#447700>!<a name='1283'></font>
ENDIF<a name='1284'>
<font color=#447700>!<a name='1285'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1286'></font>
<font color=#447700>!*** INITIALIZE RADTN VARIABLES<a name='1287'></font>
<font color=#447700>!*** CALCULATE THE NUMBER OF STEPS AT EACH POINT.<a name='1288'></font>
<font color=#447700>!*** THE ARRAY 'LVL' WILL COORDINATE VERTICAL LOCATIONS BETWEEN<a name='1289'></font>
<font color=#447700>!*** THE LIFTED WORKING ARRAYS AND THE FUNDAMENTAL MODEL ARRAYS.<a name='1290'></font>
<font color=#447700>!*** LVL HOLDS THE HEIGHT (IN MODEL LAYERS) OF THE TOPOGRAPHY AT<a name='1291'></font>
<font color=#447700>!*** EACH GRID POINT.<a name='1292'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1293'></font>
<font color=#447700>! <a name='1294'></font>
DO J=JFS,JFE<a name='1295'>
DO I=IFS,IFE<a name='1296'>
LVL(I,J)=LM-LMH(I,J)<a name='1297'>
ENDDO<a name='1298'>
ENDDO<a name='1299'>
<font color=#447700>!***<a name='1300'></font>
<font color=#447700>!*** DETERMINE MODEL LAYER LIMITS FOR HIGH(3), MIDDLE(2),<a name='1301'></font>
<font color=#447700>!*** AND LOW(1) CLOUDS. ALSO FIND MODEL LAYER THAT IS JUST BELOW<a name='1302'></font>
<font color=#447700>!*** (HEIGHT-WISE) 400 MB. (K400)<a name='1303'></font>
<font color=#447700>!*** <a name='1304'></font>
K400=0<a name='1305'>
PSUM=PT<a name='1306'>
SLPM=101325.<a name='1307'>
PDIF=SLPM-PT<a name='1308'>
DO K=1,LM<a name='1309'>
PSUM=PSUM+DETA(K)*PDIF<a name='1310'>
IF(LPTOP(3).EQ.0)THEN<a name='1311'>
IF(PSUM.GT.PHITP)LPTOP(3)=K<a name='1312'>
ELSEIF(LPTOP(2).EQ.0)THEN<a name='1313'>
IF(PSUM.GT.PMDHI)LPTOP(2)=K<a name='1314'>
ELSEIF(K400.EQ.0)THEN<a name='1315'>
IF(PSUM.GT.P400)K400=K<a name='1316'>
ELSEIF(LPTOP(1).EQ.0)THEN<a name='1317'>
IF(PSUM.GT.PLOMD)LPTOP(1)=K<a name='1318'>
ENDIF<a name='1319'>
ENDDO<a name='1320'>
<font color=#447700>!***<a name='1321'></font>
<font color=#447700>!*** CALL GRADFS ONCE TO CALC. CONSTANTS AND GET O3 DATA<a name='1322'></font>
<font color=#447700>!***<a name='1323'></font>
KCCO2=0<a name='1324'>
<font color=#447700>!***<a name='1325'></font>
<font color=#447700>!*** CALCULATE THE MIDLAYER PRESSURES IN THE STANDARD ATMOSPHERE<a name='1326'></font>
<font color=#447700>!***<a name='1327'></font>
PSS=101325.<a name='1328'>
PDIF=PSS-PT<a name='1329'>
<font color=#447700>!<a name='1330'></font>
ALLOCATE(PHALF(LM+1),STAT=I)<a name='1331'>
<font color=#447700>!<a name='1332'></font>
DO K=KPS,KPE-1<a name='1333'>
PHALF(K+1)=AETA(K)*PDIF+PT<a name='1334'>
ENDDO<a name='1335'>
<a name='1336'>
<font color=#447700>!<a name='1337'></font>
PHALF(1)=0.<a name='1338'>
PHALF(LM+1)=PSS<a name='1339'>
<font color=#447700>!***<a name='1340'></font>
<font color=#447700>!!! CALL GRADFS(PHALF,KCCO2,NUNIT_CO2)<a name='1341'></font>
<font color=#447700>!***<a name='1342'></font>
<font color=#447700>!*** CALL SOLARD TO CALCULATE NON-DIMENSIONAL SUN-EARTH DISTANCE<a name='1343'></font>
<font color=#447700>!***<a name='1344'></font>
<font color=#447700>!!! IF(MYPE.EQ.0)CALL SOLARD(SUN_DIST)<a name='1345'></font>
<font color=#447700>!!! CALL MPI_BCAST(SUN_DIST,1,MPI_REAL,0,MPI_COMM_COMP,IRTN)<a name='1346'></font>
<a name='1347'>
<font color=#447700>!***<a name='1348'></font>
<font color=#447700>!*** CALL ZENITH SIMPLY TO GET THE DAY OF THE YEAR FOR<a name='1349'></font>
<font color=#447700>!*** THE SETUP OF THE OZONE DATA<a name='1350'></font>
<font color=#447700>!***<a name='1351'></font>
TIME=(NTSD-1)*GRID%DT<a name='1352'>
<font color=#447700>!<a name='1353'></font>
<font color=#447700>!!! CALL ZENITH(TIME,DAYI,HOUR)<a name='1354'></font>
<font color=#447700>!<a name='1355'></font>
ADDL=0.<a name='1356'>
IF(MOD(IDAT(3),4).EQ.0)ADDL=1.<a name='1357'>
<font color=#447700>!<a name='1358'></font>
<font color=#447700>!!! CALL O3CLIM<a name='1359'></font>
<font color=#447700>!<a name='1360'></font>
<font color=#447700>!<a name='1361'></font>
DEALLOCATE(PHALF)<a name='1362'>
<font color=#447700>!----------------------------------------------------------------------<a name='1363'></font>
<font color=#447700>!*** SOME INITIAL VALUES RELATED TO TURBULENCE SCHEME<a name='1364'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1365'></font>
<font color=#447700>!<a name='1366'></font>
DO J=JFS,JFE<a name='1367'>
DO I=IFS,IFE<a name='1368'>
<font color=#447700>!***<a name='1369'></font>
<font color=#447700>!*** TRY A SIMPLE LINEAR INTERP TO GET 2/10 M VALUES<a name='1370'></font>
<font color=#447700>!***<a name='1371'></font>
PDSL(I,J)=PD(I,J)*RES(I,J)<a name='1372'>
LMHK=LMH(I,J)<a name='1373'>
LMVK=LMV(I,J)<a name='1374'>
<font color=#447700>!<a name='1375'></font>
KOFF=KPE-1-LMHK<a name='1376'>
KOFV=KPE-1-LMVK<a name='1377'>
<font color=#447700>!<a name='1378'></font>
ULM=U(I,KOFV+1,J)<a name='1379'>
VLM=V(I,KOFV+1,J)<a name='1380'>
TLM=T(I,KOFF+1,J)<a name='1381'>
QLM=Q(I,KOFF+1,J)<a name='1382'>
PLM=AETA1(KOFF+1)*PDTOP+AETA2(KOFF+1)*PDSL(I,J)+PT<a name='1383'>
APELM=(1.0E5/PLM)**CAPA<a name='1384'>
APELMNW=(1.0E5/PSHLTR(I,J))**CAPA<a name='1385'>
THLM=TLM*APELM<a name='1386'>
DPLM=(DETA1(KOFF+1)*PDTOP+DETA2(KOFF+1)*PDSL(I,J))*0.5<a name='1387'>
DZLM=R_D*DPLM*TLM*(1.+P608*QLM)/(G*PLM)<a name='1388'>
FAC1=10./DZLM<a name='1389'>
FAC2=(DZLM-10.)/DZLM<a name='1390'>
IF(DZLM.LE.10.)THEN<a name='1391'>
FAC1=1.<a name='1392'>
FAC2=0.<a name='1393'>
ENDIF<a name='1394'>
<font color=#447700>!<a name='1395'></font>
IF(.NOT.RESTRT)THEN<a name='1396'>
TH10(I,J)=FAC2*THS(I,J)+FAC1*THLM<a name='1397'>
Q10(I,J)=FAC2*QS(I,J)+FAC1*QLM<a name='1398'>
U10(I,J)=ULM<a name='1399'>
V10(I,J)=VLM<a name='1400'>
ENDIF<a name='1401'>
<font color=#447700>!<a name='1402'></font>
FAC1=2./DZLM<a name='1403'>
FAC2=(DZLM-2.)/DZLM<a name='1404'>
IF(DZLM.LE.2.)THEN<a name='1405'>
FAC1=1.<a name='1406'>
FAC2=0.<a name='1407'>
ENDIF<a name='1408'>
<font color=#447700>!<a name='1409'></font>
IF(.NOT.RESTRT.OR.NEST)THEN<a name='1410'>
<font color=#447700>!mp TSHLTR(I,J)=FAC2*THS(I,J)+FAC1*THLM<a name='1411'></font>
TSHLTR(I,J)=0.2*THS(I,J)+0.8*THLM<a name='1412'>
<font color=#447700>!mp QSHLTR(I,J)=FAC2*QS(I,J)+FAC1*QLM<a name='1413'></font>
QSHLTR(I,J)=0.2*QS(I,J)+0.8*QLM<a name='1414'>
ENDIF<a name='1415'>
<font color=#447700>!***<a name='1416'></font>
<font color=#447700>!*** NEED TO CONVERT TO THETA IF IS THE RESTART CASE<a name='1417'></font>
<font color=#447700>!*** AS CHKOUT.f WILL CONVERT TO TEMPERATURE<a name='1418'></font>
<font color=#447700>!***<a name='1419'></font>
IF(RESTRT)THEN<a name='1420'>
TSHLTR(I,J)=TSHLTR(I,J)*APELMNW<a name='1421'>
ENDIF<a name='1422'>
ENDDO<a name='1423'>
ENDDO<a name='1424'>
<font color=#447700>!<a name='1425'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1426'></font>
<font color=#447700>!*** INITIALIZE TAU-1 VALUES FOR ADAMS-BASHFORTH <a name='1427'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1428'></font>
<font color=#447700>!<a name='1429'></font>
DO J=jfs,jfe<a name='1430'>
DO K=KPS,KPE<a name='1431'>
DO I=ifs,ife<a name='1432'>
TOLD(I,K,J)=T(I,K,J) <font color=#447700>! T AT TAU-1<a name='1433'></font>
UOLD(I,K,J)=U(I,K,J) <font color=#447700>! U AT TAU-1<a name='1434'></font>
VOLD(I,K,J)=V(I,K,J) <font color=#447700>! V AT TAU-1<a name='1435'></font>
ENDDO<a name='1436'>
ENDDO<a name='1437'>
ENDDO<a name='1438'>
<font color=#447700>!<a name='1439'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1440'></font>
<font color=#447700>!*** INITIALIZE NONHYDROSTATIC QUANTITIES<a name='1441'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1442'></font>
<font color=#447700>!<a name='1443'></font>
<font color=#447700>!!!! SHOULD DWDT BE REDEFINED IF RESTRT?<a name='1444'></font>
<a name='1445'>
IF(.NOT.RESTRT.OR.NEST)THEN<a name='1446'>
DO J=jfs,jfe<a name='1447'>
DO K=KPS,KPE<a name='1448'>
DO I=ifs,ife<a name='1449'>
DWDT(I,K,J)=1.<a name='1450'>
ENDDO<a name='1451'>
ENDDO<a name='1452'>
ENDDO<a name='1453'>
ENDIF<a name='1454'>
<font color=#447700>!***<a name='1455'></font>
IF(GRID%SIGMA.EQ.1)THEN<a name='1456'>
DO J=jfs,jfe<a name='1457'>
DO I=ifs,ife<a name='1458'>
PDSL(I,J)=PD(I,J)<a name='1459'>
ENDDO<a name='1460'>
ENDDO<a name='1461'>
ELSE<a name='1462'>
DO J=jfs,jfe<a name='1463'>
DO I=ifs,ife<a name='1464'>
PDSL(I,J)=RES(I,J)*PD(I,J)<a name='1465'>
ENDDO<a name='1466'>
ENDDO<a name='1467'>
ENDIF<a name='1468'>
<font color=#447700>!<a name='1469'></font>
<font color=#447700>!***<a name='1470'></font>
<font color=#447700>!<a name='1471'></font>
<font color=#447700>!<a name='1472'></font>
<font color=#447700>!!!! SHOULD PINT,Z,W BE REDEFINED IF RESTRT?<a name='1473'></font>
<a name='1474'>
write(0,*)' restrt=',restrt,' nest=',nest<a name='1475'>
write(0,*)' ifs=',ifs,' ife=',ife<a name='1476'>
write(0,*)' jfs=',jfs,' jfe=',jfe<a name='1477'>
write(0,*)' kps=',kps,' kpe=',kpe<a name='1478'>
write(0,*)' pdtop=',pdtop,' pt=',pt<a name='1479'>
IF(.NOT.RESTRT.OR.NEST)THEN<a name='1480'>
DO J=jfs,jfe<a name='1481'>
DO K=KPS,KPE<a name='1482'>
DO I=ifs,ife<a name='1483'>
PINT(I,K,J)=ETA1(K)*PDTOP+ETA2(K)*PDSL(I,J)+PT<a name='1484'>
Z(I,K,J)=PINT(I,K,J)<a name='1485'>
W(I,K,J)=0.<a name='1486'>
ENDDO<a name='1487'>
ENDDO<a name='1488'>
ENDDO<a name='1489'>
ENDIF<a name='1490'>
<a name='1491'>
#ifndef NO_RESTRICT_ACCEL<a name='1492'>
<font color=#447700>!----------------------------------------------------------------------<a name='1493'></font>
<font color=#447700>!*** RESTRICTING THE ACCELERATION ALONG THE BOUNDARIES<a name='1494'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1495'></font>
<font color=#447700>!<a name='1496'></font>
DO J=jfs,jfe<a name='1497'>
DO I=ifs,ife<a name='1498'>
DWDTMN(I,J)=-EPSIN<a name='1499'>
DWDTMX(I,J)= EPSIN<a name='1500'>
ENDDO<a name='1501'>
ENDDO<a name='1502'>
<a name='1503'>
<a name='1504'>
<font color=#447700>!<a name='1505'></font>
<font color=#447700>!***<a name='1506'></font>
IF(JHL.GT.1)THEN<a name='1507'>
JHH=JDE-1-JHL+1 <font color=#447700>! JM-JHL+1<a name='1508'></font>
IHL=JHL/2+1<a name='1509'>
<font color=#447700>!<a name='1510'></font>
DO J=1,JHL<a name='1511'>
IF(J.GE.MY_JS_GLB-JBPAD2.AND.J.LE.MY_JE_GLB+JTPAD2)THEN<a name='1512'>
JX=J <font color=#447700>! -MY_JS_GLB+1<a name='1513'></font>
DO I=1,IDE-1 <font color=#447700>! IM<a name='1514'></font>
IF(I.GE.MY_IS_GLB-ILPAD2.AND.I.LE.MY_IE_GLB+IRPAD2)THEN<a name='1515'>
IX=I <font color=#447700>! -MY_IS_GLB+1<a name='1516'></font>
DWDTMN(IX,JX)=-EPSB<a name='1517'>
DWDTMX(IX,JX)= EPSB<a name='1518'>
ENDIF<a name='1519'>
ENDDO<a name='1520'>
ENDIF<a name='1521'>
ENDDO<a name='1522'>
<font color=#447700>!<a name='1523'></font>
DO J=JHH,JDE-1 <font color=#447700>! JM<a name='1524'></font>
IF(J.GE.MY_JS_GLB-JBPAD2.AND.J.LE.MY_JE_GLB+JTPAD2)THEN<a name='1525'>
JX=J <font color=#447700>! -MY_JS_GLB+1<a name='1526'></font>
DO I=1,IDE-1 <font color=#447700>! IM<a name='1527'></font>
IF(I.GE.MY_IS_GLB-ILPAD2.AND.I.LE.MY_IE_GLB+IRPAD2)THEN<a name='1528'>
IX=I <font color=#447700>! -MY_IS_GLB+1<a name='1529'></font>
DWDTMN(IX,JX)=-EPSB<a name='1530'>
DWDTMX(IX,JX)= EPSB<a name='1531'>
ENDIF<a name='1532'>
ENDDO<a name='1533'>
ENDIF<a name='1534'>
ENDDO<a name='1535'>
<font color=#447700>!<a name='1536'></font>
DO J=1,JDE-1 <font color=#447700>! JM<a name='1537'></font>
IF(J.GE.MY_JS_GLB-JBPAD2.AND.J.LE.MY_JE_GLB+JTPAD2)THEN<a name='1538'>
JX=J <font color=#447700>! -MY_JS_GLB+1<a name='1539'></font>
DO I=1,IHL<a name='1540'>
IF(I.GE.MY_IS_GLB-ILPAD2.AND.I.LE.MY_IE_GLB+IRPAD2)THEN<a name='1541'>
IX=I <font color=#447700>! -MY_IS_GLB+1<a name='1542'></font>
DWDTMN(IX,JX)=-EPSB<a name='1543'>
DWDTMX(IX,JX)= EPSB<a name='1544'>
ENDIF<a name='1545'>
ENDDO<a name='1546'>
ENDIF<a name='1547'>
ENDDO<a name='1548'>
<font color=#447700>!<a name='1549'></font>
DO J=1,JDE-1 <font color=#447700>! JM<a name='1550'></font>
IF(J.GE.MY_JS_GLB-JBPAD2.AND.J.LE.MY_JE_GLB+JTPAD2)THEN<a name='1551'>
JX=J <font color=#447700>! -MY_JS_GLB+1<a name='1552'></font>
<font color=#447700>! moved this line to inside the J-loop, 20030429, jm<a name='1553'></font>
IHH=IDE-1-IHL+MOD(j,2) <font color=#447700>! IM-IHL+MOD(J,2)<a name='1554'></font>
DO I=IHH,IDE-1 <font color=#447700>! IM<a name='1555'></font>
IF(I.GE.MY_IS_GLB-ILPAD2.AND.I.LE.MY_IE_GLB+IRPAD2)THEN<a name='1556'>
IX=I <font color=#447700>! -MY_IS_GLB+1<a name='1557'></font>
DWDTMN(IX,JX)=-EPSB<a name='1558'>
DWDTMX(IX,JX)= EPSB<a name='1559'>
ENDIF<a name='1560'>
ENDDO<a name='1561'>
ENDIF<a name='1562'>
ENDDO<a name='1563'>
<font color=#447700>!<a name='1564'></font>
ENDIF<a name='1565'>
<a name='1566'>
#else<a name='1567'>
CALL <A href='../../html_code/frame/module_wrf_error.F.html#WRF_MESSAGE'>wrf_message</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="WRF_MESSAGE_24">('start_domain_nmm: NO_RESTRICT_ACCEL')<a name='1568'>
#endif<a name='1569'>
<a name='1570'>
<font color=#447700>!-----------------------------------------------------------------------<a name='1571'></font>
<font color=#447700>!*** CALL THE GENERAL PHYSICS INITIALIZATION<a name='1572'></font>
<font color=#447700>!-----------------------------------------------------------------------<a name='1573'></font>
<font color=#447700>!<a name='1574'></font>
<a name='1575'>
ALLOCATE(SFULL(KMS:KME),STAT=I) ; SFULL = 0.<a name='1576'>
ALLOCATE(SMID(KMS:KME),STAT=I) ; SMID = 0.<a name='1577'>
ALLOCATE(EMISS(IMS:IME,JMS:JME),STAT=I) ; EMISS = 0.<a name='1578'>
ALLOCATE(GLW(IMS:IME,JMS:JME),STAT=I) ; GLW = 0.<a name='1579'>
ALLOCATE(GSW(IMS:IME,JMS:JME),STAT=I) ; GSW = 0.<a name='1580'>
ALLOCATE(HFX(IMS:IME,JMS:JME),STAT=I) ; HFX = 0.<a name='1581'>
ALLOCATE(LOWLYR(IMS:IME,JMS:JME),STAT=I) ; LOWLYR = 0.<a name='1582'>
ALLOCATE(MAVAIL(IMS:IME,JMS:JME),STAT=I) ; MAVAIL = 0.<a name='1583'>
ALLOCATE(NCA(IMS:IME,JMS:JME),STAT=I) ; NCA = 0.<a name='1584'>
ALLOCATE(QFX(IMS:IME,JMS:JME),STAT=I) ; QFX = 0.<a name='1585'>
ALLOCATE(RAINBL(IMS:IME,JMS:JME),STAT=I) ; RAINBL = 0.<a name='1586'>
ALLOCATE(RAINC(IMS:IME,JMS:JME),STAT=I) ; RAINC = 0.<a name='1587'>
ALLOCATE(RAINCV(IMS:IME,JMS:JME),STAT=I) ; RAINCV = 0.<a name='1588'>
ALLOCATE(RAINNC(IMS:IME,JMS:JME),STAT=I) ; RAINNC = 0.<a name='1589'>
ALLOCATE(RAINNCV(IMS:IME,JMS:JME),STAT=I) ; RAINNCV = 0.<a name='1590'>
<a name='1591'>
ALLOCATE(ZS(KMS:KME),STAT=I) ; ZS = 0.<a name='1592'>
ALLOCATE(SNOWC(IMS:IME,JMS:JME),STAT=I) ; SNOWC = 0.<a name='1593'>
ALLOCATE(THC(IMS:IME,JMS:JME),STAT=I) ; THC = 0.<a name='1594'>
ALLOCATE(TMN(IMS:IME,JMS:JME),STAT=I) ; TMN = 0.<a name='1595'>
ALLOCATE(TSFC(IMS:IME,JMS:JME),STAT=I) ; TSFC = 0.<a name='1596'>
ALLOCATE(XLAND(IMS:IME,JMS:JME),STAT=I) ; XLAND = 0.<a name='1597'>
ALLOCATE(XLAT(IMS:IME,JMS:JME),STAT=I) ; XLAT = 0.<a name='1598'>
ALLOCATE(XLONG(IMS:IME,JMS:JME),STAT=I) ; XLONG = 0.<a name='1599'>
ALLOCATE(Z0_DUM(IMS:IME,JMS:JME),STAT=I) ; Z0_DUM = 0.<a name='1600'>
<a name='1601'>
ALLOCATE(DZS(KMS:KME),STAT=I) ; DZS = 0.<a name='1602'>
ALLOCATE(RQCBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQCBLTEN = 0.<a name='1603'>
ALLOCATE(RQIBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQIBLTEN = 0.<a name='1604'>
ALLOCATE(RQVBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQVBLTEN = 0.<a name='1605'>
ALLOCATE(RTHBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RTHBLTEN = 0.<a name='1606'>
ALLOCATE(RUBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RUBLTEN = 0.<a name='1607'>
ALLOCATE(RVBLTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RVBLTEN = 0.<a name='1608'>
ALLOCATE(RQCCUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQCCUTEN = 0.<a name='1609'>
ALLOCATE(RQICUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQICUTEN = 0.<a name='1610'>
ALLOCATE(RQRCUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQRCUTEN = 0.<a name='1611'>
ALLOCATE(RQSCUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQSCUTEN = 0.<a name='1612'>
ALLOCATE(RQVCUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RQVCUTEN = 0.<a name='1613'>
ALLOCATE(RTHCUTEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RTHCUTEN = 0.<a name='1614'>
ALLOCATE(RTHRATEN(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RTHRATEN = 0.<a name='1615'>
ALLOCATE(RTHRATENLW(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RTHRATENLW = 0.<a name='1616'>
ALLOCATE(RTHRATENSW(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; RTHRATENSW = 0.<a name='1617'>
ALLOCATE(TSLB(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; TSLB = 0.<a name='1618'>
#if 0<a name='1619'>
ALLOCATE(W0AVG(IMS:IME,KMS:KME,JMS:JME),STAT=I) ; W0AVG = 0.<a name='1620'>
#endif<a name='1621'>
<font color=#447700>!-----------------------------------------------------------------------<a name='1622'></font>
<font color=#447700>!jm added set of g_inv<a name='1623'></font>
G_INV=1./G<a name='1624'>
ROG=R_D*G_INV<a name='1625'>
GRID%RADT=GRID%NRADS*GRID%DT/60.<a name='1626'>
GRID%BLDT=GRID%NPHS*GRID%DT/60.<a name='1627'>
GRID%CUDT=GRID%NCNVC*GRID%DT/60.<a name='1628'>
GRID%GSMDT=GRID%NPHS*GRID%DT/60.<a name='1629'>
<font color=#447700>!<a name='1630'></font>
DO J=MYJS,MYJE<a name='1631'>
DO I=MYIS,MYIE<a name='1632'>
SFCZ=FIS(I,J)*G_INV<a name='1633'>
<font color=#447700>!!!! ZINT(I,KTS,J)=SFCZ<a name='1634'></font>
PDSL(I,J)=PD(I,J)*RES(I,J)<a name='1635'>
PSURF=PINT(I,KTS,J)<a name='1636'>
EXNSFC=(1.E5/PSURF)**CAPA<a name='1637'>
XLAND(I,J)=SM(I,J)+1.<a name='1638'>
THSIJ=(SST(I,J)*EXNSFC)*(XLAND(I,J)-1.) &<a name='1639'>
& +THS(I,J)*(2.-SM(I,J))<a name='1640'>
TSFC(I,J)=THSIJ/EXNSFC<a name='1641'>
<font color=#447700>!<a name='1642'></font>
DO K=KTS,KTE-1<a name='1643'>
PLYR=(PINT(I,K,J)+PINT(I,K+1,J))*0.5<a name='1644'>
TL=T(I,K,J)<a name='1645'>
CWML=CWM(I,K,J)<a name='1646'>
<font color=#447700>!!! ZINT(I,K+1,J)=ZINT(I,K,J)+TL/PLYR &<a name='1647'></font>
<font color=#447700>!!! & *(DETA1(K)*PDTOP+DETA2(K)*PDSL(I,J))*ROG &<a name='1648'></font>
<font color=#447700>!!! & *(Q(I,K,J)*P608-CWML+1.)<a name='1649'></font>
ENDDO<a name='1650'>
<font color=#447700>!<a name='1651'></font>
DO K=KTS,KTE<a name='1652'>
<font color=#447700>!!! ZMID(I,K,J)=0.5*(ZINT(I,K,J)+ZINT(I,K+1,J))<a name='1653'></font>
ENDDO<a name='1654'>
ENDDO<a name='1655'>
ENDDO<a name='1656'>
<font color=#447700>!<a name='1657'></font>
<font color=#447700>!-----------------------------------------------------------------------<a name='1658'></font>
<font color=#447700>!*** RECREATE SIGMA VALUES AT LAYER INTERFACES FOR THE FULL VERTICAL<a name='1659'></font>
<font color=#447700>!*** DOMAIN FROM THICKNESS VALUES FOR THE TWO SUBDOMAINS.<a name='1660'></font>
<font color=#447700>!*** NOTE: KTE=NUMBER OF LAYERS PLUS ONE<a name='1661'></font>
<font color=#447700>!-----------------------------------------------------------------------<a name='1662'></font>
<font color=#447700>!<a name='1663'></font>
write(0,*)' start_domain kte=',kte<a name='1664'>
PDTOT=101325.-PT<a name='1665'>
RPDTOT=1./PDTOT<a name='1666'>
PDBOT=PDTOT-PDTOP<a name='1667'>
SFULL(KTS)=1.<a name='1668'>
SFULL(KTE)=0.<a name='1669'>
dsigsum = 0.<a name='1670'>
DO K=KTS+1,KTE<a name='1671'>
DSIG=(DETA1(K-1)*PDTOP+DETA2(K-1)*PDBOT)*RPDTOT<a name='1672'>
dsigsum=dsigsum+dsig<a name='1673'>
SFULL(K)=SFULL(K-1)-DSIG<a name='1674'>
SMID(K-1)=0.5*(SFULL(K-1)+SFULL(K))<a name='1675'>
ENDDO<a name='1676'>
dsig=(deta1(kte-1)*pdtop+deta2(kte-1)*pdbot)*rpdtot<a name='1677'>
dsigsum=dsigsum+dsig<a name='1678'>
SMID(KTE-1)=0.5*(SFULL(KTE-1)+SFULL(KTE))<a name='1679'>
<font color=#447700>!<a name='1680'></font>
<font color=#447700>!-----------------------------------------------------------------------<a name='1681'></font>
<a name='1682'>
LU_INDEX=IVGTYP<a name='1683'>
<a name='1684'>
DO J=MYJS,MYJE<a name='1685'>
DO I=MYIS,MYIE<a name='1686'>
Z0_DUM(I,J)=Z0(I,J) <font color=#447700>! hold topography component<a name='1687'></font>
ENDDO<a name='1688'>
ENDDO<a name='1689'>
<a name='1690'>
<font color=#447700>! Several arguments are RCONFIG entries in Registry.NMM. Registry no longer<a name='1691'></font>
<font color=#447700>! includes these as dummy arguments or declares them. Access them from <a name='1692'></font>
<font color=#447700>! GRID. JM 20050819<a name='1693'></font>
CALL <A href='../../html_code/phys/module_physics_init.F.html#PHY_INIT'>PHY_INIT</A><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><A NAME="PHY_INIT_2">(GRID,GRID%ID,CONFIG_FLAGS,GRID%DT,sfull,smid &<a name='1694'>
& ,PT,TSFC,GRID%RADT,GRID%BLDT,GRID%CUDT,GRID%GSMDT &<a name='1695'>
& ,RTHCUTEN, RQVCUTEN, RQRCUTEN &<a name='1696'>
& ,RQCCUTEN, RQSCUTEN, RQICUTEN &<a name='1697'>
& ,RUBLTEN,RVBLTEN,RTHBLTEN &<a name='1698'>
& ,RQVBLTEN,RQCBLTEN,RQIBLTEN &<a name='1699'>
& ,RTHRATEN,RTHRATENLW,RTHRATENSW &<a name='1700'>
& ,STEPBL,STEPRA,STEPCU &<a name='1701'>
& ,W0AVG, RAINNC, RAINC, RAINCV, RAINNCV &<a name='1702'>
& ,NCA &<a name='1703'>
& ,CLDEFI,LOWLYR &<a name='1704'>
& ,MASS_FLUX &<a name='1705'>
& ,RTHFTEN, RQVFTEN &<a name='1706'>
& ,CLDFRA,GLW,GSW,EMISS,LU_INDEX &<a name='1707'>
& ,XLAT,XLONG,ALBEDO,ALBBCK &<a name='1708'>
& ,GRID%GMT,GRID%JULYR,GRID%JULDAY &<a name='1709'>
& ,TMN,XLAND,ZNT,Z0,USTAR,MOL,PBLH,TKE_MYJ &<a name='1710'>
& ,EXCH_H,THC,SNOWC,MAVAIL,HFX,QFX,RAINBL &<a name='1711'>
& ,TSLB,ZS,DZS,GRID%num_soil_layers,warm_rain &<a name='1712'>
& ,APR_GR,APR_W,APR_MC,APR_ST,APR_AS &<a name='1713'>
& ,APR_CAPMA,APR_CAPME,APR_CAPMI &<a name='1714'>
& ,XICE,VEGFRA,SNOW,CANWAT,SMSTAV &<a name='1715'>
& ,SMSTOT, SFCRUNOFF,UDRUNOFF,GRDFLX,ACSNOW &<a name='1716'>
& ,ACSNOM,IVGTYP,ISLTYP,SFCEVP,SMOIS &<a name='1717'>
& ,SH2O, SNOWH, SMFR3D & <font color=#447700>! temporary<a name='1718'></font>
& ,GRID%DX,GRID%DY,F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY &<a name='1719'>
& ,MP_RESTART_STATE,TBPVS_STATE,TBPVS0_STATE &<a name='1720'>
& ,.TRUE.,.FALSE. &<a name='1721'>
& ,ids, ide, jds, jde, kds, kde &<a name='1722'>
& ,ims, ime, jms, jme, kms, kme &<a name='1723'>
& ,its, ite, jts, jte, kts, kte)<a name='1724'>
<font color=#447700>!-----------------------------------------------------------------------<a name='1725'></font>
<font color=#447700>!<a name='1726'></font>
DO J=JMS,JME<a name='1727'>
DO I=IMS,IME<a name='1728'>
Z0(I,J)=Z0_DUM(I,J)+0.5*Z0(I,J) <font color=#447700>! add 1/2 of veg Z0 component, <a name='1729'></font>
<font color=#447700>! expecting this code to be called<a name='1730'></font>
<font color=#447700>! both by real and by the model.<a name='1731'></font>
APREC(I,J)=RAINNC(I,J)*1.E-3<a name='1732'>
CUPREC(I,J)=RAINCV(I,J)*1.E-3<a name='1733'>
ENDDO<a name='1734'>
ENDDO<a name='1735'>
<font color=#447700>!<a name='1736'></font>
DEALLOCATE(SFULL)<a name='1737'>
DEALLOCATE(SMID)<a name='1738'>
DEALLOCATE(DZS)<a name='1739'>
DEALLOCATE(EMISS)<a name='1740'>
DEALLOCATE(GLW)<a name='1741'>
DEALLOCATE(GSW)<a name='1742'>
DEALLOCATE(HFX)<a name='1743'>
DEALLOCATE(LOWLYR)<a name='1744'>
DEALLOCATE(MAVAIL)<a name='1745'>
DEALLOCATE(NCA)<a name='1746'>
DEALLOCATE(QFX)<a name='1747'>
DEALLOCATE(RAINBL)<a name='1748'>
DEALLOCATE(RAINC)<a name='1749'>
DEALLOCATE(RAINCV)<a name='1750'>
DEALLOCATE(RAINNC)<a name='1751'>
DEALLOCATE(RAINNCV)<a name='1752'>
DEALLOCATE(RQCBLTEN)<a name='1753'>
DEALLOCATE(RQIBLTEN)<a name='1754'>
DEALLOCATE(RQVBLTEN)<a name='1755'>
DEALLOCATE(RTHBLTEN)<a name='1756'>
DEALLOCATE(RUBLTEN)<a name='1757'>
DEALLOCATE(RVBLTEN)<a name='1758'>
DEALLOCATE(RQCCUTEN)<a name='1759'>
DEALLOCATE(RQICUTEN)<a name='1760'>
DEALLOCATE(RQRCUTEN)<a name='1761'>
DEALLOCATE(RQSCUTEN)<a name='1762'>
DEALLOCATE(RQVCUTEN)<a name='1763'>
DEALLOCATE(RTHCUTEN)<a name='1764'>
DEALLOCATE(RTHRATEN)<a name='1765'>
DEALLOCATE(RTHRATENLW)<a name='1766'>
DEALLOCATE(RTHRATENSW)<a name='1767'>
DEALLOCATE(SNOWC)<a name='1768'>
DEALLOCATE(THC)<a name='1769'>
DEALLOCATE(TMN)<a name='1770'>
DEALLOCATE(TSFC)<a name='1771'>
DEALLOCATE(TSLB)<a name='1772'>
DEALLOCATE(XLAND)<a name='1773'>
DEALLOCATE(XLAT)<a name='1774'>
DEALLOCATE(XLONG)<a name='1775'>
DEALLOCATE(ZS)<a name='1776'>
#if 0<a name='1777'>
DEALLOCATE(W0AVG)<a name='1778'>
#endif<a name='1779'>
<font color=#447700>!-----------------------------------------------------------------------<a name='1780'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1781'></font>
DO J=jfs,jfe<a name='1782'>
DO I=ifs,ife<a name='1783'>
DWDTMN(I,J)=DWDTMN(I,J)*HBM3(I,J)<a name='1784'>
DWDTMX(I,J)=DWDTMX(I,J)*HBM3(I,J)<a name='1785'>
ENDDO<a name='1786'>
ENDDO<a name='1787'>
<font color=#447700>!----------------------------------------------------------------------<a name='1788'></font>
<font color=#447700>!*** INITIALIZE 3RD INDEX IN WORKING ARRAYS USED IN PFDHT, DDAMP, AND <a name='1789'></font>
<font color=#447700>!*** HZADV. THESE ARRAYS MUST HAVE AN EXTENT OF MORE THAN 1 IN J DUE <a name='1790'></font>
<font color=#447700>!*** TO THE MANY DIFFERENCES AND AVERAGES THAT ARE COMPUTED IN J<a name='1791'></font>
<font color=#447700>!*** OR BECAUSE THE ARRAY IS SIMPLY REFERENCED AT MORE THAN ONE J.<a name='1792'></font>
<font color=#447700>!*** THE WORKING "SPACE" SPANS FROM 3 ROWS SOUTH TO 3 ROWS NORTH<a name='1793'></font>
<font color=#447700>!*** OF THE ROW FOR WHICH THE PRIMARY COMPUTATION IS BEING DONE<a name='1794'></font>
<font color=#447700>!*** THUS THE 3RD DIMENSION CAN VARY FROM -3 TO +3 ALTHOUGH ALL OF<a name='1795'></font>
<font color=#447700>!*** THESE ARRAYS DO NOT NEED TO SPAN THAT MANY ROWS. FOR INSTANCE,<a name='1796'></font>
<font color=#447700>!*** SOME OF THE ARRAYS ARE ONLY USED FROM 2 ROWS SOUTH TO 1 ROW <a name='1797'></font>
<font color=#447700>!*** NORTH, OR FROM 1 ROW SOUTH TO THE CENTRAL ROW. AS THE INTEGRATION<a name='1798'></font>
<font color=#447700>!*** MOVES NORTHWARD, THE SOUTHERNMOST I,K SLAB IS DROPPED FOR EACH<a name='1799'></font>
<font color=#447700>!*** WORKING ARRAY AND THE NORTHERNMOST IS GENERATED. SO AS NOT TO<a name='1800'></font>
<font color=#447700>!*** HAVE TO ACTUALLY MOVE ANY OF THE I,K SLABS NORTHWARD, THE 3RD<a name='1801'></font>
<font color=#447700>!*** INDEX IS CYCLED THROUGH THE EXTENT OF EACH ARRAY'S J DIMENSION.<a name='1802'></font>
<font color=#447700>!*** THE FOLLOWING WILL FILL AN ARRAY WITH THE VALUES OF THE 3RD<a name='1803'></font>
<font color=#447700>!*** INDEX FOR EACH THESE VARIATIONS OF J EXTENTS FOR ALL J's IN<a name='1804'></font>
<font color=#447700>!*** THE LOCAL DOMAIN.<a name='1805'></font>
<font color=#447700>!----------------------------------------------------------------------<a name='1806'></font>
<font color=#447700>!<a name='1807'></font>
<font color=#447700>!*** CASE 0: J EXTENT IS -3 TO 3<a name='1808'></font>
<font color=#447700>!<a name='1809'></font>
KNT=0<a name='1810'>
DO J=MYJS2_P2,MYJE2_P2<a name='1811'>
KNT=KNT+1<a name='1812'>
JP3=KNT+2-7*((KNT+5)/7)<a name='1813'>
JP2=JP3-1+7*((4-JP3)/7)<a name='1814'>
JP1=JP2-1+7*((4-JP2)/7)<a name='1815'>
J00=JP1-1+7*((4-JP1)/7)<a name='1816'>
JM1=J00-1+7*((4-J00)/7)<a name='1817'>
JM2=JM1-1+7*((4-JM1)/7)<a name='1818'>
JM3=JM2-1+7*((4-JM2)/7)<a name='1819'>
INDX3_WRK(3,KNT,0)=JP3<a name='1820'>
INDX3_WRK(2,KNT,0)=JP2<a name='1821'>
INDX3_WRK(1,KNT,0)=JP1<a name='1822'>
INDX3_WRK(0,KNT,0)=J00<a name='1823'>
INDX3_WRK(-1,KNT,0)=JM1<a name='1824'>
INDX3_WRK(-2,KNT,0)=JM2<a name='1825'>
INDX3_WRK(-3,KNT,0)=JM3<a name='1826'>
ENDDO<a name='1827'>
<font color=#447700>!<a name='1828'></font>
<font color=#447700>!*** CASE 1: J EXTENT IS -2 TO 2<a name='1829'></font>
<font color=#447700>!<a name='1830'></font>
KNT=0<a name='1831'>
DO J=MYJS2_P2,MYJE2_P2<a name='1832'>
KNT=KNT+1<a name='1833'>
JP2=KNT+1-5*((KNT+3)/5)<a name='1834'>
JP1=JP2-1+5*((3-JP2)/5)<a name='1835'>
J00=JP1-1+5*((3-JP1)/5)<a name='1836'>
JM1=J00-1+5*((3-J00)/5)<a name='1837'>
JM2=JM1-1+5*((3-JM1)/5)<a name='1838'>
INDX3_WRK(3,KNT,1)=999<a name='1839'>
INDX3_WRK(2,KNT,1)=JP2<a name='1840'>
INDX3_WRK(1,KNT,1)=JP1<a name='1841'>
INDX3_WRK(0,KNT,1)=J00<a name='1842'>
INDX3_WRK(-1,KNT,1)=JM1<a name='1843'>
INDX3_WRK(-2,KNT,1)=JM2<a name='1844'>
INDX3_WRK(-3,KNT,1)=999<a name='1845'>
ENDDO<a name='1846'>
<font color=#447700>!<a name='1847'></font>
<font color=#447700>!*** CASE 2: J EXTENT IS -2 TO 1<a name='1848'></font>
<font color=#447700>!<a name='1849'></font>
KNT=0<a name='1850'>
DO J=MYJS2_P2,MYJE2_P2<a name='1851'>
KNT=KNT+1<a name='1852'>
JP1=KNT-4*((KNT+2)/4)<a name='1853'>
J00=JP1-1+4*((2-JP1)/4)<a name='1854'>
JM1=J00-1+4*((2-J00)/4)<a name='1855'>
JM2=JM1-1+4*((2-JM1)/4)<a name='1856'>
INDX3_WRK(3,KNT,2)=999<a name='1857'>
INDX3_WRK(2,KNT,2)=999<a name='1858'>
INDX3_WRK(1,KNT,2)=JP1<a name='1859'>
INDX3_WRK(0,KNT,2)=J00<a name='1860'>
INDX3_WRK(-1,KNT,2)=JM1<a name='1861'>
INDX3_WRK(-2,KNT,2)=JM2<a name='1862'>
INDX3_WRK(-3,KNT,2)=999<a name='1863'>
ENDDO<a name='1864'>
<font color=#447700>!<a name='1865'></font>
<font color=#447700>!*** CASE 3: J EXTENT IS -1 TO 2<a name='1866'></font>
<font color=#447700>!<a name='1867'></font>
KNT=0<a name='1868'>
DO J=MYJS2_P2,MYJE2_P2<a name='1869'>
KNT=KNT+1<a name='1870'>
JP2=KNT+1-4*((KNT+2)/4)<a name='1871'>
JP1=JP2-1+4*((3-JP2)/4)<a name='1872'>
J00=JP1-1+4*((3-JP1)/4)<a name='1873'>
JM1=J00-1+4*((3-J00)/4)<a name='1874'>
INDX3_WRK(3,KNT,3)=999<a name='1875'>
INDX3_WRK(2,KNT,3)=JP2<a name='1876'>
INDX3_WRK(1,KNT,3)=JP1<a name='1877'>
INDX3_WRK(0,KNT,3)=J00<a name='1878'>
INDX3_WRK(-1,KNT,3)=JM1<a name='1879'>
INDX3_WRK(-2,KNT,3)=999<a name='1880'>
INDX3_WRK(-3,KNT,3)=999<a name='1881'>
ENDDO<a name='1882'>
<font color=#447700>!<a name='1883'></font>
<font color=#447700>!*** CASE 4: J EXTENT IS -1 TO 1<a name='1884'></font>
<font color=#447700>!<a name='1885'></font>
KNT=0<a name='1886'>
DO J=MYJS2_P2,MYJE2_P2<a name='1887'>
KNT=KNT+1<a name='1888'>
JP1=KNT-3*((KNT+1)/3)<a name='1889'>
J00=JP1-1+3*((2-JP1)/3)<a name='1890'>
JM1=J00-1+3*((2-J00)/3)<a name='1891'>
INDX3_WRK(3,KNT,4)=999<a name='1892'>
INDX3_WRK(2,KNT,4)=999<a name='1893'>
INDX3_WRK(1,KNT,4)=JP1<a name='1894'>
INDX3_WRK(0,KNT,4)=J00<a name='1895'>
INDX3_WRK(-1,KNT,4)=JM1<a name='1896'>
INDX3_WRK(-2,KNT,4)=999<a name='1897'>
INDX3_WRK(-3,KNT,4)=999<a name='1898'>
ENDDO<a name='1899'>
<font color=#447700>!<a name='1900'></font>
<font color=#447700>!*** CASE 5: J EXTENT IS -1 TO 0<a name='1901'></font>
<font color=#447700>!<a name='1902'></font>
KNT=0<a name='1903'>
DO J=MYJS2_P2,MYJE2_P2<a name='1904'>
KNT=KNT+1<a name='1905'>
J00=-MOD(KNT+1,2)<a name='1906'>
JM1=-1-J00<a name='1907'>
INDX3_WRK(3,KNT,5)=999<a name='1908'>
INDX3_WRK(2,KNT,5)=999<a name='1909'>
INDX3_WRK(1,KNT,5)=999<a name='1910'>
INDX3_WRK(0,KNT,5)=J00<a name='1911'>
INDX3_WRK(-1,KNT,5)=JM1<a name='1912'>
INDX3_WRK(-2,KNT,5)=999<a name='1913'>
INDX3_WRK(-3,KNT,5)=999<a name='1914'>
ENDDO<a name='1915'>
<font color=#447700>!<a name='1916'></font>
<font color=#447700>!*** CASE 6: J EXTENT IS 0 TO 1<a name='1917'></font>
<font color=#447700>!<a name='1918'></font>
KNT=0<a name='1919'>
DO J=MYJS2_P2,MYJE2_P2<a name='1920'>
KNT=KNT+1<a name='1921'>
JP1=MOD(KNT,2)<a name='1922'>
J00=1-JP1<a name='1923'>
INDX3_WRK(3,KNT,6)=999<a name='1924'>
INDX3_WRK(2,KNT,6)=999<a name='1925'>
INDX3_WRK(1,KNT,6)=JP1<a name='1926'>
INDX3_WRK(0,KNT,6)=J00<a name='1927'>
INDX3_WRK(-1,KNT,6)=999<a name='1928'>
INDX3_WRK(-2,KNT,6)=999<a name='1929'>
INDX3_WRK(-3,KNT,6)=999<a name='1930'>
ENDDO<a name='1931'>
<a name='1932'>
#ifdef DM_PARALLEL<a name='1933'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_1.inc.html'>HALO_NMM_INIT_1.inc</A>><A NAME="HALO_NMM_INIT_1.inc_45"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1934'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_2.inc.html'>HALO_NMM_INIT_2.inc</A>><A NAME="HALO_NMM_INIT_2.inc_46"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1935'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_3.inc.html'>HALO_NMM_INIT_3.inc</A>><A NAME="HALO_NMM_INIT_3.inc_47"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1936'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_4.inc.html'>HALO_NMM_INIT_4.inc</A>><A NAME="HALO_NMM_INIT_4.inc_48"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1937'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_5.inc.html'>HALO_NMM_INIT_5.inc</A>><A NAME="HALO_NMM_INIT_5.inc_49"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1938'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_6.inc.html'>HALO_NMM_INIT_6.inc</A>><A NAME="HALO_NMM_INIT_6.inc_50"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1939'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_7.inc.html'>HALO_NMM_INIT_7.inc</A>><A NAME="HALO_NMM_INIT_7.inc_51"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1940'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_8.inc.html'>HALO_NMM_INIT_8.inc</A>><A NAME="HALO_NMM_INIT_8.inc_52"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1941'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_9.inc.html'>HALO_NMM_INIT_9.inc</A>><A NAME="HALO_NMM_INIT_9.inc_53"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1942'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_10.inc.html'>HALO_NMM_INIT_10.inc</A>><A NAME="HALO_NMM_INIT_10.inc_54"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1943'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_11.inc.html'>HALO_NMM_INIT_11.inc</A>><A NAME="HALO_NMM_INIT_11.inc_55"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1944'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_12.inc.html'>HALO_NMM_INIT_12.inc</A>><A NAME="HALO_NMM_INIT_12.inc_56"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1945'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_13.inc.html'>HALO_NMM_INIT_13.inc</A>><A NAME="HALO_NMM_INIT_13.inc_57"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1946'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_14.inc.html'>HALO_NMM_INIT_14.inc</A>><A NAME="HALO_NMM_INIT_14.inc_58"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1947'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_15.inc.html'>HALO_NMM_INIT_15.inc</A>><A NAME="HALO_NMM_INIT_15.inc_59"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1948'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_16.inc.html'>HALO_NMM_INIT_16.inc</A>><A NAME="HALO_NMM_INIT_16.inc_60"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1949'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_17.inc.html'>HALO_NMM_INIT_17.inc</A>><A NAME="HALO_NMM_INIT_17.inc_61"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1950'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_18.inc.html'>HALO_NMM_INIT_18.inc</A>><A NAME="HALO_NMM_INIT_18.inc_62"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1951'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_19.inc.html'>HALO_NMM_INIT_19.inc</A>><A NAME="HALO_NMM_INIT_19.inc_63"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1952'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_20.inc.html'>HALO_NMM_INIT_20.inc</A>><A NAME="HALO_NMM_INIT_20.inc_64"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1953'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_21.inc.html'>HALO_NMM_INIT_21.inc</A>><A NAME="HALO_NMM_INIT_21.inc_65"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1954'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_22.inc.html'>HALO_NMM_INIT_22.inc</A>><A NAME="HALO_NMM_INIT_22.inc_66"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1955'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_23.inc.html'>HALO_NMM_INIT_23.inc</A>><A NAME="HALO_NMM_INIT_23.inc_67"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1956'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_24.inc.html'>HALO_NMM_INIT_24.inc</A>><A NAME="HALO_NMM_INIT_24.inc_68"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1957'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_25.inc.html'>HALO_NMM_INIT_25.inc</A>><A NAME="HALO_NMM_INIT_25.inc_69"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1958'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_26.inc.html'>HALO_NMM_INIT_26.inc</A>><A NAME="HALO_NMM_INIT_26.inc_70"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1959'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_27.inc.html'>HALO_NMM_INIT_27.inc</A>><A NAME="HALO_NMM_INIT_27.inc_71"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1960'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_28.inc.html'>HALO_NMM_INIT_28.inc</A>><A NAME="HALO_NMM_INIT_28.inc_72"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1961'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_29.inc.html'>HALO_NMM_INIT_29.inc</A>><A NAME="HALO_NMM_INIT_29.inc_73"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1962'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_30.inc.html'>HALO_NMM_INIT_30.inc</A>><A NAME="HALO_NMM_INIT_30.inc_74"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1963'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_31.inc.html'>HALO_NMM_INIT_31.inc</A>><A NAME="HALO_NMM_INIT_31.inc_75"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1964'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_32.inc.html'>HALO_NMM_INIT_32.inc</A>><A NAME="HALO_NMM_INIT_32.inc_76"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1965'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_33.inc.html'>HALO_NMM_INIT_33.inc</A>><A NAME="HALO_NMM_INIT_33.inc_77"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1966'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_34.inc.html'>HALO_NMM_INIT_34.inc</A>><A NAME="HALO_NMM_INIT_34.inc_78"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1967'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_35.inc.html'>HALO_NMM_INIT_35.inc</A>><A NAME="HALO_NMM_INIT_35.inc_79"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1968'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_36.inc.html'>HALO_NMM_INIT_36.inc</A>><A NAME="HALO_NMM_INIT_36.inc_80"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1969'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_37.inc.html'>HALO_NMM_INIT_37.inc</A>><A NAME="HALO_NMM_INIT_37.inc_81"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1970'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_38.inc.html'>HALO_NMM_INIT_38.inc</A>><A NAME="HALO_NMM_INIT_38.inc_82"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1971'>
# include <<A href='../../html_code/include/HALO_NMM_INIT_39.inc.html'>HALO_NMM_INIT_39.inc</A>><A NAME="HALO_NMM_INIT_39.inc_83"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1972'>
#endif<a name='1973'>
<a name='1974'>
#define COPY_OUT<a name='1975'>
#include <<A href='../../html_code/include/nmm_scalar_derefs.inc.html'>nmm_scalar_derefs.inc</A>><A NAME="nmm_scalar_derefs.inc_84"><A href='../../html_code/dyn_nmm/start_domain_nmm.F.html#START_DOMAIN_NMM' TARGET='bottom_target'><IMG SRC="../../gif/cyan.gif" border=0></a><a name='1976'>
<a name='1977'>
RETURN<a name='1978'>
<a name='1979'>
<a name='1980'>
END SUBROUTINE start_domain_nmm<a name='1981'>
<a name='1982'>
</pre></body></html>