<p><b>laura@ucar.edu</b> 2013-01-15 16:00:59 -0700 (Tue, 15 Jan 2013)</p><p>added entire physics_gfs directory<br>
</p><hr noshade><pre><font color="gray">Deleted: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.ibm
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.ibm        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.ibm        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1,824 +0,0 @@
-####   make with mpxlf_r (threading) 
-#EXEC    =  f.applic.x   
-#F77     = mpxlf
-#
-#FINCS    = -I/global/save/wx23sm/BASELIBS/esmf_3_1_0rp2_src/esmf/mod/modO/AIX.default.64.mpi.default
-FINCS    = -I/nwprod/lib/incmod/esmf_3_1_0rp2
- FINCM    = -I/nwprod/lib/incmod/w3_d
-#FINCM    = -I/global/save/wx23sm/GFS/f2010/trunk/para/lib/incmod/w3_d
-
-ARCH    = -qarch=pwr6 -qtune=pwr6 -qcache=auto -qnohot
-PGSZ    = -bdatapsize:64K -bstackpsize:64K -btextpsize:4K
-
-OPTS     = -g -qsuffix=cpp=f -O3 -qrealsize=8 -qstrict -qxlf77=leadzero -qmaxmem=-1 -qnolm -qsmp=noauto -qnosave $(ARCH)
-OPTS90   = -qsuffix=cpp=f -O3 -qrealsize=8 -qstrict -qxlf77=leadzero -qmaxmem=-1 -qnolm -qsmp=noauto -qnosave $(ARCH)
-OPTS90A  = -qsuffix=cpp=f -O3 -qrealsize=8 -qstrict -qxlf77=leadzero -qmaxmem=-1 -qnolm -qsmp=noauto -qnosave $(ARCH)
-#TRAPS    = -qinitauto=7FBFFFFF -qflttrap=enable:overflow:zerodivide:invalid
-#TRAPS    = -qinitauto=7FBFFFFF -qflttrap=enable:overflow:zerodivide:invalid -qsigtrap=xl__trcedump
-#DEBUG    = -g -qfullpath -qcheck -qextchk
-#
-#FFLAGM  = -g -qsuffix=cpp=f -O0 -qrealsize=8 -qstrict -qxlf77=leadzero -qmaxmem=-1 -qnolm -qsmp=noauto -qnosave -q fixed -qinitauto=7FBFFFFF -qflttrap=enable:overflow:zerodivide:invalid -qsigtrap=xl__trcedump -qcheck -qextchk -NS2048
-#FFLAGM2 = -g $(FINCM) -qsuffix=cpp=f -O0 -qrealsize=8 -qstrict -qxlf77=leadzero -qmaxmem=-1 -qnolm -qsmp=noauto -qnosave -q fixed -qinitauto=7FBFFFFF -qflttrap=enable:overflow:zerodivide:invalid -qsigtrap=xl__trcedump -qcheck -qextchk -NS2048
- FFLAG90  = $(OPTS90) $(FINCS) $(FINCM) -qfree -NS2048
- FFLAG90A = $(OPTS90A) $(FINCS) -qfree -NS2048
- FFLAGS   = $(OPTS) $(TRAPS) -qfixed
- FFLAGX   = $(OPTS) $(TRAPS) -qfixed
- FFLAGIO  = $(OPTS) $(TRAPS) -qfixed
- FFLAGY   = $(OPTS) -qfixed
- FFLAGM   = $(OPTS) $(FINCS) $(TRAPS) $(DEBUG) -NS2048 -qfixed
- FFLAGM2  = $(OPTS) $(FINCS) $(FINCM) $(TRAPS) $(DEBUG) -NS2048 -qfixed
- FFLAGM3  = $(OPTS) $(FINCS) $(TRAPS) $(DEBUG) -NS2048 -qfree
- FFLAGSF  = -g -O3 -qnosave -qfree=f90 -qcheck
- FFLAGSI  = -g -O3 -qnosave -qfree=f90
- FFLAGB   = -g -O3 -qnosave -qfixed
-
- LDR     = mpxlf95_r -qsmp=noauto 
-LDFLAGS = -lessl_r -lmass -qsmp=noauto
-
-#LIBS = -lC ${ESMFLIB} -L /nwprod/lib/ -l w3_d -l bacio_4
-#ESMFLIB  = /global/save/wx20wy/BASELIBS/esmf_3_1_0_src/AIX/lib/libO
-#LDFLAGS = -lessl_r -lmass -qsmp=noauto ${ESMFLIB}/libesmf.a ${ESMFLIB}/libnetcdf_stubs.a
-
-#ESMFLIB  = /global/save/wx23sm/BASELIBS/esmf_3_1_0rp2_src/esmf/lib/libO/AIX.default.64.mpi.default
-#LDFLAGS = -lessl_r -lmass -qsmp=noauto ${ESMFLIB}/libesmf.a $(PGSZ)
-ESMFLIB  = /nwprod/lib
-LDFLAGS = -lessl_r -lmass -qsmp=noauto ${ESMFLIB}/libesmf_3_1_0rp2.a $(PGSZ)
-LIBS = -lC -L /nwprod/lib/ -l w3_d -l bacio_4 -lsp_d
-#LIBS = -lC -L /global/save/wx23sm/GFS/f2010/trunk/para/lib -l bacio_4 -l w3_d -l sp_d
-
-.SUFFIXES: .o .f .F .h
-#
-# *****************************************************************
-#
-OBJS0        =                               \
-          NameList_ESMFMod.o          \
-          Sfc_Var_ESMFMod.o           \
-          Nst_Var_ESMFMod.o           \
-          GFS_ErrMsgMod.o             \
-          Lib_ESMFStateAddGetMod.o    \
-          GFS_InternalState_ESMFMod.o \
-          GFS_ESMFStateMod.o          \
-          GFS_GetCf_ESMFMod.o         \
-          GFS_Initialize_ESMFMod.o    \
-          GFS_Run_ESMFMod.o           \
-          GFS_Finalize_ESMFMod.o      \
-          GFS_AddParameterToStateMod.o \
-          GFS_GridComp_ESMFMod.o      \
-          Grid_ESMFCreate.o           \
-          StartTimeGet_ESMF.o         \
-          Ensemble_sub.o              \
-          GEFS_Cpl_InternalState_ESMFMod.o         \
-          GEFS_CplState_ESMFMod.o                  \
-          GEFS_four2fln_gg.o                        \
-          GEFS_uvoedz.o                            \
-          GEFS_uveodz.o                            \
-          GEFS_sumfln_slg_gg.o                     \
-          GEFS_dozeuv.o                            \
-          GEFS_dezouv.o                            \
-          GEFS_Sto_Per_Scheme_Step1.o              \
-          GEFS_Sto_Per_Scheme_Step2.o              \
-          GEFS_Cpl_Run_ESMFMod.o                   \
-          GEFS_Cpl_ESMFMod.o                       \
-          GEFS_GetParameterFromStateMod.o          \
-          GEFS_CplComp_ESMFMod.o                   \
-          GFS_ESMF.o                               \
-          Cal_Sto_Coef.o                           \
-          GEFS_bcst_global.o
-
-#          GFS_ESMFStateAddGetMod.o
-#          GFS_Standalone_ESMF_ENS.o
-#          GFS_InputFile2ImportState.o 
-
-
-OBJ_MOD        = machine.o             \
-          num_parthds.o         \
-          module.o              \
-          iounitdef.o           \
-          physcons.o            \
-          funcphys.o            \
-          progtm_module.o       \
-          rascnvv2.o            \
-          resol_def.o           \
-          gg_def.o              \
-          vert_def.o            \
-          sig_io.o              \
-          date_def.o            \
-          layout1.o             \
-          layout_grid_tracers.o \
-          namelist_def.o        \
-          namelist_soilveg.o    \
-          coordinate_def.o      \
-          tracer_const_h-new.o  \
-          reduce_lons_grid_module.o \
-          matrix_sig_def.o      \
-          mpi_def.o             \
-          semi_lag_def.o        \
-          deldifs_def.o         \
-          deldifs_hyb_slg_def.o       \
-          bfilt_def.o           \
-          sfcio_module.o        \
-          sigio_module.o        \
-          sigio_r_module.o      \
-          reduce_lons_grid_module.o \
-          gfsio_module.o        \
-          gfsio_def.o           \
-           d3d_def.o             \
-          nstio_module.o           \
-          module_nst_parameters.o  \
-          module_nst_water_prop.o  \
-          module_nst_model.o       \
-          calpreciptype.o          \
-          module_bfmicrophysics.o  \
-           kinds.o                  \
-           peuc.o                   \
-           pran.o                   \
-          prana.o                  \
-          layout_lag.o             \
-          slglat.o                 \
-          slgshr.o                 \
-          pmgrid.o                 \
-          set_longgs.o             \
-          get_cd_hyb_slg.o         \
-          sltini_all.o             \
-          set_halos_lag_sn.o       \
-          getcon_lag.o             \
-          gfidi_speed.o            \
-          plot85_a_h.o             \
-          setlats_lag.o            \
-          gfidi_tracers.o          \
-          slgscan_all_redgg.o      \
-          get_topo_grid_grad.o     \
-          mpi_tracers_a_to_b.o     \
-          gloopa_tracers.o         \
-          ndsl_pack.o                   \
-          ndsl_para.o                   \
-          sicdif_hyb_slg.o
-
-
-OBJS        = \
-        get_cd_hyb.o get_am_bm_hyb.o\
-        get_cd_hyb_gc.o get_am_bm_hyb_gc.o\
-        wrt3d_hyb.o\
-        wrtg3d_hyb.o\
-        setlats.o\
-        getysk.o\
-        getysk_gc.o\
-        ifshuff.o\
-        shuffle_grid.o\
-        digifilt.o\
-        gcycle.o\
-        wrtout.o\
-        grid_collect.o\
-        info.o\
-        four_to_grid.o\
-        hpmdummy.o\
-        coundummy.o\
-        damp_speed.o\
-        deldifs.o\
-        deldifs_tracers.o\
-        deldifs_hyb_ecm.o\
-        get_tc_ecm.o\
-        fix_fields.o\
-        spect_fields.o\
-        delnpe.o\
-        delnpo.o \
-        dezouv.o\
-        impadj_sig.locl.o\
-        impadj_slg.locl.o\
-        impadj_hyb.locl.o\
-        impadj_hyb.locl_gc.o\
-        dozeuv.o\
-        epslon.o\
-        filtr1eo.o \
-        filtr2eo.o\
-        get_ls_node.o\
-        get_lats_node_shuff_a_fix.o\
-        get_lats_node_shuff_r_fix.o\
-        getcon.o \
-        gfidi_sig.o\
-        gfidi_hyb.o\
-        gfidi_hyb_gc.o\
-        gfidi_hyb_gc_h-new.o\
-        gfidi_hyb_gc_ndsl.o\
-        ndsl_advecth.o\
-        glats.o\
-        gloopa_opt_gg.o\
-        lon_lat_para.o\
-        compns.o\
-        gozrineo.o\
-        pln2eo.o \
-        rms_spect.o \
-        uvoedz.o                     \
-        mpi_quit.o                   \
-        sicdif_sig.o                 \
-        sicdif_hyb.o                 \
-        sicdif_hyb_gc.o              \
-        dotstep_tracers.o            \
-        ensemble_wrt.o               \
-        gather_times.o               \
-        four2fln_gg.o                \
-        sortrx.o                     \
-        sumfln_slg_gg.o              \
-        sumder2_slg_gg.o             \
-        triseori.o                   \
-        uveodz.o                     \
-        ysminv.o                     \
-        ver_sig.o       
-
-
-##        wrtout_hyb.not.sure.o
-
-##        bartrieo.o\
-OBJS_PORT        = \
-        fftpack.o \
-        four2grid.fftpack.o \
-        noblas.o\
-        funcphys_subsx.o\
-
-OBJS_RAD        =              \
-        radlw_param.o          \
-        radlw_datatb.o         \
-        radlw_main.o           \
-        radsw_param.o          \
-        radsw_datatb.o         \
-        radsw_main.o           \
-        radiation_astronomy.o  \
-        radiation_aerosols.o   \
-        radiation_gases.o      \
-        radiation_clouds.o     \
-        radiation_surface.o    \
-        grrad.o                \
-        gloopr.o
-
-#
-#        astronomy.o    \
-##        funcphys_subs.o
-
-OBJS_PHY= \
-        ozinterp.o           \
-        ozphys.o             \
-        gloopb.o             \
-        gbphys.o             \
-        gbphys_adv_hyb_gc.o  \
-        gbphys_adv_hyb_gc_h-new.o \
-        dcyc2.o              \
-        dcyc2.pre.rad.o      \
-        set_soilveg.o        \
-        sfc_drv.o            \
-        sfc_land.o           \
-        progt2.o             \
-        sfc_sice.o           \
-        sfc_ocean.o          \
-        sfc_nst.o            \
-        sfc_diff.o           \
-        sfc_diag.o           \
-        sflx.o               \
-        moninp.o             \
-        moninp1.o            \
-        moninq.o             \
-        moninq1.o            \
-        tridi2t3.o           \
-        gwdps.o              \
-        gwdc.o               \
-        sascnv.o             \
-        sascnvn.o            \
-        cnvc90.o             \
-        shalcv.o             \
-        shalcv_opr.o         \
-        shalcnv.o            \
-        lrgsclr.o            \
-        gscond.o             \
-        precpd.o             \
-        mstadb.o             \
-        mstadbtn.o           \
-        mstcnv.o             \
-        get_prs.o            \
-        gsmddrive.o          \
-        omegtes.o            \
-        omegtes_gc.o         \
-        omegas.o             \
-        hyb2sig.o            \
-        hyb2press.o          \
-        hyb2press_gc.o       \
-        sig2press.o
-
-#        gbphys_call.o       \
-###funcphys_subsx.o   only srt gpxs was called in fix_fields -table not used
-
-OBJS_IO= \
-        read_fix.o            \
-        gribit.o              \
-        wrt3d.o               \
-        wrt3d_hyb.o           \
-        wrtg3d.o              \
-        wrtg3d_hyb.o          \
-        wrtsfc.o              \
-        para_fixio_w.o        \
-        para_nstio_w.o        \
-        sfcsub.o              \
-        treadeo.io.o          \
-        treadeo.gfsio.o       \
-        grid_to_spec.o        \
-        spect_to_grid.o       \
-        spect_tv_enthalpy_ps.o\
-        setsig.o              \
-        twriteeo.o            \
-        bafrio.o              \
-        spect_send.o          \
-        spect_write.o
-
-#        spect_tv_enthalpy_ps.o\   # To be used later
-
-#        para_fixio_all_w.o \
-#        para_fixio_iop_w.o \
-#        wrtsfc_all.o       \
-#        wrtsfc_iop.o       \
-##         getaer.o           \
-##         conrad.o           \
-##         crhtab.o           \
-
-OBJS_CC= \
-        mpi_more.o  \
-        cmp.comm.o  \
-        atm.comm.o  \
-        tiles.o
-
-SRC        = $(OBJS0:.o=.f) $(OBJ_MOD:.o=.f) $(OBJS:.o=.f) $(OBJS_RAD:.o=.f) $(OBJS_PHY:.o=.f) $(OBJS_IO:.o=.f) $(OBJS_CC:.o=.f)
-#
-INCS = f_hpm.h mpi_inc.h function2
-
-#
-# *****************************************************************
-#
-all: model-mpi
-
-model-mpi: $(OBJ_MOD) $(OBJS_CC) $(OBJS0) $(OBJS) $(OBJS_PHY) $(OBJS_RAD) $(OBJS_IO) 
-        $(LDR) $(LDFLAGS) -o $(EXEC) $(OBJ_MOD) $(OBJS_CC) $(OBJS0) $(OBJS) $(OBJS_PHY) $(OBJS_RAD) $(OBJS_IO) $(LIBS)
-
-clean:
-        rm -f $(OBJ_MOD) $(OBJS0) $(OBJS) $(OBJS_RAD) $(OBJS_PHY) $(OBJS_IO) *.mod
-
-tar:
-        tar -cvf tar.gfs.r4r8 $(SRC) $(INCS) $(COMS) $(OBJS_PORT:.o=.f) lonsper* res* xx* Makefile* ini.* scr.* m*real_?
-
-.F.o:
-        $(F77) $(FFLAGS) -c -d $&lt; 
-        #$(F77) $(FFLAGS) -c -d -WF,-DCLR:${RASV} $&lt; 
-.f.o:
-        $(F77) $(FFLAGS) -c $&lt; 
-
-
-omegas.o:        omegas.f
-                $(F77) $(FFLAGM) -c omegas.f
-
-#
-# *****************************************************************
-#
-cnvc90.o:        cnvc90.f
-                $(F77) $(FFLAGM) -c cnvc90.f
-
-calpreciptype.o:        calpreciptype.f
-                $(F77) $(FFLAGM3) -c calpreciptype.f
-
-dcyc2.o:        dcyc2.f
-                $(F77) $(FFLAGM) -c dcyc2.f
-
-dcyc2.pre.rad.o:        dcyc2.pre.rad.f
-                $(F77) $(FFLAGM) -c dcyc2.pre.rad.f
-
-digifilt.o:        digifilt.f
-                $(F77) $(FFLAGX) -c digifilt.f
-
-funcphys_subsx.o:        funcphys_subsx.f
-                $(F77) $(FFLAGM) -c funcphys_subsx.f
-
-gbphys_adv_hyb_gc.o:        gbphys_adv_hyb_gc.f
-                $(F77) $(FFLAGM) -c gbphys_adv_hyb_gc.f
-
-gbphys_adv_hyb_gc_h-new.o:        gbphys_adv_hyb_gc_h-new.f
-                $(F77) $(FFLAGM) -c gbphys_adv_hyb_gc_h-new.f
-
-#gbphys_call.o:        gbphys_call.f
-#                $(F77) $(FFLAGM) -c gbphys_call.f
-
-gbphys.o:        gbphys.f
-                $(F77) $(FFLAGM) -c gbphys.f
-
-get_prs.o:        get_prs.f
-                $(F77) $(FFLAGM) -c get_prs.f
-
-gscond.o:        gscond.f
-                $(F77) $(FFLAGM) -c gscond.f
-
-gsmddrive.o:        gsmddrive.f
-                $(F77) $(FFLAGM) -c gsmddrive.f
-
-reduce_lons_grid_module.o:        reduce_lons_grid_module.f
-                $(F77) $(FFLAGM) -c reduce_lons_grid_module.f
-
-module_bfmicrophysics.o:        module_bfmicrophysics.f
-                $(F77) $(FFLAGM) -c module_bfmicrophysics.f
-
-gwdps.o:        gwdps.f
-                $(F77) $(FFLAGM) -c gwdps.f
-
-gwdc.o:                gwdc.f
-                $(F77) $(FFLAGM) -c gwdc.f
-
-hyb2press_gc.o:        hyb2press_gc.f
-                $(F77) $(FFLAGM) -c hyb2press_gc.f
-
-hyb2press.o:        hyb2press.f
-                $(F77) $(FFLAGM) -c hyb2press.f
-
-hyb2sig.o:        hyb2sig.f
-                $(F77) $(FFLAGM) -c hyb2sig.f
-
-lrgsclr.o:        lrgsclr.f
-                $(F77) $(FFLAGM) -c lrgsclr.f
-
-moninp.o:        moninp.f
-                $(F77) $(FFLAGM) -c moninp.f
-
-moninp1.o:        moninp1.f
-                $(F77) $(FFLAGM) -c moninp1.f
-
-moninq.o:        moninq.f
-                $(F77) $(FFLAGM) -c moninq.f
-
-moninq1.o:        moninq1.f
-                $(F77) $(FFLAGM) -c moninq1.f
-
-mstadb.o:        mstadb.f
-                $(F77) $(FFLAGM) -c mstadb.f
-
-mstadbtn.o:        mstadbtn.f
-                $(F77) $(FFLAGM) -c mstadbtn.f
-
-mstcnv.o:        mstcnv.f
-                $(F77) $(FFLAGM) -c mstcnv.f
-
-omegtes_gc.o:        omegtes_gc.f
-                $(F77) $(FFLAGM) -c omegtes_gc.f
-
-omegtes.o:        omegtes.f
-                $(F77) $(FFLAGM) -c omegtes.f
-
-ozinterp.o:        ozinterp.f
-                $(F77) $(FFLAGM) -c ozinterp.f
-
-ozphys.o:        ozphys.f
-                $(F77) $(FFLAGM) -c ozphys.f
-
-precpd.o:        precpd.f
-                $(F77) $(FFLAGM) -c precpd.f
-
-rascnvv2.o:        rascnvv2.f
-                $(F77) $(FFLAGM) -c rascnvv2.f
-
-sfc_sice.o:        sfc_sice.f
-                $(F77) $(FFLAGM) -c sfc_sice.f
-
-set_soilveg.o:        set_soilveg.f
-                $(F77) $(FFLAGM) -c set_soilveg.f
-
-namelist_soilveg.o:        namelist_soilveg.f
-                $(F77) $(FFLAGM) -c namelist_soilveg.f
-
-sfc_land.o:        sfc_land.f
-                $(F77) $(FFLAGM) -c sfc_land.f
-                                                                            
-progt2.o:        progt2.f
-                $(F77) $(FFLAGM) -c progt2.f
-
-sfc_drv.o:        sfc_drv.f
-                $(F77) $(FFLAGM) -c sfc_drv.f
-
-sflx.o:         sflx.f
-                $(F77) $(FFLAGM) -c sflx.f
-
-sfc_ocean.o:        sfc_ocean.f
-                $(F77) $(FFLAGM) -c sfc_ocean.f
-
-sfc_nst.o:        sfc_nst.f
-                $(F77) $(FFLAGM) -c sfc_nst.f
-
-sfc_diff.o:        sfc_diff.f
-                $(F77) $(FFLAGM) -c sfc_diff.f
-
-sfc_diag.o:        sfc_diag.f
-                $(F77) $(FFLAGM) -c sfc_diag.f
-
-sascnv.o:        sascnv.f
-                $(F77) $(FFLAGM) -c sascnv.f
-
-sascnvn.o:        sascnvn.f
-                $(F77) $(FFLAGM) -c sascnvn.f
-
-tridi2t3.o:        tridi2t3.f
-                $(F77) $(FFLAGM) -c tridi2t3.f
-
-shalcv.o:        shalcv.f
-                $(F77) $(FFLAGM) -c shalcv.f
-
-shalcv_opr.o:        shalcv_opr.f
-                $(F77) $(FFLAGM) -c shalcv_opr.f
-
-shalcnv.o:        shalcnv.f
-                $(F77) $(FFLAGM) -c shalcnv.f
-
-sig2press.o:        sig2press.f
-                $(F77) $(FFLAGM) -c sig2press.f
-
-# *****************************************************************
-
-radlw_param.o:        radlw_param.f
-                $(F77) $(FFLAGM) -c radlw_param.f
-
-
-radlw_datatb.o:        radlw_datatb.f
-                $(F77) $(FFLAGM) -c radlw_datatb.f
-
-
-radlw_main.o:        radlw_main.f
-                $(F77) $(FFLAGM2) -c radlw_main.f
-
-radsw_param.o:        radsw_param.f
-                $(F77) $(FFLAGM) -c radsw_param.f
-
-radsw_datatb.o:        radsw_datatb.f
-                $(F77) $(FFLAGM) -c radsw_datatb.f
-
-radsw_main.o:        radsw_main.f
-                $(F77) $(FFLAGM2) -c radsw_main.f
-
-radiation_astronomy.o:        radiation_astronomy.f
-                $(F77) $(FFLAGM) -c radiation_astronomy.f
-
-radiation_aerosols.o:        radiation_aerosols.f
-                $(F77) $(FFLAGM) -c radiation_aerosols.f
-
-radiation_gases.o:        radiation_gases.f
-                $(F77) $(FFLAGM) -c radiation_gases.f
-
-radiation_clouds.o:        radiation_clouds.f
-                $(F77) $(FFLAGM) -c radiation_clouds.f
-
-radiation_surface.o:        radiation_surface.f
-                $(F77) $(FFLAGM) -c radiation_surface.f
-
-grrad.o:        grrad.f
-                $(F77) $(FFLAGM) -c grrad.f
-
-progtm_module.o:        progtm_module.f
-                $(F77) $(FFLAGM) -c progtm_module.f
-
-machine.o:        machine.f
-                $(F77) $(FFLAGM) -c machine.f
-
-num_parthds.o:        num_parthds.f
-                $(F77) $(FFLAGM) -c num_parthds.f
-
-kinds.o:        GEFS_Cpl_Cal_Sto_Coef.fd/kinds.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/kinds.f90
-
-peuc.o:                GEFS_Cpl_Cal_Sto_Coef.fd/peuc.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/peuc.f90
-
-pran.o:                GEFS_Cpl_Cal_Sto_Coef.fd/pran.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/pran.f90
-
-prana.o:        GEFS_Cpl_Cal_Sto_Coef.fd/prana.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/prana.f90
-
-#
-# *****************************************************************
-#
-physcons.o:        physcons.f
-                $(F77) $(FFLAG90) -c physcons.f
-
-iounitdef.o:        iounitdef.f
-                $(F77) $(FFLAG90) -c iounitdef.f
-
-funcphys.o:        funcphys.f
-                $(F77) $(FFLAG90) -c funcphys.f
-
-sfcio_module.o:        sfcio_module.f
-                $(F77) $(FFLAGSF) -c sfcio_module.f
-
-sigio_module.o:        sigio_module.f
-                $(F77) $(FFLAGSI) -c sigio_module.f
-
-gfsio_def.o:        gfsio_def.f
-                $(F77) $(FFLAGSI) -c gfsio_def.f
-
-gfsio_module.o:        gfsio_module.f
-                $(F77) $(FFLAGSI) -c gfsio_module.f
-
-sigio_r_module.o:        sigio_r_module.f
-                $(F77) $(FFLAGSI) -c sigio_r_module.f
-
-bafrio.o:        bafrio.f
-                $(F77) $(FFLAGB) -c bafrio.f
-
-#
-#sigdas.io.o:        sigdas.io.f
-##                $(F77) $(FFLAGIO) -c sigdas.io.f
-
-read_fix.o:        read_fix.f
-                $(F77) $(FFLAGX) -c read_fix.f
-
-softcount.o:        softcount.f
-                $(F77) $(FFLAGX) -c softcount.f
-
-gloopr.o:        gloopr.f
-                $(F77) $(FFLAGM2) -c gloopr.f
-
-gloopb.o:        gloopb.f
-                $(F77) $(FFLAGM2) -c gloopb.f
-
-sfcsub.o:        sfcsub.f
-                $(F77) $(FFLAGM) -c sfcsub.f
-
-gcycle.o:        gcycle.f
-                $(F77) $(FFLAGM) -c gcycle.f
-
-getaer.o:        getaer.f
-                $(F77) $(FFLAGIO) -c getaer.f
-
-wrt3d.o:        wrt3d.f
-                $(F77) $(FFLAGIO) -c wrt3d.f
-
-gribit.o:        gribit.f
-                $(F77) $(FFLAGIO) -c gribit.f
-
-wrtsfc.o:        wrtsfc.f
-                $(F77) $(FFLAGIO) -c wrtsfc.f
-
-para_fixio_w.o:                para_fixio_w.f
-                $(F77) $(FFLAGIO) -c para_fixio_w.f
-
-para_nstio_w.o:                para_nstio_w.f
-                $(F77) $(FFLAGIO) -c para_nstio_w.f
-
-#para_fixio_all_w.o:                para_fixio_all_w.f
-#                $(F77) $(FFLAGIO) -c para_fixio_all_w.f
-
-#para_fixio_iop_w.o:        para_fixio_iop_w.f
-#                $(F77) $(FFLAGIO) -c para_fixio_iop_w.f
-
-#conrad.o:        conrad.f
-#                $(F77) $(FFLAGIO) -c conrad.f
-
-#crhtab.o:        crhtab.f
-#                $(F77) $(FFLAGIO) -c crhtab.f
-
-treadeo.io.o:        treadeo.io.f
-                $(F77) $(FFLAGIO) -c treadeo.io.f
-
-treadeo.gfsio.o:        treadeo.gfsio.f
-                $(F77) $(FFLAGIO) -c treadeo.gfsio.f
-
-twriteeo.o:        twriteeo.f
-                $(F77) $(FFLAGIO) -c twriteeo.f
-
-spect_send.o:        spect_send.f
-                $(F77) $(FFLAGIO) -c spect_send.f
-
-spect_write.o:        spect_write.f
-                $(F77) $(FFLAGIO) -c spect_write.f
-
-spect_to_grid.o:        spect_to_grid.f
-                $(F77) $(FFLAGIO) -c spect_to_grid.f
-
-spect_tv_enthalpy_ps.o:        spect_tv_enthalpy_ps.f
-                $(F77) $(FFLAGIO) -c spect_tv_enthalpy_ps.f
-
-grid_to_spec.o:        grid_to_spec.f
-                $(F77) $(FFLAGIO) -c grid_to_spec.f
-
-wrtout.o:        wrtout.f
-                $(F77) $(FFLAGIO) -c wrtout.f
-
-nstio_module.o:         nstio_module.f
-                $(F77) $(FFLAG90) -c nstio_module.f
-
-module_nst_parameters.o:        module_nst_parameters.f
-                $(F77) $(FFLAG90) -c module_nst_parameters.f
-
-module_nst_water_prop.o:         module_nst_water_prop.f
-                $(F77) $(FFLAG90) -c module_nst_water_prop.f
-
-module_nst_model.o:                module_nst_model.f
-                $(F77) $(FFLAG90) -c module_nst_model.f
-
-NameList_ESMFMod.o:         NameList_ESMFMod.f
-                $(F77) $(FFLAG90) -c NameList_ESMFMod.f
-
-Sfc_Var_ESMFMod.o:        Sfc_Var_ESMFMod.f
-                $(F77) $(FFLAG90) -c Sfc_Var_ESMFMod.f
-
-Nst_Var_ESMFMod.o:        Nst_Var_ESMFMod.f
-                $(F77) $(FFLAG90) -c Nst_Var_ESMFMod.f
-
-Lib_ESMFStateAddGetMod.o:        Lib_ESMFStateAddGetMod.f
-                $(F77) $(FFLAG90) -c Lib_ESMFStateAddGetMod.f
-
-#GFS_ESMFStateAddGetMod.o:         GFS_ESMFStateAddGetMod.f
-#                $(F77) $(FFLAG90) -c GFS_ESMFStateAddGetMod.f
-
-GFS_InternalState_ESMFMod.o:         GFS_InternalState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_InternalState_ESMFMod.f
-
-GFS_ESMFStateMod.o:         GFS_ESMFStateMod.f
-                $(F77) $(FFLAG90) -c GFS_ESMFStateMod.f
-
-GFS_ErrMsgMod.o:         GFS_ErrMsgMod.f
-                $(F77) $(FFLAG90) -c GFS_ErrMsgMod.f
-
-GFS_GetCf_ESMFMod.o:         GFS_GetCf_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_GetCf_ESMFMod.f
-
-GFS_ESMFMod.o:         GFS_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_ESMFMod.f
-
-#GFS_Grid_fnl_ESMFMod.o:         GFS_Grid_fnl_ESMFMod.f
-        #$(F77) $(FFLAG90A) -c GFS_Grid_fnl_ESMFMod.f
-
-GFS_GridComp_ESMFMod.o:         GFS_GridComp_ESMFMod.f
-                $(F77) $(FFLAG90A) -c GFS_GridComp_ESMFMod.f
-
-GFS_Initialize_ESMFMod.o:         GFS_Initialize_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Initialize_ESMFMod.f
-
-GFS_Run_ESMFMod.o:        GFS_Run_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Run_ESMFMod.f
-
-GFS_Finalize_ESMFMod.o:        GFS_Finalize_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Finalize_ESMFMod.f
-
-GFS_InputFile2ImportState.o:         GFS_InputFile2ImportState.f
-                $(F77) $(FFLAG90) -c GFS_InputFile2ImportState.f
-
-#GFS_Standalone_ESMF_ENS.o:         GFS_Standalone_ESMF_ENS.f
-#                $(F77) $(FFLAG90) -c GFS_Standalone_ESMF_ENS.f
-
-GFS_ESMF.o:        GFS_ESMF.f
-                $(F77) $(FFLAG90) -c GFS_ESMF.f
-
-Grid_ESMFCreate.o:         Grid_ESMFCreate.f
-                $(F77) $(FFLAG90) -c Grid_ESMFCreate.f
-
-StartTimeGet_ESMF.o:         StartTimeGet_ESMF.f
-                $(F77) $(FFLAG90) -c StartTimeGet_ESMF.f
-
-Ensemble_sub.o:         Ensemble_sub.f
-                $(F77) $(FFLAG90) -c Ensemble_sub.f
-
-mpi_more.o:        mpi_more.f
-                $(F77) $(FFLAGS) -c mpi_more.f
-
-cmp.comm.o:        cmp.comm.f
-                $(F77) $(FFLAGS) -c cmp.comm.f
-
-atm.comm.o:        atm.comm.f
-                $(F77) $(FFLAGS) -c atm.comm.f
-
-tiles.o:        tiles.f
-                $(F77) $(FFLAGS) -c tiles.f
-
-GEFS_Cpl_InternalState_ESMFMod.o:        GEFS_Cpl_InternalState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_InternalState_ESMFMod.f
-
-GEFS_CplState_ESMFMod.o:        GEFS_CplState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_CplState_ESMFMod.f
-
-GEFS_Sto_Per_Scheme.o:        GEFS_Sto_Per_Scheme.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme.f
-
-GEFS_Cpl_Run_ESMFMod.o:        GEFS_Cpl_Run_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_Run_ESMFMod.f
-
-GEFS_Cpl_ESMFMod.o:        GEFS_Cpl_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_ESMFMod.f
-
-GEFS_CplComp_ESMFMod.o:        GEFS_CplComp_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_CplComp_ESMFMod.f
-
-GFS_AddParameterToStateMod.o:        GFS_AddParameterToStateMod.f
-                $(F77) $(FFLAG90) -c GFS_AddParameterToStateMod.f
-
-GEFS_Sto_Per_Scheme_Step1.o:        GEFS_Sto_Per_Scheme_Step1.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme_Step1.f
-
-GEFS_Sto_Per_Scheme_Step2.o:        GEFS_Sto_Per_Scheme_Step2.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme_Step2.f
-
-GEFS_GetParameterFromStateMod.o:        GEFS_GetParameterFromStateMod.f
-                $(F77) $(FFLAG90) -c GEFS_GetParameterFromStateMod.f
-
-Cal_Sto_Coef.o:                GEFS_Cpl_Cal_Sto_Coef.fd/Cal_Sto_Coef.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/Cal_Sto_Coef.f90
-
-GEFS_bcst_global.o:        GEFS_bcst_global.f
-                $(F77) $(FFLAG90) -c GEFS_bcst_global.f

Deleted: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.jet
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.jet        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/Makefile_GFS.jet        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1,829 +0,0 @@
-
-EXEC = global_fcst
-F77= mpif90
-F90= mpif90
-LIBDIR=/mnt/pan2/projects/emcda/dkleist/nwprod/lib
-
-#
- FINCS    = -I /mnt/pan2/projects/emcda/dkleist/nwprod/incmod/esmf_3_1_0rp2
- FINCM    = -I /mnt/pan2/projects/emcda/dkleist/nwprod/incmod/w3lib-2.0_d
-## FINCG   = -I /mnt/pan2/projects/emcda/dkleist/nwprod/incmod/sigio_4
-## FINCC   = -I /mnt/pan2/projects/emcda/dkleist/nwprod/incmod/sfcio_4
-
- ARCH    = 
- PGSZ    = 
-
- OPTS     = -O3 -convert big_endian -traceback -r8
- OPTS90   = -O3 -convert big_endian -traceback -r8
- OPTS90A  = -O3 -convert big_endian -traceback -r8
-
- FFLAG90  = $(OPTS90) $(FINCS) -free
- FFLAG90A = $(OPTS90A) $(FINCS) -free
- FFLAGS   = $(OPTS) $(TRAPS) 
- FFLAGX   = $(OPTS) $(TRAPS) 
- FFLAGIO  = $(OPTS) $(TRAPS)
- FFLAGY   = $(OPTS)
- FFLAGM   = $(OPTS) $(FINCS) $(TRAPS) $(DEBUG)
- FFLAGM2  = $(OPTS) $(FINCS) $(FINCM) $(TRAPS) $(DEBUG)
- FFLAGM3  = $(OPTS) $(FINCS) $(TRAPS) $(DEBUG) -free
- FFLAG_SER = -O3 -convert big_endian   -traceback -r8
-
- FFLAGSF  = -O3 -convert big_endian  -traceback -FR
- FFLAGSI  = -O3 -convert big_endian  -traceback -FR
- FFLAGB   = -O3 -convert big_endian  -traceback
-
- LDR     = mpif90
-
- LDFLAGS = 
- LIBDIR = /mnt/pan2/projects/emcda/dkleist/nwprod/lib
-##LIBS = -L$(MKL) -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lstdc++ -limf -lm -lrt -ldl -L$(LIBDIR) -lsp_d -lw3lib-2.0_d -lbacio_4 -lesmf_3_1_0rp2 -threads
-##LIBS = -lstdc++ -limf -lm -lrt -ldl -threads -L$(LIBDIR) -lsp_d -lw3lib-2.0_d -lbacio_4 -lesmf_3_1_0rp2 -lsfcio_4 -lsigio_4 -L${MPICH}/lib -lmpichcxx -mkl=sequential \
-##      -L$(MKL) -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lguide
-
-LIBS = -L$(MKL) -lmkl_intel_lp64 -lmkl_core -lmkl_intel_thread -lguide -L$(LIBDIR) -lesmf_3_1_0rp2 -lbacio_4 -lsp_d -lw3lib-2.0_d -lrt -lstdc++
-
-
-.SUFFIXES: .o .f .F .h
-#
-# *****************************************************************
-#
-OBJS0        =                               \
-          NameList_ESMFMod.o          \
-          Sfc_Var_ESMFMod.o           \
-          Nst_Var_ESMFMod.o           \
-          GFS_ErrMsgMod.o             \
-          Lib_ESMFStateAddGetMod.o    \
-          GFS_InternalState_ESMFMod.o \
-          GFS_ESMFStateMod.o          \
-          GFS_GetCf_ESMFMod.o         \
-          GFS_Initialize_ESMFMod.o    \
-          GFS_Run_ESMFMod.o           \
-          GFS_Finalize_ESMFMod.o      \
-          GFS_AddParameterToStateMod.o \
-          GFS_GridComp_ESMFMod.o      \
-          Grid_ESMFCreate.o           \
-          StartTimeGet_ESMF.o         \
-          Ensemble_sub.o              \
-          GEFS_Cpl_InternalState_ESMFMod.o         \
-          GEFS_CplState_ESMFMod.o                  \
-          GEFS_four2fln_gg.o                        \
-          GEFS_uvoedz.o                            \
-          GEFS_uveodz.o                            \
-          GEFS_sumfln_slg_gg.o                     \
-          GEFS_dozeuv.o                            \
-          GEFS_dezouv.o                            \
-          GEFS_Sto_Per_Scheme_Step1.o              \
-          GEFS_Sto_Per_Scheme_Step2.o              \
-          GEFS_Cpl_Run_ESMFMod.o                   \
-          GEFS_Cpl_ESMFMod.o                       \
-          GEFS_GetParameterFromStateMod.o          \
-          GEFS_CplComp_ESMFMod.o                   \
-          GFS_ESMF.o                               \
-          Cal_Sto_Coef.o                           \
-          GEFS_bcst_global.o
-
-#          GFS_ESMFStateAddGetMod.o
-#          GFS_Standalone_ESMF_ENS.o
-#          GFS_InputFile2ImportState.o 
-
-
-OBJ_MOD        = machine.o             \
-          num_parthds.o         \
-          module.o              \
-          iounitdef.o           \
-          physcons.o            \
-          funcphys.o            \
-          progtm_module.o       \
-          rascnvv2.o            \
-          mersenne_twister.o    \
-          resol_def.o           \
-          gg_def.o              \
-          vert_def.o            \
-          sig_io.o              \
-          date_def.o            \
-          layout1.o             \
-          layout_grid_tracers.o \
-          namelist_def.o        \
-          namelist_soilveg.o    \
-          coordinate_def.o      \
-          tracer_const_h-new.o  \
-          reduce_lons_grid_module.o \
-          matrix_sig_def.o      \
-          mpi_def.o             \
-          semi_lag_def.o        \
-          deldifs_def.o         \
-          deldifs_hyb_slg_def.o       \
-          bfilt_def.o           \
-          reduce_lons_grid_module.o \
-          gfsio_module.o        \
-          gfsio_def.o           \
-          sigio_module.o        \
-          sigio_r_module.o      \
-          sfcio_module.o        \
-          d3d_def.o             \
-          nstio_module.o           \
-          module_nst_parameters.o  \
-          module_nst_water_prop.o  \
-          module_nst_model.o       \
-          calpreciptype.o          \
-          module_bfmicrophysics.o  \
-          kinds.o                  \
-          peuc.o                   \
-          pran.o                   \
-          prana.o                  \
-          layout_lag.o             \
-          slglat.o                 \
-          slgshr.o                 \
-          pmgrid.o                 \
-          set_longgs.o             \
-          get_cd_hyb_slg.o         \
-          sltini_all.o             \
-          set_halos_lag_sn.o       \
-          getcon_lag.o             \
-          gfidi_speed.o            \
-          plot85_a_h.o             \
-          setlats_lag.o            \
-          gfidi_tracers.o          \
-          slgscan_all_redgg.o      \
-          get_topo_grid_grad.o     \
-          mpi_tracers_a_to_b.o     \
-          gloopa_tracers.o         \
-          ndsl_pack.o                   \
-          ndsl_para.o                   \
-          sicdif_hyb_slg.o
-
-
-OBJS        = \
-        get_cd_hyb.o get_am_bm_hyb.o\
-        get_cd_hyb_gc.o get_am_bm_hyb_gc.o\
-        wnew.o\
-        wrt3d_hyb.o\
-        wrtg3d_hyb.o\
-        setlats.o\
-        getysk.o\
-        getysk_gc.o\
-        ifshuff.o\
-        shuffle_grid.o\
-        digifilt.o\
-        gcycle.o\
-        wrtout.o\
-        grid_collect.o\
-        info.o\
-        four_to_grid.o\
-        hpmdummy.o\
-        coundummy.o\
-        damp_speed.o\
-        deldifs.o\
-        deldifs_tracers.o\
-        deldifs_hyb_ecm.o\
-        get_tc_ecm.o\
-        fix_fields.o\
-        spect_fields.o\
-        delnpe.o\
-        delnpo.o \
-        dezouv.o\
-        impadj_sig.locl.o\
-        impadj_slg.locl.o\
-        impadj_hyb.locl.o\
-        impadj_hyb.locl_gc.o\
-        dozeuv.o\
-        epslon.o\
-        filtr1eo.o \
-        filtr2eo.o\
-        get_ls_node.o\
-        get_lats_node_shuff_a_fix.o\
-        get_lats_node_shuff_r_fix.o\
-        getcon.o \
-        gfidi_sig.o\
-        gfidi_hyb.o\
-        gfidi_hyb_gc.o\
-        gfidi_hyb_gc_h-new.o\
-        gfidi_hyb_gc_ndsl.o\
-        ndsl_advecth.o\
-        glats.o\
-        gloopa_opt_gg.o\
-        lon_lat_para.o\
-        compns.o\
-        gozrineo.o\
-        pln2eo.o \
-        rms_spect.o \
-        uvoedz.o                     \
-        mpi_quit.o                   \
-        sicdif_sig.o                 \
-        sicdif_hyb.o                 \
-        sicdif_hyb_gc.o              \
-        dotstep_tracers.o            \
-        ensemble_wrt.o               \
-        gather_times.o               \
-        four2fln_gg.o                \
-        sortrx.o                     \
-        sumfln_slg_gg.o              \
-        sumder2_slg_gg.o             \
-        triseori.o                   \
-        uveodz.o                     \
-        ysminv.o                     \
-        ver_sig.o       
-
-
-##        wrtout_hyb.not.sure.o
-
-##        bartrieo.o\
-OBJS_PORT        = \
-        fftpack.o \
-        four2grid.fftpack.o \
-        noblas.o\
-        funcphys_subsx.o\
-
-OBJS_RAD        =              \
-        radlw_param.o          \
-        radlw_datatb.o         \
-        radlw_main.o           \
-        radsw_param.o          \
-        radsw_datatb.o         \
-        radsw_main.o           \
-        radiation_astronomy.o  \
-        radiation_aerosols.o   \
-        radiation_gases.o      \
-        radiation_clouds.o     \
-        radiation_surface.o    \
-        grrad.o                \
-        gloopr.o
-
-#
-#        astronomy.o    \
-##        funcphys_subs.o
-
-OBJS_PHY= \
-        ozinterp.o           \
-        ozphys.o             \
-        gloopb.o             \
-        gbphys.o             \
-        gbphys_adv_hyb_gc.o  \
-        gbphys_adv_hyb_gc_h-new.o \
-        dcyc2.o              \
-        dcyc2.pre.rad.o      \
-        set_soilveg.o        \
-        sfc_drv.o            \
-        sfc_land.o           \
-        progt2.o             \
-        sfc_sice.o           \
-        sfc_ocean.o          \
-        sfc_nst.o            \
-        sfc_diff.o           \
-        sfc_diag.o           \
-        sflx.o               \
-        moninp.o             \
-        moninp1.o            \
-        moninq.o             \
-        moninq1.o            \
-        tridi2t3.o           \
-        gwdps.o              \
-        gwdc.o               \
-        sascnv.o             \
-        sascnvn.o            \
-        cnvc90.o             \
-        shalcv.o             \
-        shalcv_opr.o         \
-        shalcnv.o            \
-        lrgsclr.o            \
-        gscond.o             \
-        precpd.o             \
-        mstadb.o             \
-        mstadbtn.o           \
-        mstcnv.o             \
-        get_prs.o            \
-        gsmddrive.o          \
-        omegtes.o            \
-        omegtes_gc.o         \
-        omegas.o             \
-        hyb2sig.o            \
-        hyb2press.o          \
-        hyb2press_gc.o       \
-        sig2press.o
-
-#        gbphys_call.o       \
-###funcphys_subsx.o   only srt gpxs was called in fix_fields -table not used
-
-OBJS_IO= \
-        read_fix.o            \
-        gribit.o              \
-        wrt3d.o               \
-        wrt3d_hyb.o           \
-        wrtg3d.o              \
-        wrtg3d_hyb.o          \
-        wrtsfc.o              \
-        para_fixio_w.o        \
-        para_nstio_w.o        \
-        sfcsub.o              \
-        treadeo.io.o          \
-        treadeo.gfsio.o       \
-        grid_to_spec.o        \
-        spect_to_grid.o       \
-        spect_tv_enthalpy_ps.o\
-        setsig.o              \
-        twriteeo.o            \
-        bafrio.o              \
-        spect_send.o          \
-        spect_write.o
-
-#        spect_tv_enthalpy_ps.o\   # To be used later
-
-#        para_fixio_all_w.o \
-#        para_fixio_iop_w.o \
-#        wrtsfc_all.o       \
-#        wrtsfc_iop.o       \
-##         getaer.o           \
-##         conrad.o           \
-##         crhtab.o           \
-
-OBJS_CC= \
-        mpi_more.o  \
-        cmp.comm.o  \
-        atm.comm.o  \
-        tiles.o
-
-SRC        = $(OBJS0:.o=.f) $(OBJ_MOD:.o=.f) $(OBJS:.o=.f) $(OBJS_RAD:.o=.f) $(OBJS_PHY:.o=.f) $(OBJS_IO:.o=.f) $(OBJS_CC:.o=.f)
-#
-INCS = f_hpm.h mpi_inc.h function2
-
-#
-# *****************************************************************
-#
-all: model-mpi
-
-model-mpi: $(OBJ_MOD) $(OBJS_CC) $(OBJS0) $(OBJS) $(OBJS_PHY) $(OBJS_RAD) $(OBJS_IO) 
-        $(LDR) -o $(EXEC) $(OBJ_MOD) $(OBJS_CC) $(OBJS0) $(OBJS) $(OBJS_PHY) $(OBJS_RAD) $(OBJS_IO) $(LIBS) $(LDFLAGS)
-
-clean:
-        rm -f $(OBJ_MOD) $(OBJS0) $(OBJS) $(OBJS_RAD) $(OBJS_PHY) $(OBJS_IO) *.mod
-
-tar:
-        tar -cvf tar.gfs.r4r8 $(SRC) $(INCS) $(COMS) $(OBJS_PORT:.o=.f) lonsper* res* xx* Makefile* ini.* scr.* m*real_?
-
-.F.o:
-        $(F77) $(FFLAGS) -c -d $&lt; 
-        #$(F77) $(FFLAGS) -c -d -WF,-DCLR:${RASV} $&lt; 
-.f.o:
-        $(F77) $(FFLAGS) -c $&lt; 
-
-
-omegas.o:        omegas.f
-                $(F77) $(FFLAGM) -c omegas.f
-
-#
-# *****************************************************************
-#
-cnvc90.o:        cnvc90.f
-                $(F77) $(FFLAGM) -c cnvc90.f
-
-calpreciptype.o:        calpreciptype.f
-                $(F77) $(FFLAGM3) -c calpreciptype.f
-
-dcyc2.o:        dcyc2.f
-                $(F77) $(FFLAGM) -c dcyc2.f
-
-dcyc2.pre.rad.o:        dcyc2.pre.rad.f
-                $(F77) $(FFLAGM) -c dcyc2.pre.rad.f
-
-digifilt.o:        digifilt.f
-                $(F77) $(FFLAGX) -c digifilt.f
-
-funcphys_subsx.o:        funcphys_subsx.f
-                $(F77) $(FFLAGM) -c funcphys_subsx.f
-
-gbphys_adv_hyb_gc.o:        gbphys_adv_hyb_gc.f
-                $(F77) $(FFLAGM) -c gbphys_adv_hyb_gc.f
-
-gbphys_adv_hyb_gc_h-new.o:        gbphys_adv_hyb_gc_h-new.f
-                $(F77) $(FFLAGM) -c gbphys_adv_hyb_gc_h-new.f
-
-#gbphys_call.o:        gbphys_call.f
-#                $(F77) $(FFLAGM) -c gbphys_call.f
-
-gbphys.o:        gbphys.f
-                $(F77) $(FFLAGM) -c gbphys.f
-
-get_prs.o:        get_prs.f
-                $(F77) $(FFLAGM) -c get_prs.f
-
-filtr1eo.o:     filtr1eo.f
-                $(F77) $(FFLAG_SER)   -c filtr1eo.f
-
-filtr2eo.o:     filtr2eo.f
-                $(F77) $(FFLAG_SER)   -c filtr2eo.f
-
-gscond.o:        gscond.f
-                $(F77) $(FFLAGM) -c gscond.f
-
-gsmddrive.o:        gsmddrive.f
-                $(F77) $(FFLAGM) -c gsmddrive.f
-
-reduce_lons_grid_module.o:        reduce_lons_grid_module.f
-                $(F77) $(FFLAGM) -c reduce_lons_grid_module.f
-
-module_bfmicrophysics.o:        module_bfmicrophysics.f
-                $(F77) $(FFLAGM) -c module_bfmicrophysics.f
-
-gwdps.o:        gwdps.f
-                $(F77) $(FFLAGM) -c gwdps.f
-
-gwdc.o:                gwdc.f
-                $(F77) $(FFLAGM) -c gwdc.f
-
-hyb2press_gc.o:        hyb2press_gc.f
-                $(F77) $(FFLAGM) -c hyb2press_gc.f
-
-hyb2press.o:        hyb2press.f
-                $(F77) $(FFLAGM) -c hyb2press.f
-
-hyb2sig.o:        hyb2sig.f
-                $(F77) $(FFLAGM) -c hyb2sig.f
-
-lrgsclr.o:        lrgsclr.f
-                $(F77) $(FFLAGM) -c lrgsclr.f
-
-moninp.o:        moninp.f
-                $(F77) $(FFLAGM) -c moninp.f
-
-moninp1.o:        moninp1.f
-                $(F77) $(FFLAGM) -c moninp1.f
-
-moninq.o:        moninq.f
-                $(F77) $(FFLAGM) -c moninq.f
-
-moninq1.o:        moninq1.f
-                $(F77) $(FFLAGM) -c moninq1.f
-
-mstadb.o:        mstadb.f
-                $(F77) $(FFLAGM) -c mstadb.f
-
-mstadbtn.o:        mstadbtn.f
-                $(F77) $(FFLAGM) -c mstadbtn.f
-
-mstcnv.o:        mstcnv.f
-                $(F77) $(FFLAGM) -c mstcnv.f
-
-omegtes_gc.o:        omegtes_gc.f
-                $(F77) $(FFLAGM) -c omegtes_gc.f
-
-omegtes.o:        omegtes.f
-                $(F77) $(FFLAGM) -c omegtes.f
-
-ozinterp.o:        ozinterp.f
-                $(F77) $(FFLAGM) -c ozinterp.f
-
-ozphys.o:        ozphys.f
-                $(F77) $(FFLAGM) -c ozphys.f
-
-precpd.o:        precpd.f
-                $(F77) $(FFLAGM) -c precpd.f
-
-rascnvv2.o:        rascnvv2.f
-                $(F77) $(FFLAGM) -c rascnvv2.f
-
-sfc_sice.o:        sfc_sice.f
-                $(F77) $(FFLAGM) -c sfc_sice.f
-
-set_soilveg.o:        set_soilveg.f
-                $(F77) $(FFLAGM) -c set_soilveg.f
-
-namelist_soilveg.o:        namelist_soilveg.f
-                $(F77) $(FFLAGM) -c namelist_soilveg.f
-
-sfc_land.o:        sfc_land.f
-                $(F77) $(FFLAGM) -c sfc_land.f
-                                                                            
-progt2.o:        progt2.f
-                $(F77) $(FFLAGM) -c progt2.f
-
-sfc_drv.o:        sfc_drv.f
-                $(F77) $(FFLAGM) -c sfc_drv.f
-
-sflx.o:         sflx.f
-                $(F77) $(FFLAGM) -c sflx.f
-
-sfc_ocean.o:        sfc_ocean.f
-                $(F77) $(FFLAGM) -c sfc_ocean.f
-
-sfc_nst.o:        sfc_nst.f
-                $(F77) $(FFLAGM) -c sfc_nst.f
-
-sfc_diff.o:        sfc_diff.f
-                $(F77) $(FFLAGM) -c sfc_diff.f
-
-sfc_diag.o:        sfc_diag.f
-                $(F77) $(FFLAGM) -c sfc_diag.f
-
-sascnv.o:        sascnv.f
-                $(F77) $(FFLAGM) -c sascnv.f
-
-sascnvn.o:        sascnvn.f
-                $(F77) $(FFLAGM) -c sascnvn.f
-
-tridi2t3.o:        tridi2t3.f
-                $(F77) $(FFLAGM) -c tridi2t3.f
-
-shalcv.o:        shalcv.f
-                $(F77) $(FFLAGM) -c shalcv.f
-
-shalcv_opr.o:        shalcv_opr.f
-                $(F77) $(FFLAGM) -c shalcv_opr.f
-
-shalcnv.o:        shalcnv.f
-                $(F77) $(FFLAGM) -c shalcnv.f
-
-sig2press.o:        sig2press.f
-                $(F77) $(FFLAGM) -c sig2press.f
-
-# *****************************************************************
-
-radlw_param.o:        radlw_param.f
-                $(F77) $(FFLAGM) -c radlw_param.f
-
-
-radlw_datatb.o:        radlw_datatb.f
-                $(F77) $(FFLAGM) -c radlw_datatb.f
-
-
-radlw_main.o:        radlw_main.f
-                $(F77) $(FFLAGM2) -c radlw_main.f
-
-radsw_param.o:        radsw_param.f
-                $(F77) $(FFLAGM) -c radsw_param.f
-
-radsw_datatb.o:        radsw_datatb.f
-                $(F77) $(FFLAGM) -c radsw_datatb.f
-
-radsw_main.o:        radsw_main.f
-                $(F77) $(FFLAGM2) -c radsw_main.f
-
-radiation_astronomy.o:        radiation_astronomy.f
-                $(F77) $(FFLAGM) -c radiation_astronomy.f
-
-radiation_aerosols.o:        radiation_aerosols.f
-                $(F77) $(FFLAGM) -c radiation_aerosols.f
-
-radiation_gases.o:        radiation_gases.f
-                $(F77) $(FFLAGM) -c radiation_gases.f
-
-radiation_clouds.o:        radiation_clouds.f
-                $(F77) $(FFLAGM) -c radiation_clouds.f
-
-radiation_surface.o:        radiation_surface.f
-                $(F77) $(FFLAGM) -c radiation_surface.f
-
-grrad.o:        grrad.f
-                $(F77) $(FFLAGM) -c grrad.f
-
-progtm_module.o:        progtm_module.f
-                $(F77) $(FFLAGM) -c progtm_module.f
-
-machine.o:        machine.f
-                $(F77) $(FFLAGM) -c machine.f
-
-num_parthds.o:        num_parthds.f
-                $(F77) $(FFLAGM) -c num_parthds.f
-
-kinds.o:        GEFS_Cpl_Cal_Sto_Coef.fd/kinds.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/kinds.f90
-
-peuc.o:                GEFS_Cpl_Cal_Sto_Coef.fd/peuc.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/peuc.f90
-
-pran.o:                GEFS_Cpl_Cal_Sto_Coef.fd/pran.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/pran.f90
-
-prana.o:        GEFS_Cpl_Cal_Sto_Coef.fd/prana.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/prana.f90
-
-#
-# *****************************************************************
-#
-physcons.o:        physcons.f
-                $(F77) $(FFLAG90) -c physcons.f
-
-iounitdef.o:        iounitdef.f
-                $(F77) $(FFLAG90) -c iounitdef.f
-
-funcphys.o:        funcphys.f
-                $(F77) $(FFLAG90) -c funcphys.f
-
-sfcio_module.o:        sfcio_module.f
-                $(F77) $(FFLAGSF) -c sfcio_module.f
-
-sigio_module.o:        sigio_module.f
-                $(F77) $(FFLAGSI) -c sigio_module.f
-
-gfsio_def.o:        gfsio_def.f
-                $(F77) $(FFLAGSI) -c gfsio_def.f
-
-gfsio_module.o:        gfsio_module.f
-                $(F77) $(FFLAGSI) -c gfsio_module.f
-
-sigio_r_module.o:        sigio_r_module.f
-                $(F77) $(FFLAGSI) -c sigio_r_module.f
-
-bafrio.o:        bafrio.f
-                $(F77) $(FFLAGB) -c bafrio.f
-
-#
-#sigdas.io.o:        sigdas.io.f
-##                $(F77) $(FFLAGIO) -c sigdas.io.f
-
-read_fix.o:        read_fix.f
-                $(F77) $(FFLAGX) -c read_fix.f
-
-softcount.o:        softcount.f
-                $(F77) $(FFLAGX) -c softcount.f
-
-gloopr.o:        gloopr.f
-                $(F77) $(FFLAG_SER) -c gloopr.f
-
-gloopb.o:        gloopb.f
-                $(F77) $(FFLAGM2) -c gloopb.f
-
-sfcsub.o:        sfcsub.f
-                $(F77) $(FFLAG_SER) -c sfcsub.f
-
-gcycle.o:        gcycle.f
-                $(F77) $(FFLAGM) -c gcycle.f
-
-getaer.o:        getaer.f
-                $(F77) $(FFLAGIO) -c getaer.f
-
-wrt3d.o:        wrt3d.f
-                $(F77) $(FFLAGIO) -c wrt3d.f
-
-gribit.o:        gribit.f
-                $(F77) $(FFLAGIO) -c gribit.f
-
-wrtsfc.o:        wrtsfc.f
-                $(F77) $(FFLAGIO) -c wrtsfc.f
-
-para_fixio_w.o:                para_fixio_w.f
-                $(F77) $(FFLAGIO) -c para_fixio_w.f
-
-para_nstio_w.o:                para_nstio_w.f
-                $(F77) $(FFLAGIO) -c para_nstio_w.f
-
-#para_fixio_all_w.o:                para_fixio_all_w.f
-#                $(F77) $(FFLAGIO) -c para_fixio_all_w.f
-
-#para_fixio_iop_w.o:        para_fixio_iop_w.f
-#                $(F77) $(FFLAGIO) -c para_fixio_iop_w.f
-
-#conrad.o:        conrad.f
-#                $(F77) $(FFLAGIO) -c conrad.f
-
-#crhtab.o:        crhtab.f
-#                $(F77) $(FFLAGIO) -c crhtab.f
-
-treadeo.io.o:        treadeo.io.f
-                $(F77) $(FFLAGIO) -c treadeo.io.f
-
-treadeo.gfsio.o:        treadeo.gfsio.f
-                $(F77) $(FFLAGIO) -c treadeo.gfsio.f
-
-twriteeo.o:        twriteeo.f
-                $(F77) $(FFLAGIO) -c twriteeo.f
-
-spect_send.o:        spect_send.f
-                $(F77) $(FFLAGIO) -c spect_send.f
-
-spect_write.o:        spect_write.f
-                $(F77) $(FFLAGIO) -c spect_write.f
-
-spect_to_grid.o:        spect_to_grid.f
-                $(F77) $(FFLAGIO) -c spect_to_grid.f
-
-spect_tv_enthalpy_ps.o:        spect_tv_enthalpy_ps.f
-                $(F77) $(FFLAGIO) -c spect_tv_enthalpy_ps.f
-
-grid_to_spec.o:        grid_to_spec.f
-                $(F77) $(FFLAGIO) -c grid_to_spec.f
-
-wrtout.o:        wrtout.f
-                $(F77) $(FFLAGIO) -c wrtout.f
-
-nstio_module.o:         nstio_module.f
-                $(F77) $(FFLAG90) -c nstio_module.f
-
-module_nst_parameters.o:        module_nst_parameters.f
-                $(F77) $(FFLAG90) -c module_nst_parameters.f
-
-module_nst_water_prop.o:         module_nst_water_prop.f
-                $(F77) $(FFLAG90) -c module_nst_water_prop.f
-
-module_nst_model.o:                module_nst_model.f
-                $(F77) $(FFLAG90) -c module_nst_model.f
-
-NameList_ESMFMod.o:         NameList_ESMFMod.f
-                $(F77) $(FFLAG90) -c NameList_ESMFMod.f
-
-Sfc_Var_ESMFMod.o:        Sfc_Var_ESMFMod.f
-                $(F77) $(FFLAG90) -c Sfc_Var_ESMFMod.f
-
-Nst_Var_ESMFMod.o:        Nst_Var_ESMFMod.f
-                $(F77) $(FFLAG90) -c Nst_Var_ESMFMod.f
-
-Lib_ESMFStateAddGetMod.o:        Lib_ESMFStateAddGetMod.f
-                $(F77) $(FFLAG90) -c Lib_ESMFStateAddGetMod.f
-
-#GFS_ESMFStateAddGetMod.o:         GFS_ESMFStateAddGetMod.f
-#                $(F77) $(FFLAG90) -c GFS_ESMFStateAddGetMod.f
-
-GFS_InternalState_ESMFMod.o:         GFS_InternalState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_InternalState_ESMFMod.f
-
-GFS_ESMFStateMod.o:         GFS_ESMFStateMod.f
-                $(F77) $(FFLAG90) -c GFS_ESMFStateMod.f
-
-GFS_ErrMsgMod.o:         GFS_ErrMsgMod.f
-                $(F77) $(FFLAG90) -c GFS_ErrMsgMod.f
-
-GFS_GetCf_ESMFMod.o:         GFS_GetCf_ESMFMod.f cmp.comm.o atm.comm.o
-                $(F77) $(FFLAG90) -c GFS_GetCf_ESMFMod.f
-
-GFS_ESMFMod.o:         GFS_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_ESMFMod.f
-
-#GFS_Grid_fnl_ESMFMod.o:         GFS_Grid_fnl_ESMFMod.f
-        #$(F77) $(FFLAG90A) -c GFS_Grid_fnl_ESMFMod.f
-
-GFS_GridComp_ESMFMod.o:         GFS_GridComp_ESMFMod.f
-                $(F77) $(FFLAG90A) -c GFS_GridComp_ESMFMod.f
-
-GFS_Initialize_ESMFMod.o:         GFS_Initialize_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Initialize_ESMFMod.f
-
-GFS_Run_ESMFMod.o:        GFS_Run_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Run_ESMFMod.f
-
-GFS_Finalize_ESMFMod.o:        GFS_Finalize_ESMFMod.f
-                $(F77) $(FFLAG90) -c GFS_Finalize_ESMFMod.f
-
-GFS_InputFile2ImportState.o:         GFS_InputFile2ImportState.f
-                $(F77) $(FFLAG90) -c GFS_InputFile2ImportState.f
-
-#GFS_Standalone_ESMF_ENS.o:         GFS_Standalone_ESMF_ENS.f
-#                $(F77) $(FFLAG90) -c GFS_Standalone_ESMF_ENS.f
-
-GFS_ESMF.o:        GFS_ESMF.f
-                $(F77) $(FFLAG90) -c GFS_ESMF.f
-
-Grid_ESMFCreate.o:         Grid_ESMFCreate.f
-                $(F77) $(FFLAG90) -c Grid_ESMFCreate.f
-
-StartTimeGet_ESMF.o:         StartTimeGet_ESMF.f
-                $(F77) $(FFLAG90) -c StartTimeGet_ESMF.f
-
-Ensemble_sub.o:         Ensemble_sub.f
-                $(F77) $(FFLAG90) -c Ensemble_sub.f
-
-mpi_more.o:        mpi_more.f
-                $(F77) $(FFLAGS) -c mpi_more.f
-
-cmp.comm.o:        cmp.comm.f
-                $(F77) $(FFLAGS) -c cmp.comm.f
-
-atm.comm.o:        atm.comm.f
-                $(F77) $(FFLAGS) -c atm.comm.f
-
-tiles.o:        tiles.f
-                $(F77) $(FFLAGS) -c tiles.f
-
-GEFS_Cpl_InternalState_ESMFMod.o:        GEFS_Cpl_InternalState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_InternalState_ESMFMod.f
-
-GEFS_CplState_ESMFMod.o:        GEFS_CplState_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_CplState_ESMFMod.f
-
-GEFS_Sto_Per_Scheme.o:        GEFS_Sto_Per_Scheme.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme.f
-
-GEFS_Cpl_Run_ESMFMod.o:        GEFS_Cpl_Run_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_Run_ESMFMod.f
-
-GEFS_Cpl_ESMFMod.o:        GEFS_Cpl_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_Cpl_ESMFMod.f
-
-GEFS_CplComp_ESMFMod.o:        GEFS_CplComp_ESMFMod.f
-                $(F77) $(FFLAG90) -c GEFS_CplComp_ESMFMod.f
-
-GFS_AddParameterToStateMod.o:        GFS_AddParameterToStateMod.f
-                $(F77) $(FFLAG90) -c GFS_AddParameterToStateMod.f
-
-GEFS_Sto_Per_Scheme_Step1.o:        GEFS_Sto_Per_Scheme_Step1.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme_Step1.f
-
-GEFS_Sto_Per_Scheme_Step2.o:        GEFS_Sto_Per_Scheme_Step2.f
-                $(F77) $(FFLAG90) -c GEFS_Sto_Per_Scheme_Step2.f
-
-GEFS_GetParameterFromStateMod.o:        GEFS_GetParameterFromStateMod.f
-                $(F77) $(FFLAG90) -c GEFS_GetParameterFromStateMod.f
-
-Cal_Sto_Coef.o:                GEFS_Cpl_Cal_Sto_Coef.fd/Cal_Sto_Coef.f90
-                $(F90) $(FFLAG90) -c GEFS_Cpl_Cal_Sto_Coef.fd/Cal_Sto_Coef.f90
-
-GEFS_bcst_global.o:        GEFS_bcst_global.f
-                $(F77) $(FFLAG90) -c GEFS_bcst_global.f

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/gbphys.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/gbphys.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/gbphys.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -420,10 +420,10 @@
 !    &amp; dxmax=ln(1.0/(5000.0*2500.0)),  dxmin=ln(1.0/(192.0*94.0)
 !    &amp; dxmax=-16.118095651,dxmin=-9.800790154,dxinv=1.0/(dxmax-dxmin)
      &amp; dxmaxs=-16.118095651,dxmins=-9.800790154,                        &amp;
-     &amp; dxinvs=1.0/(dxmaxs-dxmins),
+     &amp; dxinvs=1.0/(dxmaxs-dxmins),                                      &amp;
 !    &amp; dxmax=ln(1.0/(2000.0*1000.0)),  dxmin=ln(1.0/(192.0*94.0)
      &amp; dxmaxr=-14.5086577385,dxminr=-9.800790154,                       &amp;
-     &amp; dxinvr=1.0/(dxmaxr-dxminr)
+     &amp; dxinvr=1.0/(dxmaxr-dxminr)                                       
 !    &amp; dxmax=ln(1.0/(2500.0*1250.0)),  dxmin=ln(1.0/(192.0*94.0)
 !    &amp; dxmax=-14.95494484115,dxmin=-9.800790154, dxinv=1.0/(dxmax-dxmin)
        real(kind=kind_phys) dxmax, dxmin, dxinv
@@ -913,7 +913,7 @@
      &amp;          .and. (.not. invrsn(i))) then
               tem = (tgrs(i,k+1)-tgrs(i,k)) / (prsl(i,k)-prsl(i,k+1))
 
-              if ((tem &gt; 0.00010 .and. tx1(i) &lt; 0.0) .or. 
+              if ((tem &gt; 0.00010 .and. tx1(i) &lt; 0.0) .or.               &amp; 
      &amp;            (tem-abs(tx1(i)) &gt; 0.0 .and. tx2(i) &lt; 0.0)) then
                 invrsn(i) = .true.
 
@@ -1072,7 +1072,7 @@
           if ( nst_fcst &gt; 1 ) then
             do i = 1, im
               if ( slmsk(i) == 0.0 ) then
-                tsea(i) = max(271.0, tref(i) + (xt(i)+xt(i))/xz(i)
+                tsea(i) = max(271.0, tref(i) + (xt(i)+xt(i))/xz(i)      &amp;
      &amp;                                       - dt_cool(i))              &amp;
      &amp;                                       - oro(i) * rlapse
               endif
@@ -2219,7 +2219,7 @@
 !           print *,' gq0b=',gq0(ipr,:,1)
 !         endif
 
-          call mstcnv(im,ix,levs,dtp,gt0,gq0,prsl,del,prslk,rain1
+          call mstcnv(im,ix,levs,dtp,gt0,gq0,prsl,del,prslk,rain1 &amp;
      &amp;,                            gq0(1,1,ntcw), rhc, lprnt,ipr)
 
 !         if (lprnt) then
@@ -2237,7 +2237,7 @@
             if (ldiag3d) then
               do k=1,levs
                 do i=1,im
-                  dt3dt(i,k,4) = dt3dt(i,k,4) + (gt0(i,k)-dtdt(i,k))
+                  dt3dt(i,k,4) = dt3dt(i,k,4) + (gt0(i,k)-dtdt(i,k)) &amp;
      &amp;                                        * frain
 !                 dq3dt(i,k,2) = dq3dt(i,k,2) + (gq0(i,k,1)-dqdt(i,k,1))
 !    &amp;                                        * frain
@@ -2247,7 +2247,7 @@
             if (ldiag3d .or. lggfs3d) then
               do k=1,levs
                 do i=1,im
-                  dq3dt(i,k,2) = dq3dt(i,k,2) + (gq0(i,k,1)-dqdt(i,k,1))
+                  dq3dt(i,k,2) = dq3dt(i,k,2) + (gq0(i,k,1)-dqdt(i,k,1)) &amp;
      &amp;                                        * frain
                 enddo
               enddo
@@ -2351,9 +2351,9 @@
       if (cal_pre) then
 !       HCHUANG: add dominant precipitation type algorithm
 
-        call calpreciptype(kdt,nrcm,im,ix,LEVS,LEVS+1,rann,
-     &amp;    xlat,xlon,gt0,gq0,prsl,prsi,RAIN,
-     &amp;    phii,num_p3d,tsea,sr,phy_f3d(1,1,3),    ! input
+        call calpreciptype(kdt,nrcm,im,ix,LEVS,LEVS+1,rann, &amp;
+     &amp;    xlat,xlon,gt0,gq0,prsl,prsi,RAIN,                 &amp;
+     &amp;    phii,num_p3d,tsea,sr,phy_f3d(1,1,3), &amp;  ! input
      &amp;    DOMR,DOMZR,DOMIP,DOMS)                  ! output
 
 !

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/machine.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/machine.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/machine.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1,4 +1,5 @@
       MODULE MACHINE
+      use mpas_kind_types
 
       IMPLICIT NONE
       SAVE
@@ -6,10 +7,12 @@
       integer, parameter :: kind_io4  = 4, kind_io8  = 8 , kind_ior = 8
      &amp;,                     kind_evod = 8, kind_dbl_prec = 8
      &amp;,                     kind_rad  = selected_real_kind(13,60) ! the '60' maps to 64-bit real
-     &amp;,                     kind_phys = selected_real_kind(13,60) ! the '60' maps to 64-bit real
+!    &amp;,                     kind_phys = selected_real_kind(13,60) ! the '60' maps to 64-bit real
      &amp;,                     kind_REAL = 8                         ! used in cmp_comm
      &amp;,                     kind_INTEGER = 4                      ! -,,-
 !
+      integer, parameter:: kind_phys = RKIND
+
       real(kind=kind_evod), parameter :: mprec = 1.e-12           ! machine precision to restrict dep
 
       END MODULE MACHINE

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/module_bfmicrophysics.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/module_bfmicrophysics.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/module_bfmicrophysics.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -228,13 +228,13 @@
 !       read(1) my_growth    ! Applicable only for DTPH=180 s for offline testing
         CLOSE (1)
       else
-        etime1=timef()
+!       etime1=timef()
         CALL ICE_LOOKUP                   ! Lookup tables for ice
-        etime2=timef()
+!       etime2=timef()
         if (mype == 0)                                                  &amp;
      &amp;  print *,'CPU time (sec) in ICE_LOOKUP = ',(etime2-etime1)*0.001
         CALL RAIN_LOOKUP                  ! Lookup tables for rain
-        etime1=timef()
+!       etime1=timef()
         if (mype == 0)                                                  &amp;
      &amp;  print *,'CPU time (sec) in RAIN_LOOKUP = ',(etime1-etime2)*0.001
         if (write_lookup) then

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/para_nstio_w.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/para_nstio_w.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/para_nstio_w.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -52,9 +52,9 @@
         head%idate   = idate
 !
         call nstio_aldbta(head,data,iret)
-        PRINT 99,nw,xhour,IDATE,iret
-99      FORMAT(1H ,'in para_nstio_w nw=',i7,2x,'HOUR=',f8.2,3x,'IDATE=',
-     &amp;  4(1X,I4)),' iret=',I2
+!       PRINT 99,nw,xhour,IDATE,iret
+!99     FORMAT(1H ,'in para_nstio_w nw=',i7,2x,'HOUR=',f8.2,3x,'IDATE=',
+!    &amp;  4(1X,I4)),' iret=',I2
       ENDIF
 !!
       kmsk= nint(nst_fld%slmsk)

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radiation_clouds.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/radiation_clouds.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radiation_clouds.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -216,24 +216,34 @@
 !===&gt; ...  begin here
 !
       if (icld == 0) then
+!ldf begin (2012-11-09):
         if (me == 0) print *,' - Using Diagnostic Cloud Method'
+        if (me == 0) write(0,*)'--- Using Diagnostic Cloud Method'
+!ldf end.
 
 !  ---  set up tuned rh table
 
         call rhtable( me, ier )
 
         if (ier &lt; 0) then
-          write(6,99) ier
+!ldf begin (2012-11-09):
+!         write(6,99) ier
+          write(0,99) ier
+!ldf end.
   99      format(3x,' *** Error in finding tuned RH table ***'          &amp;
      &amp;,         /3x,'     STOP at calling subroutine RHTABLE !!'/)
           stop 99
         endif
       else
         if (me == 0) then
-          print *,' - Using Prognostic Cloud Method'
-          if (np3d == 3) print *,'   --- Ferrier cloud microphysics'
-          if (np3d == 4) print *,                                       &amp;
-     &amp;                      '   --- Zhao/Carr/Sundqvist microphysics'
+!ldf begin (2012-11-09):
+!         print *,' - Using Prognostic Cloud Method'
+!         if (np3d == 3) print *,'   --- Ferrier cloud microphysics'
+!         if (np3d == 4) print *,                                       &amp;
+!    &amp;                      '   --- Zhao/Carr/Sundqvist microphysics'
+          if (np3d == 3) write(0,*) '--- Using Ferrier cloud microphysics'
+          if (np3d == 4) write(0,*) '--- Using Zhao/Carr/Sundqvist microphysics'
+!ldf end.
         endif
       endif
 

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radiation_gases.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/radiation_gases.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radiation_gases.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -220,8 +220,11 @@
         co2_glb = co2vmr_def
 
         if ( me == 0 ) then
-          print *,' - Using prescribed co2 global mean value=',         &amp;
-     &amp;              co2vmr_def
+!ldf (2012-11-09) begin:
+!        print *,' - Using prescribed co2 global mean value=',         &amp;
+!    &amp;              co2vmr_def
+         write(0,*) '---- Using prescribed co2 global mean value=', co2vmr_def
+!ldf end.
         endif
 
         return

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radlw_main.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/radlw_main.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radlw_main.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1174,41 +1174,75 @@
       iovrlw  = iovr     ! assign module variable of overlap flag
 
       if ( iovrlw&lt;0 .or. iovrlw&gt;1 ) then
-        print *,'  *** Error in specification of cloud overlap flag',   &amp;
-     &amp;          ' IOVRLW=',iovrlw,' in RLWINIT !!'
+!ldf begin (2012-11-09):
+!       print *,'  *** Error in specification of cloud overlap flag',   &amp;
+!    &amp;          ' IOVRLW=',iovrlw,' in RLWINIT !!'
+        write(0,*) '  *** Error in specification of cloud overlap flag',   &amp;
+     &amp;          '         IOVRLW=',iovrlw,' in RLWINIT !!'
+!ldf end.
         stop
       endif
 
       if (me == 0) then
-        print *,' - Using AER Longwave Radiation, Version: ', VTAGLW
-
+!ldf begin (2012-11-09):
+!       print *,' - Using AER Longwave Radiation, Version: ', VTAGLW
+        write(0,*)'--- Using AER Longwave Radiation, Version: ', VTAGLW
+!ldf end.
         if (iaerlw &gt; 0) then
-          print *,'   --- Using input aerosol parameters for LW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Using input aerosol parameters for LW'
+          write(0,*) '--- Using input aerosol parameters for LW'
+!ldf end.
         else
-          print *,'   --- Aerosol effect is NOT included in LW, all'    &amp;
-     &amp;           ,' internal aerosol parameters are reset to zeros'
+!ldf begin (2012-11-09):
+!         print *,'   --- Aerosol effect is NOT included in LW, all'    &amp;
+!    &amp;           ,' internal aerosol parameters are reset to zeros'
+          write(0,*) '--- Aerosol effect is NOT included in LW, all'    &amp;
+     &amp;               ,' internal aerosol parameters are reset to zeros'
+!ldf end.
         endif
 
         if (irgaslw == 1) then
-          print *,'   --- Include rare gases N2O, CH4, O2, absorptions',&amp;
-     &amp;            ' in LW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Include rare gases N2O, CH4, O2, absorptions',&amp;
+!    &amp;            ' in LW'
+          write(0,*) '--- Include rare gases N2O, CH4, O2, absorptions in LW'
+!ldf end.
         else
-          print *,'   --- Rare gases effect is NOT included in LW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Rare gases effect is NOT included in LW'
+          write(0,*) '--- Rare gases effect is NOT included in LW'
+!ldf end.
         endif
 
         if (icfclw == 1) then
+!ldf begin (2012-11-09):
           print *,'   --- Include CFC gases absorptions in LW'
+          write(0,*)'--- Include CFC gases absorptions in LW'
+!ldf end.
         else
-          print *,'   --- CFC gases effect is NOT included in LW'
+!ldf begin (2012-11-09):
+!         print *,'   --- CFC gases effect is NOT included in LW'
+          write(0,*) '--- CFC gases effect is NOT included in LW'
+!ldf end.
         endif
 
         if ( isubc == 0 ) then
-          print *,'   --- Using standard grid average clouds, no sub-', &amp;
+!ldf begin (2012-11-09):
+!         print *,'   --- Using standard grid average clouds, no sub-', &amp;
+!    &amp;            'column clouds approximation'
+          write(0,*) '--- Using standard grid average clouds, no sub-', &amp;
      &amp;            'column clouds approximation'
+!ldf end.
         else
-          print *,'   --- Sub-column cloud scheme is not available in', &amp;
+!ldf begin (2012-11-09):
+!         print *,'   --- Sub-column cloud scheme is not available in', &amp;
+!    &amp;            ' this version of code.  Using standard grid',        &amp;
+!    &amp;            ' average of clouds'
+          write(0,*) '--- Sub-column cloud scheme is not available in', &amp;
      &amp;            ' this version of code.  Using standard grid',        &amp;
      &amp;            ' average of clouds'
+!ldf end.
         endif
       endif
 
@@ -1216,8 +1250,12 @@
 
       if ((icwp == 0 .and. iflagliq /= 0) .or.                          &amp;
      &amp;    (icwp == 1 .and. iflagliq == 0)) then
+!ldf begin (2012-11-09): 
         print *, ' *** Model cloud scheme inconsistent with LW',        &amp;
      &amp;           ' radiation cloud radiative property setup !!'
+        write(0,*) ' *** Model cloud scheme inconsistent with LW',        &amp;
+     &amp;             ' radiation cloud radiative property setup !!'
+!ldf end.
         stop
       endif
 

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radsw_main.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/radsw_main.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/radsw_main.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1179,43 +1179,79 @@
       iovrsw  = iovr     ! assign module variable of overlap flag
 
       if ( iovrsw&lt;0 .or. iovrsw&gt;2 ) then
-        print *,'  *** Error in specification of cloud overlap flag',   &amp;
-     &amp;          ' IOVRSW=',iovrsw,' in RSWINIT !!'
+!ldf begin (2012-11-09):
+!       print *,'  *** Error in specification of cloud overlap flag',   &amp;
+!    &amp;          ' IOVRSW=',iovrsw,' in RSWINIT !!'
+        write(0,*) '  *** Error in specification of cloud overlap flag',   &amp;
+     &amp;             ' IOVRSW=',iovrsw,' in RSWINIT !!'
+!ldf end.
         stop
       endif
 
       if (me == 0) then
-        print *,' - Using AER Shortwave Radiation, Version: ',VTAGSW
+!ldf begin (2012-11-09):
+!       print *,' - Using AER Shortwave Radiation, Version: ',VTAGSW
+        write(0,*) '--- Using AER Shortwave Radiation, Version: ',VTAGSW
+!ldf end.
 
         if (imodsw == 1) then
-          print *,'   --- Delta-eddington 2-stream transfer scheme'
+!ldf begin (2012-11-09):
+!         print *,'   --- Delta-eddington 2-stream transfer scheme'
+          write(0,*) '--- Delta-eddington 2-stream transfer scheme'
+!end
         else if (imodsw == 2) then
-          print *,'   --- PIFM 2-stream transfer scheme'
+!ldf begin (2012-11-09):
+!         print *,'   --- PIFM 2-stream transfer scheme'
+          write(0,*) '--- PIFM 2-stream transfer scheme'
+!ldf end.
         else if (imodsw == 3) then
-          print *,'   --- Discrete ordinates 2-stream transfer scheme'
+!ldf begin (2012-11-09):
+!         print *,'   --- Discrete ordinates 2-stream transfer scheme'
+          write(0,*) '--- Discrete ordinates 2-stream transfer scheme'
         endif
 
         if (iaersw == 0) then
-          print *,'   --- Aerosol effect is NOT included in SW, all'    &amp;
-     &amp;           ,' internal aerosol parameters are reset to zeros'
+!ldf begin (2012-11-09):
+!         print *,'   --- Aerosol effect is NOT included in SW, all'    &amp;
+!    &amp;           ,' internal aerosol parameters are reset to zeros'
+          write(0,*) '--- Aerosol effect is NOT included in SW, all'    &amp;
+     &amp;              ,' internal aerosol parameters are reset to zeros'
         else
-          print *,'   --- Using input aerosol parameters for SW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Using input aerosol parameters for SW'
+          write(0,*) '--- Using input aerosol parameters for SW'
+!ldf end.
         endif
 
         if (irgassw == 0) then
-          print *,'   --- Rare gases absorption is NOT included in SW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Rare gases absorption is NOT included in SW'
+          write(0,*) '--- Rare gases absorption is NOT included in SW'
+!ldf end.
         else
-          print *,'   --- Include rare gases N2O, CH4, O2, absorptions',&amp;
-     &amp;            ' in SW'
+!ldf begin (2012-11-09):
+!         print *,'   --- Include rare gases N2O, CH4, O2, absorptions',&amp;
+!    &amp;            ' in SW'
+          write(0,*) '--- Include rare gases N2O, CH4, O2, absorptions in SW'
+!ldf end.
         endif
 
         if ( isubc == 0 ) then
-          print *,'   --- Using standard grid average clouds, no sub-', &amp;
-     &amp;            'column clouds approximation'
+!ldf begin (2012-11-09):
+!         print *,'   --- Using standard grid average clouds, no sub-', &amp;
+!    &amp;            'column clouds approximation'
+          write(0,*) '--- Using standard grid average clouds, no sub-', &amp;
+     &amp;               'column clouds approximation'
+!ldf end.
         else
-          print *,'   --- Sub-column cloud scheme is not available in', &amp;
-     &amp;            ' this version of code.  Using standard grid',        &amp;
-     &amp;            ' average of clouds'
+!ldf begin (2012-11-09):
+!         print *,'   --- Sub-column cloud scheme is not available in', &amp;
+!    &amp;            ' this version of code.  Using standard grid',        &amp;
+!    &amp;            ' average of clouds'
+          write(0,*) '--- Sub-column cloud scheme is not available in', &amp;
+     &amp;               ' this version of code.  Using standard grid',        &amp;
+     &amp;               ' average of clouds'
+!ldf end.
         endif
       endif
 
@@ -1223,8 +1259,12 @@
 
       if ((icwp == 0 .and. iflagliq /= 0) .or.                          &amp;
      &amp;    (icwp == 1 .and. iflagliq == 0)) then
-        print *, ' *** Model cloud scheme inconsistent with SW',        &amp;
+!ldf begin (2012-11-09):
+!       print *, ' *** Model cloud scheme inconsistent with SW',        &amp;
+!    &amp;           ' radiation cloud radiative property setup !!'
+        write(0,*)' *** Model cloud scheme inconsistent with SW',        &amp;
      &amp;           ' radiation cloud radiative property setup !!'
+!ldf end.
         stop
       endif
 
@@ -1971,7 +2011,8 @@
       implicit none
 
       real (kind=kind_phys), parameter :: zcrit = 0.9995   ! thresold for conservative scattering
-      real (kind=kind_phys), parameter :: zsr3  = sqrt(3.0)
+!     real (kind=kind_phys), parameter :: zsr3  = sqrt(3.0)
+      real (kind=kind_phys), parameter :: zsr3  = 1.732050807568877
 
 !  ---  inputs:
       integer, intent(in) :: ib

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/rascnvv2.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/rascnvv2.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/rascnvv2.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -275,7 +275,7 @@
      &amp;,                    RAINC(im),     CDRAG(im),  DDVEL(im)         &amp;
      &amp;,                    rannum(ix,nrcm),dlqfac                       &amp;
      &amp;,                    ud_mf(im,k), dd_mf(im,k), det_mf(im,k)
-      real(kind=kind_phys) DT, facmb, garea(im), dtf, rhc(im,k)         &amp;
+      real(kind=kind_phys) DT, facmb, garea(im), dtf, rhc(im,k)
 !
 !     Added for aerosol scavenging for GOCART
 !
@@ -297,7 +297,7 @@
       PARAMETER (ICM=100, DAYLEN=86400.0, PFAC=1.0/450.0,clwmin=1.0e-10)
       Integer  IC(ICM)
 !
-      real(kind=kind_phys), allocatable ::  ALFINT(:,:), uvi(:,:)
+      real(kind=kind_phys), allocatable ::  ALFINT(:,:), uvi(:,:)       &amp;
      &amp;,                                     trcfac(:,:), rcu(:,:)
       real(kind=kind_phys)            ALFINQ(K),   PRSM(K),  PSJM(K)    &amp;
      &amp;,                    alfind(K), rhc_l(k), dtvd(2,4)
@@ -1047,7 +1047,7 @@
      &amp;                     dhdpmn, dhdp(kd:k)
       parameter (ERRMIN=0.0001, ERRMI2=0.1*ERRMIN)
       INTEGER I, L,  N,  KD1, II                                        &amp;
-     &amp;,       KP1, IT, KM1, KTEM, KK, KK1, LM1, LL, LP1, kbls, kmxh
+     &amp;,       KP1, IT, KM1, KTEM, KK, KK1, LM1, LL, LP1, kbls, kmxh     &amp;
      &amp;,       kblh, kblm, kblpmn, kmax, kmaxm1, kmaxp1, klcl, kmin, kmxb
 
       real avt, avq, avr, avh

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sfc_land.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/sfc_land.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sfc_land.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -12,7 +12,7 @@
      &amp;       qsurf, snowmt, gflux, zsoil, rhscnpy, rhsmc,               &amp;
      &amp;       aim, bim, cim, drain, evap, hflx, ep, cmm, chh,            &amp;
      &amp;       evbs, evcw, trans, sbsno, snowc, stm, snohf,               &amp;
-     &amp;       twilt, tref  
+     &amp;       twilt, tref                                                &amp; 
      &amp;     )
 
 ! ===================================================================== !

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sfc_nst.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/sfc_nst.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sfc_nst.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -1,16 +1,16 @@
 
-      SUBROUTINE SFC_NST                                                &amp;
+      SUBROUTINE SFC_NST
 !...................................
 !  ---  inputs:
-     &amp;     ( IM, KM, PS, U1, V1, T1, Q1, tref, cm, ch,                  &amp;
-     &amp;       prsl1, prslki, slimsk, xlon, sinlat, stress,               &amp;
-     &amp;       sfcemis, dlwflx, sfcnsw, rain, timestep, kdt,              &amp;
-     &amp;       ddvel, flag_iter, flag_guess, nst_fcst, lprnt, ipr,        &amp;
+     &amp;     ( IM, KM, PS, U1, V1, T1, Q1, tref, cm, ch,
+     &amp;       prsl1, prslki, slimsk, xlon, sinlat, stress,
+     &amp;       sfcemis, dlwflx, sfcnsw, rain, timestep, kdt,
+     &amp;       ddvel, flag_iter, flag_guess, nst_fcst, lprnt, ipr,
 !  --- Input/output
-     &amp;       tskin, tsurf, xt, xs, xu, xv, xz, zm, xtts, xzts, dt_cool, &amp;
-     &amp;       z_c,   c_0,   c_d,   w_0, w_d, d_conv, ifd, Qrain,         &amp;
+     &amp;       tskin, tsurf, xt, xs, xu, xv, xz, zm, xtts, xzts, dt_cool,
+     &amp;       z_c,   c_0,   c_d,   w_0, w_d, d_conv, ifd, Qrain,
 !  ---  outputs:
-     &amp;       qsurf, gflux, CMM, CHH, EVAP, HFLX, EP                     &amp;
+     &amp;       qsurf, gflux, CMM, CHH, EVAP, HFLX, EP
      &amp;      )
 !
 ! ===================================================================== !
@@ -143,8 +143,8 @@
 
 !  ---  inputs:
       integer, intent(in) :: im, km, kdt, ipr, nst_fcst
-      real (kind=kind_phys), dimension(im), intent(in) :: ps, u1, v1,   &amp;
-     &amp;       t1, q1, tref, cm, ch, prsl1, prslki, slimsk, xlon,         &amp;
+      real (kind=kind_phys), dimension(im), intent(in) :: ps, u1, v1,
+     &amp;       t1, q1, tref, cm, ch, prsl1, prslki, slimsk, xlon,
      &amp;       sinlat, stress, sfcemis, dlwflx, sfcnsw, rain, ddvel
 
       real (kind=kind_phys), intent(in) :: timestep
@@ -153,12 +153,12 @@
 
 !  ---  input/outputs:
 ! control variables of DTL system (5+2) and SL (2) and coefficients for d(Tz)/d(Ts) calculation
-      real (kind=kind_phys), dimension(im), intent(inout) :: tskin,     &amp;
-     &amp;      tsurf, xt, xs, xu, xv, xz, zm, xtts, xzts, dt_cool,         &amp;
+      real (kind=kind_phys), dimension(im), intent(inout) :: tskin,
+     &amp;      tsurf, xt, xs, xu, xv, xz, zm, xtts, xzts, dt_cool,
      &amp;      z_c, c_0, c_d, w_0, w_d, d_conv, ifd, qrain
 
 !  ---  outputs:
-      real (kind=kind_phys), dimension(im), intent(out) ::              &amp;
+      real (kind=kind_phys), dimension(im), intent(out) ::
      &amp;       qsurf, gflux, cmm, chh, evap, hflx, ep
 
 !
@@ -198,7 +198,6 @@
 !  external functions called: iw3jdn
       integer :: iw3jdn
 !======================================================================================================
-cc
       PARAMETER (elocp=hvap/cp)
 
       sss = 34.0             ! temporarily, when sea surface salinity data is not ready

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sflx.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/sflx.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/sflx.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -4609,7 +4609,7 @@
 
         parameter (cvfrz = 3)
 
-c ----------------------------------------------------------------------
+! ----------------------------------------------------------------------
 !  --- ...  determine rainfall infiltration rate and runoff.  include
 !           the infiltration formule from schaake and koren model.
 !           modified by q duan

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/tracer_const_h-new.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/tracer_const_h-new.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/tracer_const_h-new.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -15,7 +15,7 @@
       integer ntrac,me,nlunit
       namelist /tracer_constant/ ri,cpi
 
-c
+!
       if( ntrac.ne.num_tracer ) then
         if( me.eq.0 ) then
           write(*,*) ' Error ; inconsistent number of tracer '

Modified: branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/wrtnst.f
===================================================================
--- branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/core_atmos_physics_gfs/physics_gfs/wrtnst.f        2012-10-30 17:55:34 UTC (rev 2283)
+++ branches/atmphys_gfs_mmm/src/core_atmos_physics_gfs/physics_gfs/wrtnst.f        2013-01-15 23:00:59 UTC (rev 2386)
@@ -112,7 +112,8 @@
      $               IHCBL , IHCTL , IMCLYR, IMCTL , IMCBL , IMCTL ,
      $               ILCLYR, ILCTL , ILCBL , ILCTL , ISFC  , ISFC /
 !    $               ILCLYR, ILCTL , ILCBL , ILCTL /
-      INTEGER     IDATE(4), IDS(255),IENS(5)
+!     INTEGER     IDATE(4), IDS(255),IENS(5)
+      INTEGER     IDATE(4), IDS(255)
       real (kind=kind_io8) SI(LEVP1)
 C
 csela..................................................................

</font>
</pre>