Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/WRFV2-gfortran-x86_64-2.patch =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/WRFV2-gfortran-x86_64-2.patch (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/WRFV2-gfortran-x86_64-2.patch (revision 1185) @@ -0,0 +1,361 @@ +--- WRFV2.original/arch/configure.defaults 2007-10-05 16:16:41.329127000 +0200 ++++ WRFV2/arch/configure.defaults 2008-09-03 20:18:02.656065000 +0200 +@@ -10943,5 +10943,358 @@ + $(RM) $*.b + $(FC) -c $(PROMOTION) $(FCBASEOPTS) $(MODULE_DIRS) $*.f90 + ++########################################################### ++#ARCH PC Linux x86_64 (IA64 and AMD Opteron), GNU Fortran compiler, DM-Parallel (RSL, MPICH, Allows nesting) ++# ++DMPARALLEL = 1 ++MAX_PROC = 1024 ++FC = mpif90 ++LD = mpif90 ++CC = mpicc ++SCC = gcc ++SFC = gfortran ++RWORDSIZE = $(NATIVE_RWORDSIZE) ++PROMOTION = ++CFLAGS = -DDM_PARALLEL -DWRF_RSL_IO \ ++ -DMAXDOM_MAKE=$(MAX_DOMAINS) -DMAXPROC_MAKE=$(MAX_PROC) -I../external/RSL/RSL \ ++ -DLANDREAD_STUB ++FCOPTIM = #-Mipa=fast ++FCDEBUG = #-g ++FCBASEOPTS = -w $(FCDEBUG) -ffree-form -ffree-line-length-none ++FCFLAGS = $(FCOPTIM) $(FCBASEOPTS) ++ARCHFLAGS = -DGRIB1 CONFIGURE_GRIB2_FLAG -DINTIO -DWRF_RSL_IO -DRSL -DDM_PARALLEL \ ++ -DIWORDSIZE=4 -DDWORDSIZE=8 -DRWORDSIZE=$(RWORDSIZE) -DLWORDSIZE=4 CONFIGURE_NETCDF_FLAG \ ++ -DLIMIT_ARGS -DLANDREAD_STUB ++# machine-specific flags needed to link in ESMF library (C++ run-time-library, etc.) ++ESMF_LIB_FLAGS = ++ESMF_IO_LIB = ESMFIOLIB ++ESMF_IO_LIB_EXT = ESMFIOEXTLIB ++INCLUDE_MODULES = -J../main -I../main -I../external/io_netcdf -I../external/io_int -I../external/esmf_time_f90 \ ++ -I../external -I../frame -I../share -I../phys -I../inc -I../chem ++PERL = perl ++REGISTRY = Registry ++LIB = CONFIGURE_NETCDF_LIB_PATH -L../external/RSL/RSL -lrsl \ ++ -L../external/io_grib1 -lio_grib1 \ ++ CONFIGURE_GRIB2_LIB \ ++ -L../external/io_grib_share -lio_grib_share \ ++ -L../external/io_int -lwrfio_int \ ++ ../frame/module_internal_header_util.o ../frame/pack_utils.o -L../external/esmf_time_f90 -lesmf_time ++LDFLAGS = CONFIGURE_LDFLAGS ++ENVCOMPDEFS = CONFIGURE_COMPILEFLAGS ++CPP = /lib/cpp -C -P -traditional ++POUND_DEF = -DNO_RRTM_PHYSICS -traditional $(COREDEFS) -DNONSTANDARD_SYSTEM -DCONFIG_BUF_LEN=$(CONFIG_BUF_LEN) -DMAX_DOMAINS_F=$(MAX_DOMAINS) ++CPPFLAGS = -I$(LIBINCLUDE) -C -P $(ARCHFLAGS) -I../external/RSL/RSL -C -P `cat ../inc/dm_comm_cpp_flags` $(ENVCOMPDEFS) $(POUND_DEF) ++AR = ar ru ++M4 = m4 -B 14000 ++RANLIB = ranlib ++NETCDFPATH = CONFIGURE_NETCDF_PATH ++CC_TOOLS = cc ++ ++externals : wrf_ioapi_includes CONFIGURE_WRFIO_NF wrfio_grib_share wrfio_grib1 CONFIGURE_WRFIO_GRIB2 ../external/RSL/RSL/librsl.a wrfio_int gen_comms.c module_dm.F esmf_time ++ ++gen_comms.c : ../external/RSL/gen_comms.c ++ ( /bin/cp ../tools/gen_comms_warning ../tools/gen_comms.c ; \ ++ cat ../external/RSL/gen_comms.c >> ../tools/gen_comms.c ) ++ ++module_dm.F : ../external/RSL/module_dm.F ++ ( /bin/cp module_dm_warning module_dm.F ; cat ../external/RSL/module_dm.F >> module_dm.F ) ++ ++wrfio_nf : ++ ( cd ../external/io_netcdf ; \ ++ make NETCDFPATH=CONFIGURE_NETCDF_PATH RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(FC) $(PROMOTION) $(FCFLAGS) -w" TRADFLAG="-traditional" ) ++ ++wrfio_int : ++ ( cd ../external/io_int ; \ ++ make CC="$(CC)" RANLIB="$(RANLIB)" CPP="$(CPP)" SFC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" FC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" all ) ++ ++wrfio_grib_share : ++ ( cd ../external/io_grib_share ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib1 : ++ ( cd ../external/io_grib1 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib2 : ++ ( cd ../external/io_grib2 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS) CONFIGURE_GRIB2_INC" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" FIXED="-Mfixed" archive) ++ ++esmf_time : ++ ( cd ../external/esmf_time_f90 ; \ ++ make FC="$(FC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" RANLIB="$(RANLIB)" CPP="$(CPP) -I../../inc -I. $(POUND_DEF)" ) ++ ++../external/RSL/RSL/librsl.a : ++ ( cd ../external/RSL/RSL ; make CC="$(CC)" FC="$(FC) $(PROMOTION)" MAX_DOMAINS=$(MAX_DOMAINS) MAX_PROC=$(MAX_PROC) LEARN_BCAST=-DLEARN_BCAST linux ) ++ ++# compile these without high optimization to speed compile ++solve_interface.o : solve_interface.F ++shift_domain_em.o : shift_domain_em.F ++module_io_mm5.o : module_io_mm5.F ++module_si_io.o : module_si_io.F ++module_io_wrf.o : module_io_wrf.F ++module_domain.o : module_domain.F ++module_start.o : module_start.F ++module_initialize.o : module_initialize.F ++module_initialize_b_wave.o : module_initialize_b_wave.F ++module_initialize_hill2d_x.o : module_initialize_hill2d_x.F ++module_initialize_quarter_ss.o : module_initialize_quarter_ss.F ++module_initialize_squall2d_x.o : module_initialize_squall2d_x.F ++module_initialize_squall2d_y.o : module_initialize_squall2d_y.F ++module_initialize_real.o : module_initialize_real.F ++module_dm.o : module_dm.F ++start_domain.o : start_domain.F ++mediation_integrate.o : mediation_integrate.F ++module_configure.o : module_configure.F ++ ++solve_interface.o \ ++shift_domain_em.o \ ++module_io_mm5.o module_si_io.o module_io_wrf.o module_domain.o \ ++module_start.o module_initialize.o module_initialize_b_wave.o \ ++module_initialize_hill2d_x.o module_initialize_quarter_ss.o \ ++module_initialize_squall2d_x.o module_initialize_squall2d_y.o \ ++mediation_integrate.o \ ++module_configure.o \ ++module_initialize_real.o module_dm.o start_domain.o : ++ $(RM) $@ ++ $(SED_FTN) $*.F > $*.b ++ $(CPP) -I../inc $(CPPFLAGS) $*.b > $*.f90 ++ $(RM) $*.b ++ $(FC) -c $(PROMOTION) $(FCBASEOPTS) $(MODULE_DIRS) $*.f90 ++ ++ ++########################################################### ++#ARCH PC Linux x86_64 (IA64 and Opteron), GNU Fortran compiler (Single-threaded, RSL, Allows nesting) ++# ++DMPARALLEL = 1 ++MAX_PROC = 1024 ++FC = gfortran ++LD = gfortran ++CC = gcc ++SCC = $(CC) ++SFC = gfortran ++RWORDSIZE = $(NATIVE_RWORDSIZE) ++FCOPTIM = #-O3 -funroll-loops -ftree-vectorize -ftree-loop-linear ++FCDEBUG = #-g ++FCBASEOPTS = -w $(FCDEBUG) -ffree-form -ffree-line-length-none -fconvert=big-endian -frecord-marker=4 ++FCFLAGS = $(FCOPTIM) $(FCBASEOPTS) ++ARCHFLAGS = -DGRIB1 -DINTIO -DWRF_RSL_IO -DRSL -DDM_PARALLEL \ ++ -DIWORDSIZE=4 -DDWORDSIZE=8 -DRWORDSIZE=$(RWORDSIZE) -DLWORDSIZE=4 CONFIGURE_NETCDF_FLAG \ ++ -DLIMIT_ARGS -DSTUBMPI ++PROMOTION = ++CFLAGS = -DLANDREADSTUB -DDM_PARALLEL -DWRF_RSL_IO \ ++ -DMAXDOM_MAKE=$(MAX_DOMAINS) -DMAXPROC_MAKE=$(MAX_PROC) -I../external/RSL/RSL $(ARCHFLAGS) ++# machine-specific flags needed to link in ESMF library (C++ run-time-library, etc.) ++ESMF_LIB_FLAGS = ++ESMF_IO_LIB = ESMFIOLIB ++ESMF_IO_LIB_EXT = ESMFIOEXTLIB ++INCLUDE_MODULES = -J../main -I../main -I../external/io_netcdf -I../external/io_int -I../external/esmf_time_f90 \ ++ -I../external -I../frame -I../share -I../phys -I../chem -I../inc ++PERL = perl ++REGISTRY = Registry ++LIB = CONFIGURE_NETCDF_LIB_PATH -L../external/RSL/RSL -lrsl \ ++ -L../external/io_grib1 -lio_grib1 \ ++ -L../external/io_int -lwrfio_int \ ++ -L../external/io_grib_share -lio_grib_share \ ++ ../frame/module_internal_header_util.o ../frame/pack_utils.o -L../external/esmf_time_f90 -lesmf_time ++LDFLAGS = $(FCFLAGS) CONFIGURE_LDFLAGS -fconvert=big-endian ++ENVCOMPDEFS = CONFIGURE_COMPILEFLAGS ++CPP = /lib/cpp -C -P -traditional ++POUND_DEF = -DNO_RRTM_PHYSICS $(COREDEFS) -DNONSTANDARD_SYSTEM -DCONFIG_BUF_LEN=$(CONFIG_BUF_LEN) -DMAX_DOMAINS_F=$(MAX_DOMAINS) ++CPPFLAGS = -I$(LIBINCLUDE) -C -P $(ARCHFLAGS) -I../external/RSL/RSL -C -P `cat ../inc/dm_comm_cpp_flags` $(ENVCOMPDEFS) $(POUND_DEF) ++AR = ar ru ++M4 = m4 -B 14000 ++RANLIB = ranlib ++NETCDFPATH = CONFIGURE_NETCDF_PATH ++CC_TOOLS = cc ++ ++externals : wrf_ioapi_includes CONFIGURE_WRFIO_NF wrfio_grib_share wrfio_grib1 ../external/RSL/RSL/librsl.a wrfio_int gen_comms.c module_dm.F esmf_time ++ ++gen_comms.c : ../external/RSL/gen_comms.c ++ ( /bin/cp ../tools/gen_comms_warning ../tools/gen_comms.c ; \ ++ cat ../external/RSL/gen_comms.c >> ../tools/gen_comms.c ) ++ ++module_dm.F : ../external/RSL/module_dm.F ++ ( /bin/cp module_dm_warning module_dm.F ; cat ../external/RSL/module_dm.F >> module_dm.F ) ++ ++wrfio_nf : ++ ( cd ../external/io_netcdf ; \ ++ make NETCDFPATH=CONFIGURE_NETCDF_PATH RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(FC) $(PROMOTION) $(FCFLAGS) -w" TRADFLAG="-traditional" ) ++ ++wrfio_grib_share : ++ ( cd ../external/io_grib_share ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib1 : ++ ( cd ../external/io_grib1 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib2 : ++ ( cd ../external/io_grib2 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS) CONFIGURE_GRIB2_INC" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" FIXED="-Mfixed" archive) ++ ++wrfio_int : ++ ( cd ../external/io_int ; \ ++ make CC="$(CC)" RANLIB="$(RANLIB)" CPP="$(CPP)" SFC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" FC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" all ) ++ ++esmf_time : ++ ( cd ../external/esmf_time_f90 ; \ ++ make FC="$(FC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" RANLIB="$(RANLIB)" CPP="$(CPP) -I../../inc -I. $(POUND_DEF)" ) ++ ++../external/RSL/RSL/librsl.a : ++ ( cd ../external/RSL/RSL ; make CC="$(CC)" CFLAGS='-DSTUBS' FC="$(FC) $(PROMOTION)" FFLAGS='' MAX_DOMAINS=$(MAX_DOMAINS) MAX_PROC=$(MAX_PROC) LEARN_BCAST=-DLEARN_BCAST stub ) ++ ++# compile these without high optimization to speed compile ++solve_interface.o : solve_interface.F ++shift_domain_em.o : shift_domain_em.F ++module_io_mm5.o : module_io_mm5.F ++module_si_io.o : module_si_io.F ++module_io_wrf.o : module_io_wrf.F ++module_domain.o : module_domain.F ++module_start.o : module_start.F ++module_initialize.o : module_initialize.F ++module_initialize_b_wave.o : module_initialize_b_wave.F ++module_initialize_hill2d_x.o : module_initialize_hill2d_x.F ++module_initialize_quarter_ss.o : module_initialize_quarter_ss.F ++module_initialize_squall2d_x.o : module_initialize_squall2d_x.F ++module_initialize_squall2d_y.o : module_initialize_squall2d_y.F ++module_initialize_real.o : module_initialize_real.F ++module_dm.o : module_dm.F ++start_domain.o : start_domain.F ++mediation_integrate.o : mediation_integrate.F ++module_configure.o : module_configure.F ++ ++solve_interface.o \ ++shift_domain_em.o \ ++module_io_mm5.o module_si_io.o module_io_wrf.o module_domain.o \ ++module_start.o module_initialize.o module_initialize_b_wave.o \ ++module_initialize_hill2d_x.o module_initialize_quarter_ss.o \ ++module_initialize_squall2d_x.o module_initialize_squall2d_y.o \ ++module_configure.o \ ++mediation_integrate.o \ ++module_initialize_real.o module_dm.o start_domain.o : ++ $(RM) $@ ++ $(SED_FTN) $*.F > $*.b ++ $(CPP) -I../inc $(CPPFLAGS) $*.b > $*.f90 ++ $(RM) $*.b ++ $(FC) -c $(PROMOTION) $(FCBASEOPTS) $(MODULE_DIRS) $*.f90 ++ ++ ++########################################################### ++#ARCH PC Linux x86_64 (IA64 and Opteron), GNU Fortran compiler DM-Parallel (RSL_LITE, OpenMPI, Allows nesting) ++# ++DMPARALLEL = 1 ++MAX_PROC = 1024 ++FC = mpif90 ++LD = mpif90 ++CC = mpicc ++SCC = gcc ++SFC = gfortran ++RWORDSIZE = $(NATIVE_RWORDSIZE) ++PROMOTION = ++CFLAGS = -DDM_PARALLEL \ ++ -DMAXDOM_MAKE=$(MAX_DOMAINS) -DMAXPROC_MAKE=$(MAX_PROC) -I../external/RSL_LITE ++FCOPTIM = ++FCDEBUG = ++FCBASEOPTS = -w $(FCDEBUG) -ffree-form -ffree-line-length-none -fconvert=big-endian -frecord-marker=4 ++FCFLAGS = $(FCOPTIM) $(FCBASEOPTS) ++ARCHFLAGS = -DGRIB1 CONFIGURE_GRIB2_FLAG -DINTIO -DDM_PARALLEL \ ++ -DIWORDSIZE=4 -DDWORDSIZE=8 -DRWORDSIZE=$(RWORDSIZE) -DLWORDSIZE=4 CONFIGURE_NETCDF_FLAG \ ++ -DLIMIT_ARGS -DNMM_NEST=$(WRF_NMM_NEST) ++# machine-specific flags needed to link in ESMF library (C++ run-time-library, etc.) ++ESMF_LIB_FLAGS = ++ESMF_IO_LIB = ESMFIOLIB ++ESMF_IO_LIB_EXT = ESMFIOEXTLIB ++INCLUDE_MODULES = -J../main -I../main -I../external/io_netcdf -I../external/io_int -I../external/esmf_time_f90 \ ++ -I../external -I../frame -I../share -I../phys -I../chem -I../inc ++PERL = perl ++REGISTRY = Registry ++LIB = CONFIGURE_NETCDF_LIB_PATH -L../external/RSL_LITE -lrsl_lite \ ++ -L../external/io_grib1 -lio_grib1 \ ++ CONFIGURE_GRIB2_LIB \ ++ -L../external/io_grib_share -lio_grib_share \ ++ -L../external/io_int -lwrfio_int \ ++ ../frame/module_internal_header_util.o ../frame/pack_utils.o -L../external/esmf_time_f90 -lesmf_time ++LDFLAGS = $(FCFLAGS) CONFIGURE_LDFLAGS -fconvert=big-endian ++ENVCOMPDEFS = CONFIGURE_COMPILEFLAGS ++CPP = /lib/cpp -C -P -traditional ++POUND_DEF = -DNO_RRTM_PHYSICS -traditional $(COREDEFS) -DNONSTANDARD_SYSTEM -DCONFIG_BUF_LEN=$(CONFIG_BUF_LEN) -DMAX_DOMAINS_F=$(MAX_DOMAINS) ++CPPFLAGS = -I$(LIBINCLUDE) -C -P $(ARCHFLAGS) -I../external/RSL_LITE -C -P `cat ../inc/dm_comm_cpp_flags` $(ENVCOMPDEFS) $(POUND_DEF) ++AR = ar ru ++M4 = m4 -B 14000 ++RANLIB = ranlib ++NETCDFPATH = CONFIGURE_NETCDF_PATH ++CC_TOOLS = cc ++ ++externals : wrf_ioapi_includes CONFIGURE_WRFIO_NF wrfio_grib_share wrfio_grib1 CONFIGURE_WRFIO_GRIB2 ../external/RSL_LITE/librsl_lite.a wrfio_int gen_comms.c module_dm.F esmf_time ++ ++gen_comms.c : ../external/RSL_LITE/gen_comms.c ++ ( /bin/cp ../tools/gen_comms_warning ../tools/gen_comms.c ; \ ++ cat ../external/RSL_LITE/gen_comms.c >> ../tools/gen_comms.c ) ++ ++module_dm.F : ../external/RSL_LITE/module_dm.F ++ ( /bin/cp module_dm_warning module_dm.F ; cat ../external/RSL_LITE/module_dm.F >> module_dm.F ) ++ ++wrfio_nf : ++ ( cd ../external/io_netcdf ; \ ++ make NETCDFPATH=CONFIGURE_NETCDF_PATH RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(FC) $(PROMOTION) $(FCFLAGS) -w" TRADFLAG="-traditional" ) ++ ++wrfio_grib_share : ++ ( cd ../external/io_grib_share ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib1 : ++ ( cd ../external/io_grib1 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) ++ ++wrfio_grib2 : ++ ( cd ../external/io_grib2 ; \ ++ make CC="$(SCC)" CFLAGS="$(CFLAGS) CONFIGURE_GRIB2_INC" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" FIXED="-Mfixed" archive) ++ ++wrfio_int : ++ ( cd ../external/io_int ; \ ++ make CC="$(CC)" RANLIB="$(RANLIB)" CPP="$(CPP)" SFC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" FC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" all ) ++ ++esmf_time : ++ ( cd ../external/esmf_time_f90 ; \ ++ make FC="$(FC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" RANLIB="$(RANLIB)" CPP="$(CPP) -I../../inc -I. $(POUND_DEF)" ) ++ ++../external/RSL_LITE/librsl_lite.a : ++ ( cd ../external/RSL_LITE ; make CC="$(CC)" FC="$(FC) $(FCFLAGS) $(PROMOTION) -byteswapio" MAX_DOMAINS=$(MAX_DOMAINS) MAX_PROC=$(MAX_PROC) LEARN_BCAST=-DLEARN_BCAST ) ++ ++# compile these without high optimization to speed compile ++solve_interface.o : solve_interface.F ++shift_domain_em.o : shift_domain_em.F ++module_io_mm5.o : module_io_mm5.F ++module_si_io.o : module_si_io.F ++module_io_wrf.o : module_io_wrf.F ++module_domain.o : module_domain.F ++module_start.o : module_start.F ++module_initialize.o : module_initialize.F ++module_initialize_b_wave.o : module_initialize_b_wave.F ++module_initialize_hill2d_x.o : module_initialize_hill2d_x.F ++module_initialize_quarter_ss.o : module_initialize_quarter_ss.F ++module_initialize_squall2d_x.o : module_initialize_squall2d_x.F ++module_initialize_squall2d_y.o : module_initialize_squall2d_y.F ++module_initialize_real.o : module_initialize_real.F ++module_dm.o : module_dm.F ++start_domain.o : start_domain.F ++mediation_integrate.o : mediation_integrate.F ++module_configure.o : module_configure.F ++ ++solve_interface.o \ ++shift_domain_em.o \ ++module_io_mm5.o module_si_io.o module_io_wrf.o module_domain.o \ ++module_start.o module_initialize.o module_initialize_b_wave.o \ ++module_initialize_hill2d_x.o module_initialize_quarter_ss.o \ ++module_initialize_squall2d_x.o module_initialize_squall2d_y.o \ ++module_configure.o \ ++mediation_integrate.o \ ++module_initialize_real.o module_dm.o start_domain.o : ++ $(RM) $@ ++ $(SED_FTN) $*.F > $*.b ++ $(CPP) -I../inc $(CPPFLAGS) $*.b > $*.f90 ++ $(RM) $*.b ++ $(FC) -c $(PROMOTION) $(FCBASEOPTS) $(MODULE_DIRS) $*.f90 ++ + ################################################################## + #ARCH NULL Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/co2clouds_module_initialize_quarter_ss.F =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/co2clouds_module_initialize_quarter_ss.F (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/co2clouds_module_initialize_quarter_ss.F (revision 1185) @@ -0,0 +1,965 @@ +!IDEAL:MODEL_LAYER:INITIALIZATION +! + +! This MODULE holds the routines which are used to perform various initializations +! for the individual domains. + +! This MODULE CONTAINS the following routines: + +! initialize_field_test - 1. Set different fields to different constant +! values. This is only a test. If the correct +! domain is not found (based upon the "id") +! then a fatal error is issued. + +!----------------------------------------------------------------------- + +MODULE module_initialize + + USE module_domain + USE module_io_domain + USE module_state_description + USE module_model_constants + USE module_bc + USE module_timing + USE module_configure + USE module_init_utilities +#ifdef DM_PARALLEL + USE module_dm +#endif + + +CONTAINS + + +!------------------------------------------------------------------- +! this is a wrapper for the solver-specific init_domain routines. +! Also dereferences the grid variables and passes them down as arguments. +! This is crucial, since the lower level routines may do message passing +! and this will get fouled up on machines that insist on passing down +! copies of assumed-shape arrays (by passing down as arguments, the +! data are treated as assumed-size -- ie. f77 -- arrays and the copying +! business is avoided). Fie on the F90 designers. Fie and a pox. + + SUBROUTINE init_domain ( grid ) + + IMPLICIT NONE + + ! Input data. + TYPE (domain), POINTER :: grid + ! Local data. + INTEGER :: dyn_opt + INTEGER :: idum1, idum2 + + CALL nl_get_dyn_opt( 1, dyn_opt ) + + CALL set_scalar_indices_from_config ( head_grid%id , idum1, idum2 ) + + IF ( dyn_opt .eq. 1 & + .or. dyn_opt .eq. 2 & + .or. dyn_opt .eq. 3 & + ) THEN + CALL init_domain_rk( grid & +! +#include +! + ) + + ELSE + WRITE(0,*)' init_domain: unknown or unimplemented dyn_opt = ',dyn_opt + CALL wrf_error_fatal ( ' init_domain: unknown or unimplemented dyn_opt ' ) + ENDIF + + END SUBROUTINE init_domain + +!------------------------------------------------------------------- + + SUBROUTINE init_domain_rk ( grid & +! +# include +! +) + IMPLICIT NONE + + ! Input data. + TYPE (domain), POINTER :: grid + +# include + + TYPE (grid_config_rec_type) :: config_flags + + ! Local data + INTEGER :: & + ids, ide, jds, jde, kds, kde, & + ims, ime, jms, jme, kms, kme, & + its, ite, jts, jte, kts, kte, & + i, j, k + + ! Local data + + INTEGER, PARAMETER :: nl_max = 1000 + REAL, DIMENSION(nl_max) :: zk, p_in, theta, rho, u, v, qv, pd_in + INTEGER :: nl_in + + + INTEGER :: icm,jcm, ii, im1, jj, jm1, loop, error, fid, nxc, nyc + REAL :: u_mean,v_mean, f0, p_surf, p_level, qvf, z_at_v, z_at_u + REAL :: z_scale, xrad, yrad, zrad, rad, delt, cof1, cof2 +! REAL, EXTERNAL :: interp_0 + REAL :: hm + REAL :: pi + +! stuff from original initialization that has been dropped from the Registry + REAL :: vnu, xnu, xnus, dinit0, cbh, p0_temp, t0_temp, zd, zt + REAL :: qvf1, qvf2, pd_surf + INTEGER :: it + real :: thtmp, ptmp, temp(3) + + LOGICAL :: moisture_init + LOGICAL :: stretch_grid, dry_sounding + + INTEGER :: xs , xe , ys , ye + REAL :: mtn_ht + LOGICAL, EXTERNAL :: wrf_dm_on_monitor + +!!MARS + REAL :: lon_input, lat_input, alt_input, tsurf_input +!!MARS + + +#ifdef DM_PARALLEL +# include +#endif + + + SELECT CASE ( model_data_order ) + CASE ( DATA_ORDER_ZXY ) + kds = grid%sd31 ; kde = grid%ed31 ; + ids = grid%sd32 ; ide = grid%ed32 ; + jds = grid%sd33 ; jde = grid%ed33 ; + + kms = grid%sm31 ; kme = grid%em31 ; + ims = grid%sm32 ; ime = grid%em32 ; + jms = grid%sm33 ; jme = grid%em33 ; + + kts = grid%sp31 ; kte = grid%ep31 ; ! note that tile is entire patch + its = grid%sp32 ; ite = grid%ep32 ; ! note that tile is entire patch + jts = grid%sp33 ; jte = grid%ep33 ; ! note that tile is entire patch + CASE ( DATA_ORDER_XYZ ) + ids = grid%sd31 ; ide = grid%ed31 ; + jds = grid%sd32 ; jde = grid%ed32 ; + kds = grid%sd33 ; kde = grid%ed33 ; + + ims = grid%sm31 ; ime = grid%em31 ; + jms = grid%sm32 ; jme = grid%em32 ; + kms = grid%sm33 ; kme = grid%em33 ; + + its = grid%sp31 ; ite = grid%ep31 ; ! note that tile is entire patch + jts = grid%sp32 ; jte = grid%ep32 ; ! note that tile is entire patch + kts = grid%sp33 ; kte = grid%ep33 ; ! note that tile is entire patch + CASE ( DATA_ORDER_XZY ) + ids = grid%sd31 ; ide = grid%ed31 ; + kds = grid%sd32 ; kde = grid%ed32 ; + jds = grid%sd33 ; jde = grid%ed33 ; + + ims = grid%sm31 ; ime = grid%em31 ; + kms = grid%sm32 ; kme = grid%em32 ; + jms = grid%sm33 ; jme = grid%em33 ; + + its = grid%sp31 ; ite = grid%ep31 ; ! note that tile is entire patch + kts = grid%sp32 ; kte = grid%ep32 ; ! note that tile is entire patch + jts = grid%sp33 ; jte = grid%ep33 ; ! note that tile is entire patch + + END SELECT + + + stretch_grid = .true. +! stretch_grid = .false. + + delt = 3. + delt = 10. +! z_scale = .50 + z_scale = .40 + +!!!MARS +!!!MARS +! open(unit=16,file='input_vert',form='formatted',status='old') +! rewind(16) +! read(16,*) delt, z_scale +! write(6,*) 'delt, z_scale are ', delt, z_scale +! close(16) +!!!MARS +!!!MARS + + pi = 2.*asin(1.0) + write(6,*) ' pi is ',pi + nxc = (ide-ids)/2 + nyc = (jde-jds)/2 + + CALL model_to_grid_config_rec ( grid%id , model_config_rec , config_flags ) + +! here we check to see if the boundary conditions are set properly + + CALL boundary_condition_check( config_flags, bdyzone, error, grid%id ) + + moisture_init = .true. + + grid%itimestep=0 + +#ifdef DM_PARALLEL + CALL wrf_dm_bcast_bytes( icm , IWORDSIZE ) + CALL wrf_dm_bcast_bytes( jcm , IWORDSIZE ) +#endif + + CALL nl_set_mminlu(1, ' ') + CALL nl_set_iswater(1,0) + CALL nl_set_cen_lat(1,40.) + CALL nl_set_cen_lon(1,-105.) + CALL nl_set_truelat1(1,0.) + CALL nl_set_truelat2(1,0.) + CALL nl_set_moad_cen_lat (1,0.) + CALL nl_set_stand_lon (1,0.) + CALL nl_set_map_proj(1,0) + + +! here we initialize data we currently is not initialized +! in the input data + + DO j = jts, jte + DO i = its, ite + grid%msft(i,j) = 1. + grid%msfu(i,j) = 1. + grid%msfv(i,j) = 1. + grid%sina(i,j) = 0. + grid%cosa(i,j) = 1. + grid%e(i,j) = 0. + grid%f(i,j) = 0. !! MARS: put coriolis here if needed + + END DO + END DO + + DO j = jts, jte + DO k = kts, kte + DO i = its, ite + grid%em_ww(i,k,j) = 0. + END DO + END DO + END DO + + grid%step_number = 0 + +! set up the grid + + IF (stretch_grid) THEN ! exponential stretch for eta (nearly constant dz) + DO k=1, kde + grid%em_znw(k) = (exp(-(k-1)/float(kde-1)/z_scale) - exp(-1./z_scale))/ & + (1.-exp(-1./z_scale)) + ENDDO + ELSE + DO k=1, kde + grid%em_znw(k) = 1. - float(k-1)/float(kde-1) + ENDDO + ENDIF + +!!MARS +!!MARS + open(unit=12,file='levels',form='formatted',status='old') + rewind(12) + DO k=1, kde + read(12,*) grid%em_znw(k) + write(6,*) 'read level ', k,grid%em_znw(k) + ENDDO + close(12) +!!MARS +!!MARS + + +! print *, 'coucou' +! stop + + DO k=1, kde-1 + grid%em_dnw(k) = grid%em_znw(k+1) - grid%em_znw(k) + grid%em_rdnw(k) = 1./grid%em_dnw(k) + grid%em_znu(k) = 0.5*(grid%em_znw(k+1)+grid%em_znw(k)) + ENDDO + DO k=2, kde-1 + grid%em_dn(k) = 0.5*(grid%em_dnw(k)+grid%em_dnw(k-1)) + grid%em_rdn(k) = 1./grid%em_dn(k) + grid%em_fnp(k) = .5* grid%em_dnw(k )/grid%em_dn(k) + grid%em_fnm(k) = .5* grid%em_dnw(k-1)/grid%em_dn(k) + ENDDO + + cof1 = (2.*grid%em_dn(2)+grid%em_dn(3))/(grid%em_dn(2)+grid%em_dn(3))*grid%em_dnw(1)/grid%em_dn(2) + cof2 = grid%em_dn(2) /(grid%em_dn(2)+grid%em_dn(3))*grid%em_dnw(1)/grid%em_dn(3) + grid%cf1 = grid%em_fnp(2) + cof1 + grid%cf2 = grid%em_fnm(2) - cof1 - cof2 + grid%cf3 = cof2 + + grid%cfn = (.5*grid%em_dnw(kde-1)+grid%em_dn(kde-1))/grid%em_dn(kde-1) + grid%cfn1 = -.5*grid%em_dnw(kde-1)/grid%em_dn(kde-1) + grid%rdx = 1./config_flags%dx + grid%rdy = 1./config_flags%dy + +! get the sounding from the ascii sounding file, first get dry sounding and +! calculate base state + + dry_sounding = .true. + IF ( wrf_dm_on_monitor() ) THEN + write(6,*) ' getting dry sounding for base state ' + + CALL get_sounding( zk, p_in, pd_in, theta, rho, u, v, qv, dry_sounding, nl_max, nl_in ) + ENDIF + CALL wrf_dm_bcast_real( zk , nl_max ) + CALL wrf_dm_bcast_real( p_in , nl_max ) + CALL wrf_dm_bcast_real( pd_in , nl_max ) + CALL wrf_dm_bcast_real( theta , nl_max ) + CALL wrf_dm_bcast_real( rho , nl_max ) + CALL wrf_dm_bcast_real( u , nl_max ) + CALL wrf_dm_bcast_real( v , nl_max ) + CALL wrf_dm_bcast_real( qv , nl_max ) + CALL wrf_dm_bcast_integer ( nl_in , 1 ) + + write(6,*) ' returned from reading sounding, nl_in is ',nl_in + +! find ptop for the desired ztop (ztop is input from the namelist), +! and find surface pressure + + grid%p_top = interp_0( p_in, zk, config_flags%ztop, nl_in ) + +!!MARS +!!MARS + open(unit=14,file='input_coord',form='formatted',status='old') + rewind(14) + read(14,*) lon_input + read(14,*) lat_input + close(14) + write(6,*) ' lon is ',lon_input + write(6,*) ' lat is ',lat_input +!!MARS +!!MARS + +!!MARS +!!MARS + open(unit=18,file='input_more',form='formatted',status='old') + rewind(18) + read(18,*) alt_input, tsurf_input + close(18) + write(6,*) ' alt is ',alt_input + write(6,*) ' tsurf is ',tsurf_input +!!MARS +!!MARS + + DO j=jts,jte + DO i=its,ite +!!MARS + grid%ht(i,j) = alt_input + grid%tsk(i,j) = tsurf_input +!!MARS + grid%xlat(i,j) = lat_input + grid%xlong(i,j) = lon_input + grid%mars_emiss(i,j)=0.95 + grid%mars_cice(i,j)=0. + grid%slpx(i,j) = 0. + grid%slpy(i,j) = 0. +!!MARS + ENDDO + ENDDO + +! print *, 'coucou' +! stop + + + xs=ide/2 -3 + xs=ids -3 + xe=xs + 6 + ys=jde/2 -3 + ye=ys + 6 + mtn_ht = 500 +#ifdef MTN + DO j=max(ys,jds),min(ye,jde-1) + DO i=max(xs,ids),min(xe,ide-1) + grid%ht(i,j) = mtn_ht * 0.25 * & + ( 1. + COS ( 2*pi/(xe-xs) * ( i-xs ) + pi ) ) * & + ( 1. + COS ( 2*pi/(ye-ys) * ( j-ys ) + pi ) ) + ENDDO + ENDDO +#endif +#ifdef EW_RIDGE + DO j=max(ys,jds),min(ye,jde-1) + DO i=ids,ide + grid%ht(i,j) = mtn_ht * 0.50 * & + ( 1. + COS ( 2*pi/(ye-ys) * ( j-ys ) + pi ) ) + ENDDO + ENDDO +#endif +#ifdef NS_RIDGE + DO j=jds,jde + DO i=max(xs,ids),min(xe,ide-1) + grid%ht(i,j) = mtn_ht * 0.50 * & + ( 1. + COS ( 2*pi/(xe-xs) * ( i-xs ) + pi ) ) + ENDDO + ENDDO +#endif + DO j=jts,jte + DO i=its,ite + grid%em_phb(i,1,j) = g * grid%ht(i,j) + grid%em_ph0(i,1,j) = g * grid%ht(i,j) + ENDDO + ENDDO + + DO J = jts, jte + DO I = its, ite + + p_surf = interp_0( p_in, zk, grid%em_phb(i,1,j)/g, nl_in ) + grid%em_mub(i,j) = p_surf-grid%p_top + +! this is dry hydrostatic sounding (base state), so given grid%em_p (coordinate), +! interp theta (from interp) and compute 1/rho from eqn. of state + + DO K = 1, kte-1 + p_level = grid%em_znu(k)*(p_surf - grid%p_top) + grid%p_top + grid%em_pb(i,k,j) = p_level + grid%em_t_init(i,k,j) = interp_0( theta, p_in, p_level, nl_in ) - t0 + grid%em_alb(i,k,j) = (r_d/p1000mb)*(grid%em_t_init(i,k,j)+t0)*(grid%em_pb(i,k,j)/p1000mb)**cvpm + ENDDO + +! calc hydrostatic balance (alternatively we could interp the geopotential from the +! sounding, but this assures that the base state is in exact hydrostatic balance with +! respect to the model eqns. + + DO k = 2,kte + grid%em_phb(i,k,j) = grid%em_phb(i,k-1,j) - grid%em_dnw(k-1)*grid%em_mub(i,j)*grid%em_alb(i,k-1,j) + ENDDO + + ENDDO + ENDDO + + IF ( wrf_dm_on_monitor() ) THEN + write(6,*) ' ptop is ',grid%p_top + write(6,*) ' base state grid%em_mub(1,1), p_surf is ',grid%em_mub(1,1),grid%em_mub(1,1)+grid%p_top + ENDIF + +! calculate full state for each column - this includes moisture. + + write(6,*) ' getting moist sounding for full state ' + dry_sounding = .false. + CALL get_sounding( zk, p_in, pd_in, theta, rho, u, v, qv, dry_sounding, nl_max, nl_in ) + + DO J = jts, min(jde-1,jte) + DO I = its, min(ide-1,ite) + +! At this point grid%p_top is already set. find the DRY mass in the column +! by interpolating the DRY pressure. + + pd_surf = interp_0( pd_in, zk, grid%em_phb(i,1,j)/g, nl_in ) + +! compute the perturbation mass and the full mass + + grid%em_mu_1(i,j) = pd_surf-grid%p_top - grid%em_mub(i,j) + grid%em_mu_2(i,j) = grid%em_mu_1(i,j) + grid%em_mu0(i,j) = grid%em_mu_1(i,j) + grid%em_mub(i,j) + +! given the dry pressure and coordinate system, interp the potential +! temperature and qv + + do k=1,kde-1 + + p_level = grid%em_znu(k)*(pd_surf - grid%p_top) + grid%p_top + + moist(i,k,j,P_QV) = interp_0( qv, pd_in, p_level, nl_in ) + grid%em_t_1(i,k,j) = interp_0( theta, pd_in, p_level, nl_in ) - t0 + grid%em_t_2(i,k,j) = grid%em_t_1(i,k,j) + + + enddo + +! integrate the hydrostatic equation (from the RHS of the bigstep +! vertical momentum equation) down from the top to get grid%em_p. +! first from the top of the model to the top pressure + + k = kte-1 ! top level + + qvf1 = 0.5*(moist(i,k,j,P_QV)+moist(i,k,j,P_QV)) + qvf2 = 1./(1.+qvf1) + qvf1 = qvf1*qvf2 + +! grid%em_p(i,k,j) = - 0.5*grid%em_mu_1(i,j)/grid%em_rdnw(k) + grid%em_p(i,k,j) = - 0.5*(grid%em_mu_1(i,j)+qvf1*grid%em_mub(i,j))/grid%em_rdnw(k)/qvf2 + qvf = 1. + rvovrd*moist(i,k,j,P_QV) + grid%em_alt(i,k,j) = (r_d/p1000mb)*(grid%em_t_1(i,k,j)+t0)*qvf* & + (((grid%em_p(i,k,j)+grid%em_pb(i,k,j))/p1000mb)**cvpm) + grid%em_al(i,k,j) = grid%em_alt(i,k,j) - grid%em_alb(i,k,j) + +! down the column + + do k=kte-2,1,-1 + qvf1 = 0.5*(moist(i,k,j,P_QV)+moist(i,k+1,j,P_QV)) + qvf2 = 1./(1.+qvf1) + qvf1 = qvf1*qvf2 + grid%em_p(i,k,j) = grid%em_p(i,k+1,j) - (grid%em_mu_1(i,j) + qvf1*grid%em_mub(i,j))/qvf2/grid%em_rdn(k+1) + qvf = 1. + rvovrd*moist(i,k,j,P_QV) + grid%em_alt(i,k,j) = (r_d/p1000mb)*(grid%em_t_1(i,k,j)+t0)*qvf* & + (((grid%em_p(i,k,j)+grid%em_pb(i,k,j))/p1000mb)**cvpm) + grid%em_al(i,k,j) = grid%em_alt(i,k,j) - grid%em_alb(i,k,j) + enddo + +! this is the hydrostatic equation used in the model after the +! small timesteps. In the model, grid%em_al (inverse density) +! is computed from the geopotential. + + + grid%em_ph_1(i,1,j) = 0. + DO k = 2,kte + grid%em_ph_1(i,k,j) = grid%em_ph_1(i,k-1,j) - (1./grid%em_rdnw(k-1))*( & + (grid%em_mub(i,j)+grid%em_mu_1(i,j))*grid%em_al(i,k-1,j)+ & + grid%em_mu_1(i,j)*grid%em_alb(i,k-1,j) ) + + grid%em_ph_2(i,k,j) = grid%em_ph_1(i,k,j) + grid%em_ph0(i,k,j) = grid%em_ph_1(i,k,j) + grid%em_phb(i,k,j) + ENDDO + + IF ( wrf_dm_on_monitor() ) THEN + if((i==2) .and. (j==2)) then + write(6,*) ' grid%em_ph_1 calc ',grid%em_ph_1(2,1,2),grid%em_ph_1(2,2,2),& + grid%em_mu_1(2,2)+grid%em_mub(2,2),grid%em_mu_1(2,2), & + grid%em_alb(2,1,2),grid%em_al(1,2,1),grid%em_rdnw(1) + endif + ENDIF + + ENDDO + ENDDO + +!#if 0 + +! thermal perturbation to kick off convection + + write(6,*) ' nxc, nyc for perturbation ',nxc,nyc + write(6,*) ' delt for perturbation ',delt + + DO J = jts, min(jde-1,jte) +! yrad = config_flags%dy*float(j-nyc)/10000. +!! NO FUCKIN' BUBBLE !! +yrad = 0. + DO I = its, min(ide-1,ite) +! xrad = config_flags%dx*float(i-nxc)/10000. +!! NO FUCKIN' BUBBLE !! +xrad = 0. + DO K = 1, kte-1 + +! put in preturbation theta (bubble) and recalc density. note, +! the mass in the column is not changing, so when theta changes, +! we recompute density and geopotential + +! zrad = 0.5*(grid%em_ph_1(i,k,j)+grid%em_ph_1(i,k+1,j) & +! +grid%em_phb(i,k,j)+grid%em_phb(i,k+1,j))/g +! zrad = (zrad-1500.)/1500. +!! NO FUCKIN' BUBBLE !! +zrad = 0. + RAD=SQRT(xrad*xrad+yrad*yrad+zrad*zrad) + IF(RAD <= 1.) THEN +!! NO FUCKIN' BUBBLE !! +! grid%em_t_1(i,k,j)=grid%em_t_1(i,k,j)+delt*COS(.5*PI*RAD)**2 + grid%em_t_2(i,k,j)=grid%em_t_1(i,k,j) + qvf = 1. + rvovrd*moist(i,k,j,P_QV) + grid%em_alt(i,k,j) = (r_d/p1000mb)*(grid%em_t_1(i,k,j)+t0)*qvf* & + (((grid%em_p(i,k,j)+grid%em_pb(i,k,j))/p1000mb)**cvpm) + grid%em_al(i,k,j) = grid%em_alt(i,k,j) - grid%em_alb(i,k,j) + ENDIF + ENDDO + +! rebalance hydrostatically + + DO k = 2,kte + grid%em_ph_1(i,k,j) = grid%em_ph_1(i,k-1,j) - (1./grid%em_rdnw(k-1))*( & + (grid%em_mub(i,j)+grid%em_mu_1(i,j))*grid%em_al(i,k-1,j)+ & + grid%em_mu_1(i,j)*grid%em_alb(i,k-1,j) ) + + grid%em_ph_2(i,k,j) = grid%em_ph_1(i,k,j) + grid%em_ph0(i,k,j) = grid%em_ph_1(i,k,j) + grid%em_phb(i,k,j) + ENDDO + + ENDDO + ENDDO + +!#endif + + IF ( wrf_dm_on_monitor() ) THEN + write(6,*) ' grid%em_mu_1 from comp ', grid%em_mu_1(1,1) + write(6,*) ' full state sounding from comp, ph, grid%em_p, grid%em_al, grid%em_t_1, qv ' + do k=1,kde-1 + write(6,'(i3,1x,5(1x,1pe10.3))') k, grid%em_ph_1(1,k,1)+grid%em_phb(1,k,1), & + grid%em_p(1,k,1)+grid%em_pb(1,k,1), grid%em_alt(1,k,1), & + grid%em_t_1(1,k,1)+t0, moist(1,k,1,P_QV) + enddo + + write(6,*) ' pert state sounding from comp, grid%em_ph_1, pp, alp, grid%em_t_1, qv ' + do k=1,kde-1 + write(6,'(i3,1x,5(1x,1pe10.3))') k, grid%em_ph_1(1,k,1), & + grid%em_p(1,k,1), grid%em_al(1,k,1), & + grid%em_t_1(1,k,1), moist(1,k,1,P_QV) + enddo + ENDIF + +! interp v + + DO J = jts, jte + DO I = its, min(ide-1,ite) + + IF (j == jds) THEN + z_at_v = grid%em_phb(i,1,j)/g + ELSE IF (j == jde) THEN + z_at_v = grid%em_phb(i,1,j-1)/g + ELSE + z_at_v = 0.5*(grid%em_phb(i,1,j)+grid%em_phb(i,1,j-1))/g + END IF + p_surf = interp_0( p_in, zk, z_at_v, nl_in ) + + DO K = 1, kte-1 + p_level = grid%em_znu(k)*(p_surf - grid%p_top) + grid%p_top + grid%em_v_1(i,k,j) = interp_0( v, p_in, p_level, nl_in ) + grid%em_v_2(i,k,j) = grid%em_v_1(i,k,j) + ENDDO + + ENDDO + ENDDO + +! interp u + + DO J = jts, min(jde-1,jte) + DO I = its, ite + + IF (i == ids) THEN + z_at_u = grid%em_phb(i,1,j)/g + ELSE IF (i == ide) THEN + z_at_u = grid%em_phb(i-1,1,j)/g + ELSE + z_at_u = 0.5*(grid%em_phb(i,1,j)+grid%em_phb(i-1,1,j))/g + END IF + + p_surf = interp_0( p_in, zk, z_at_u, nl_in ) + + DO K = 1, kte-1 + p_level = grid%em_znu(k)*(p_surf - grid%p_top) + grid%p_top + grid%em_u_1(i,k,j) = interp_0( u, p_in, p_level, nl_in ) + grid%em_u_2(i,k,j) = grid%em_u_1(i,k,j) + ENDDO + + ENDDO + ENDDO + +! set w + + DO J = jts, min(jde-1,jte) + DO K = kts, kte + DO I = its, min(ide-1,ite) + grid%em_w_1(i,k,j) = 0. + grid%em_w_2(i,k,j) = 0. + ENDDO + ENDDO + ENDDO + +! set a few more things + + DO J = jts, min(jde-1,jte) + DO K = kts, kte-1 + DO I = its, min(ide-1,ite) + grid%h_diabatic(i,k,j) = 0. + ENDDO + ENDDO + ENDDO + + IF ( wrf_dm_on_monitor() ) THEN + DO k=1,kte-1 + grid%em_t_base(k) = grid%em_t_1(1,k,1) + grid%qv_base(k) = moist(1,k,1,P_QV) + grid%u_base(k) = grid%em_u_1(1,k,1) + grid%v_base(k) = grid%em_v_1(1,k,1) + grid%z_base(k) = 0.5*(grid%em_phb(1,k,1)+grid%em_phb(1,k+1,1)+grid%em_ph_1(1,k,1)+grid%em_ph_1(1,k+1,1))/g + ENDDO + ENDIF + CALL wrf_dm_bcast_real( grid%em_t_base , kte ) + CALL wrf_dm_bcast_real( grid%qv_base , kte ) + CALL wrf_dm_bcast_real( grid%u_base , kte ) + CALL wrf_dm_bcast_real( grid%v_base , kte ) + CALL wrf_dm_bcast_real( grid%z_base , kte ) + + DO J = jts, min(jde-1,jte) + DO I = its, min(ide-1,ite) + thtmp = grid%em_t_2(i,1,j)+t0 + ptmp = grid%em_p(i,1,j)+grid%em_pb(i,1,j) + temp(1) = thtmp * (ptmp/p1000mb)**rcp + thtmp = grid%em_t_2(i,2,j)+t0 + ptmp = grid%em_p(i,2,j)+grid%em_pb(i,2,j) + temp(2) = thtmp * (ptmp/p1000mb)**rcp + thtmp = grid%em_t_2(i,3,j)+t0 + ptmp = grid%em_p(i,3,j)+grid%em_pb(i,3,j) + temp(3) = thtmp * (ptmp/p1000mb)**rcp + +!!MARS +! grid%tsk(I,J)=grid%cf1*temp(1)+grid%cf2*temp(2)+grid%cf3*temp(3) + grid%tmn(I,J)=grid%tsk(I,J)-0.5 +!!!MARS +!!TODO: passer la valeur a partir des donnees +!grid%mars_tsoil(I,:,J)=grid%tsk(I,J) +!!!MARS + ENDDO + ENDDO + + END SUBROUTINE init_domain_rk + + SUBROUTINE init_module_initialize + END SUBROUTINE init_module_initialize + +!--------------------------------------------------------------------- + +! test driver for get_sounding +! +! implicit none +! integer n +! parameter(n = 1000) +! real zk(n),p(n),theta(n),rho(n),u(n),v(n),qv(n),pd(n) +! logical dry +! integer nl,k +! +! dry = .false. +! dry = .true. +! call get_sounding( zk, p, pd, theta, rho, u, v, qv, dry, n, nl ) +! write(6,*) ' input levels ',nl +! write(6,*) ' sounding ' +! write(6,*) ' k height(m) press (Pa) pd(Pa) theta (K) den(kg/m^3) u(m/s) v(m/s) qv(g/g) ' +! do k=1,nl +! write(6,'(1x,i3,8(1x,1pe10.3))') k, zk(k), p(k), pd(k), theta(k), rho(k), u(k), v(k), qv(k) +! enddo +! end +! +!--------------------------------------------------------------------------- + + subroutine get_sounding( zk, p, p_dry, theta, rho, & + u, v, qv, dry, nl_max, nl_in ) + implicit none + + integer nl_max, nl_in + real zk(nl_max), p(nl_max), theta(nl_max), rho(nl_max), & + u(nl_max), v(nl_max), qv(nl_max), p_dry(nl_max) + logical dry + + integer n + parameter(n=1000) + logical debug + parameter( debug = .true.) + +! input sounding data + + real p_surf, th_surf, qv_surf + real pi_surf, pi(n) + real h_input(n), th_input(n), qv_input(n), u_input(n), v_input(n) + + !! special MARS + real r_input(n) + real cp_input(n) + real cv_input(n) + real cvpm_input(n) + real pfile_input(n) + real t_input(n) + real rhofile_input(n) + !! special MARS + +! diagnostics + + real rho_surf, p_input(n), rho_input(n) + real pm_input(n) ! this are for full moist sounding + +! local data + + real p1000mb,cv,cp,r,cvpm,g +! parameter (p1000mb = 1.e+05, r = 287, cp = 1003., cv = cp-r, cvpm = -cv/cp, g=9.81 ) + parameter (p1000mb = 610., r = 192., cp = 844.6, cv = cp-r, cvpm = -cv/cp, g=3.72) + integer k, it, nl + real qvf, qvf1, dz + +! first, read the sounding + + call read_sounding( p_surf, th_surf, qv_surf, & + h_input, th_input, qv_input, u_input, v_input, r_input, cp_input, pfile_input, t_input, rhofile_input, n, nl, debug ) + + + !! special MARS + do k=1,nl + cv_input(k) = cp_input(k) - r_input(k) + cvpm_input(k) = - cv_input(k) / cp_input(k) + enddo + !! special MARS + + if(dry) then + do k=1,nl + qv_input(k) = 0. + enddo + endif + + if(debug) write(6,*) ' number of input levels = ',nl + + nl_in = nl + if(nl_in .gt. nl_max ) then + write(6,*) ' too many levels for input arrays ',nl_in,nl_max + call wrf_error_fatal ( ' too many levels for input arrays ' ) + end if + +! compute diagnostics, +! first, convert qv(g/kg) to qv(g/g) + + do k=1,nl + qv_input(k) = 0.001*qv_input(k) + enddo + + p_surf = 100.*p_surf ! convert to pascals + qvf = 1. + rvovrd*qv_input(1) +! rho_surf = 1./((r/p1000mb)*th_surf*qvf*((p_surf/p1000mb)**cvpm)) +! pi_surf = (p_surf/p1000mb)**(r/cp) + rho_surf = 1./((r_input(1)/p1000mb)*th_surf*qvf*((p_surf/p1000mb)**cvpm_input(1))) + pi_surf = (p_surf/p1000mb)**(r_input(1)/cp_input(1)) + + + if(debug) then + write(6,*) ' surface density is ',rho_surf + write(6,*) ' surface pi is ',pi_surf + end if + + +! integrate moist sounding hydrostatically, starting from the +! specified surface pressure +! -> first, integrate from surface to lowest level + + qvf = 1. + rvovrd*qv_input(1) + qvf1 = 1. + qv_input(1) + rho_input(1) = rho_surf + dz = h_input(1) + do it=1,10 +!!MARS MARS + pm_input(1) = p_surf !& +! - 0.5*dz*(rho_surf+rho_input(1))*g*qvf1 +! rho_input(1) = 1./((r/p1000mb)*th_input(1)*qvf*((pm_input(1)/p1000mb)**cvpm)) +!! special MARS + rho_input(1) = 1./((r_input(1)/p1000mb)*th_input(1)*qvf*((pm_input(1)/p1000mb)**cvpm_input(1))) + enddo + +! integrate up the column + + do k=2,nl + rho_input(k) = rho_input(k-1) +! dz = h_input(k)-h_input(k-1) +! special MARS + dz = r_input(k) * t_input(k) * (- p_input(k) + p_input(k-1)) / p_input(k) / g +! dz = - cp_input(k) * (- t_input(k) + t_input(k-1)) / g + qvf1 = 0.5*(2.+(qv_input(k-1)+qv_input(k))) + qvf = 1. + rvovrd*qv_input(k) ! qv is in g/kg here + + do it=1,10 + pm_input(k) = pm_input(k-1) & + - 0.5*dz*(rho_input(k)+rho_input(k-1))*g*qvf1 +! rho_input(k) = 1./((r/p1000mb)*th_input(k)*qvf*((pm_input(k)/p1000mb)**cvpm)) +!! special MARS + rho_input(k) = 1./((r_input(k)/p1000mb)*th_input(k)*qvf*((pm_input(k)/p1000mb)**cvpm_input(k))) + !! + print *, p_input(k), pm_input(k), ((r_input(k)/p1000mb)*th_input(k)*qvf*((pm_input(k)/p1000mb)**cvpm_input(k))), k + enddo + enddo + +! we have the moist sounding + +! next, compute the dry sounding using p at the highest level from the +! moist sounding and integrating down. + + p_input(nl) = pm_input(nl) + + do k=nl-1,1,-1 + dz = h_input(k+1)-h_input(k) + p_input(k) = p_input(k+1) + 0.5*dz*(rho_input(k)+rho_input(k+1))*g + enddo + + + do k=1,nl + + zk(k) = h_input(k) + !p(k) = pm_input(k) + !p_dry(k) = p_input(k) + p(k) = pfile_input(k) + p_dry(k) = pfile_input(k) + theta(k) = th_input(k) + !rho(k) = rho_input(k) + rho(k) = rhofile_input(k) + u(k) = u_input(k) + v(k) = v_input(k) + qv(k) = qv_input(k) + + enddo + + if(debug) then + write(6,*) ' sounding ' + write(6,*) ' k height(m) press (Pa) pd(Pa) theta (K) den(kg/m^3) u(m/s) v(m/s) qv(g/g) ' + do k=1,nl + write(6,'(1x,i3,8(1x,1pe10.3))') k, zk(k), p(k), p_dry(k), theta(k), rho(k), u(k), v(k), qv(k) + enddo + + end if + + end subroutine get_sounding + +!------------------------------------------------------- + + subroutine read_sounding( ps,ts,qvs,h,th,qv,u,v,r,cp,p,t,rho,n,nl,debug ) + implicit none + integer n,nl + real ps,ts,qvs,h(n),th(n),qv(n),u(n),v(n),r(n),cp(n),p(n),t(n),rho(n) + logical end_of_file + logical debug + + integer k + + open(unit=10,file='input_sounding',form='formatted',status='old') + rewind(10) + read(10,*) ps, ts, qvs + if(debug) then + write(6,*) ' input sounding surface parameters ' + write(6,*) ' surface pressure (mb) ',ps + write(6,*) ' surface pot. temp (K) ',ts + write(6,*) ' surface mixing ratio (g/kg) ',qvs + end if + + end_of_file = .false. + k = 0 + + do while (.not. end_of_file) + + read(10,*,end=100) h(k+1), th(k+1), qv(k+1), u(k+1), v(k+1) + k = k+1 + if(debug) write(6,'(1x,i3,5(1x,e10.3))') k, h(k), th(k), qv(k), u(k), v(k) + go to 110 + 100 end_of_file = .true. + 110 continue + enddo + + + !!! special MARS + open(unit=11,file='input_therm',form='formatted',status='old') + rewind(11) + end_of_file = .false. + k = 0 + do while (.not. end_of_file) + + read(11,*,end=101) r(k+1), cp(k+1), p(k+1), rho(k+1), t(k+1) + write(*,*) k, r(k+1), cp(k+1), p(k+1), rho(k+1), t(k+1) + k = k+1 + go to 112 + 101 end_of_file = .true. + 112 continue + enddo + !!! special MARS + + + + nl = k + + close(unit=10,status = 'keep') + + end subroutine read_sounding + +END MODULE module_initialize Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/estimations =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/estimations (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/estimations (revision 1185) @@ -0,0 +1,4 @@ +\item (base) 64 bits, 1 processeur, domaine 50x50x50, résolution 7 km, sans traceurs : 0.4 h/jour. +\item (traceurs) 64 bits, 4 processeurs, domaine 201x201x61, résolution 20km, 2 traceurs : 8 h/jour. +\item (domaines emboîtés) 64 bits, 4 processeurs, 3 domaines 155x155x61, résolutions 36 km, 9 km, 2 km, sans traceurs : 20 h/jour. +\item (LES) 64 bits, 1 processeur, domaine 150x150x71, résolution 100 m, sans traceurs : 11 h/h. Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_ideal.txt =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_ideal.txt (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_ideal.txt (revision 1185) @@ -0,0 +1,4 @@ +jusque 201 OK avec ideal.exe sequentiel +ensuite il vaut mieux utiliser +mpirun -n 4 ideal.exe +le MP_STACK_SIZE est dans le bashrc Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.les =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.les (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.les (revision 1185) @@ -0,0 +1,68 @@ +#################################################################### +# @ job_name = LMD_LES_MARS +# @ job_type = MPICH +# @ total_tasks = 8 +# @ class = AP +# @ environment = real=0;tmpdir=$(job_name).$(jobid).results +#################################################################### +#################################################################### +#################################################################### +# @ output = $(job_name).$(jobid) +# @ error = $(output) +# @ blocking = unlimited +# @ notification = never +# @ queue +grep "^# @" $0 +set -vx +#################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR/${tmpdir} +#################################################################### +ln -sf $LOADL_STEP_INITDIR/*.exe . +ln -sf $LOADL_STEP_INITDIR/*.def . +ln -sf $LOADL_STEP_INITDIR/namelist.* . +ln -sf $LOADL_STEP_INITDIR/wrfi* . ### executer sur la frontale + ### au prealable !!! +#################################################################### +if [[ "${real}" == "0" ]] +then + ln -sf $LOADL_STEP_INITDIR/input_* . +else + ln -sf $LOADL_STEP_INITDIR/wrfb* . +fi +#################################################################### +declare -x MP_STACK_SIZE=640000000 +mpirun wrf.exe +#################################################################### +#################################################################### +#################################################################### + + + +#################################################################### +## @ total_tasks =16 ## ne fonctionne pas avec 241 points +## @ total_tasks =20 ## semble mieux fonctionner... bof +## @ total_tasks =24 ## il y a un souci a l'interface + + +#################################################################### +#################################################################### +###############NE MARCHE PAS AVEC LES GRANDS DOMAINES ############## +#################################################################### +#################################################################### +#echo "#! /bin/bash" > script +#echo "ulimit -s unlimited" >> script +#echo "./ideal.exe" >> script +#echo "mv rsl.out.0000 ideal_rsl.out.0000" >> script +#chmod 755 script +#mpirun ./script +#\rm script +#################################################################### +##declare -x MP_STACK_SIZE=64000000 +##mpirun ideal.exe ### PB: ideal.exe est limite a un moment +## ### mais mpirun ideal.exe donne des instabilites +##ulimit -s unlimited ### ne resout pas PB... ok 145 pas ok 201 +##./ideal.exe +##mpirun "ulimit -s unlimited ; ./ideal.exe" +##mv rsl.out.0000 ideal_rsl.out.0000 +#################################################################### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso (revision 1185) @@ -0,0 +1,30 @@ +# @ job_name=LMD_MM_MARS +# @ output = $(job_name).$(jobid) +# @ error = $(output) +# @ job_type = MPICH +# @ total_tasks=32 +## @ total_tasks=24 +## @ total_tasks =16 +## @ total_tasks =8 +## @ total_tasks =4 +# @ blocking = unlimited +# @ class =AP +## @ class =BP +## @ class =CP +# @ environment=tmpdir=$(job_name).$(jobid).results +# @ notification = never +# @ queue +grep "^# @" $0 +set -vx +#################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR/${tmpdir} +#################################################################### +ln -sf $LOADL_STEP_INITDIR/*.def . +ln -sf $LOADL_STEP_INITDIR/wrfb* . +ln -sf $LOADL_STEP_INITDIR/wrfi* . +ln -sf $LOADL_STEP_INITDIR/namelist* . +declare -x MP_STACK_SIZE=64000000 +ln -sf /san/home/spiga/trunk//mesoscale/LMD_MM_MARS/newphys_mpifort_64/wrf_x241_y241_z61_d1_t4_p32.exe wrf.exe +mpirun wrf.exe +#################################################################### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.mpich =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.mpich (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.mpich (revision 1185) @@ -0,0 +1,32 @@ +# @ job_name=GW_MARS +# @ output = $(job_name).$(jobid) +# @ error = $(output) +# @ job_type = MPICH +# @ total_tasks =24 +## @ total_tasks =16 +## @ total_tasks =8 +## @ total_tasks =4 +# @ blocking = unlimited +# @ class =AP +## @ class =BP +## @ class =CP +# @ environment=tmpdir=$(job_name).$(jobid).results +# @ queue +grep "^# @" $0 +set -vx +#################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +\rm GW_MARS_last +cd $LOADL_STEP_INITDIR +ln -sf ${tmpdir} GW_MARS_last +cd ${tmpdir} +#################################################################### +cp $LOADL_STEP_INITDIR/*.def . +cp $LOADL_STEP_INITDIR/wrf* . +cp $LOADL_STEP_INITDIR/namelist* . +cp $LOADL_STEP_INITDIR/levels . +cp $LOADL_STEP_INITDIR/input_* . +cp $LOADL_STEP_INITDIR/ze_hill . +declare -x MP_STACK_SIZE=64000000 +mpirun /san/home/spiga/trunk//mesoscale//LMD_MM_MARS/mpifort_64/wrf_x121_y121_z181_d1_t1_p24.exe +#################################################################### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.seq =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.seq (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnome.sav/gnome_launch.meso.ideal.seq (revision 1185) @@ -0,0 +1,24 @@ +# @ job_name=GW_MARS +# @ output = $(job_name).$(jobid) +# @ error = $(output) +# @ job_type = serial +# @ class = A +# @ environment=tmpdir=$(job_name).$(jobid).results +## @ resources = ConsumableCpus(1)ConsumableMemory(24gb) +# @ queue +grep "^# @" $0 +set -vx +#################################################################### +mkdir $LOADL_STEP_INITDIR/${tmpdir} +cd $LOADL_STEP_INITDIR/${tmpdir} +#################################################################### +ln -sf $LOADL_STEP_INITDIR/*.exe . +cp $LOADL_STEP_INITDIR/*.def . +#cp $LOADL_STEP_INITDIR/wrf* . +cp $LOADL_STEP_INITDIR/namelist* . +cp $LOADL_STEP_INITDIR/levels . +cp $LOADL_STEP_INITDIR/input_* . +cp $LOADL_STEP_INITDIR/ze_hill . +./ideal.exe > log_ideal +./wrf.exe > log_wrf +#################################################################### Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure (revision 1185) @@ -0,0 +1,4 @@ + 14. PC Linux x86_64 (IA64 and AMD Opteron), GNU Fortran compiler, DM-Parallel (RSL, MPICH, Allows nesting) + 15. PC Linux x86_64 (IA64 and Opteron), GNU Fortran compiler (Single-threaded, RSL, Allows nesting) + 16. PC Linux x86_64 (IA64 and Opteron), GNU Fortran compiler DM-Parallel (RSL_LITE, OpenMPI, Allows nesting) + Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure.wrf =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure.wrf (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64.configure.wrf (revision 1185) @@ -0,0 +1,256 @@ +# configure.wrf +# +# This file was automatically generated by the configure script in the +# top level directory. You may make changes to the settings in this +# file but be aware they will be overwritten each time you run configure. +# Ordinarily, it is necessary to run configure once, when the code is +# first installed. +# +# To permanently change options, change the settings for your platform +# in the file arch/configure.defaults then rerun configure. +# +SHELL = /bin/sh +DEVTOP = `pwd` +LIBINCLUDE = . +.SUFFIXES: .F .i .o .f90 .c + +#### Get core settings from environment (set in compile script) +#### Note to add a core, this has to be added to. + +COREDEFS = -DEM_CORE=$(WRF_EM_CORE) \ + -DNMM_CORE=$(WRF_NMM_CORE) -DNMM_MAX_DIM=2600 \ + -DCOAMPS_CORE=$(WRF_COAMPS_CORE) \ + -DEXP_CORE=$(WRF_EXP_CORE) + +#### Single location for defining total number of domains. You need +#### at least 1 + 2*(number of total nests). For example, 1 coarse +#### grid + three fine grids = 1 + 2(3) = 7, so MAX_DOMAINS=7. + +MAX_DOMAINS = 21 + +#### DM buffer length for the configuration flags. + +CONFIG_BUF_LEN = 16384 + + +############################################################################## +#### The settings in this section are defaults that may be overridden by the +#### architecture-specific settings in the next section. +############################################################################## + +############################################################################## +#### NOTE: Do not modify these default values here. To override these +#### values, make changes after "Architecture specific settings". +############################################################################## + +#### Native size (in bytes) of Fortran REAL data type on this architecture #### +#### Note: to change real wordsize (for example, to promote REALs from +#### 4-byte to 8-byte), modify the definition of RWORDSIZE in the +#### section following "Architecture specific settings". Do not +#### change NATIVE_RWORDSIZE as is it architecture-specific. +NATIVE_RWORDSIZE = 4 + +#### Default sed command and script for Fortran source files #### +SED_FTN = sed -f ../arch/standard.sed + +# Hack to work around $(PWD) not changing during OSF1 build. +# $(IO_GRIB_SHARE_DIR) is reset during the OSF1 build only. +IO_GRIB_SHARE_DIR = + +#### ESMF switches #### +#### These are set up by Config.pl #### +# switch to use separately installed ESMF library for coupling: 1==true +ESMF_COUPLING = 0 +# select dependences on module_utility.o +ESMF_MOD_DEPENDENCE = ../external/esmf_time_f90/module_utility.o +# select -I options for separately installed ESMF library, if present +ESMF_MOD_INC = +# select -I options for external/io_esmf vs. external/esmf_time_f90 +ESMF_IO_INC = -I../external/esmf_time_f90 +# select cpp token for external/io_esmf vs. external/esmf_time_f90 +ESMF_IO_DEFS = +# select build target for external/io_esmf vs. external/esmf_time_f90 +ESMF_TARGET = esmf_time + +############################################################################## + + +#### Architecture specific settings #### + +# Settings for PC Linux x86_64 (IA64 and Opteron), GNU Fortran compiler (Single-threaded, RSL, Allows nesting) +# +DMPARALLEL = 1 +MAX_PROC = 1024 +FC = gfortran +LD = gfortran +CC = gcc -DFSEEKO64_OK +SCC = $(CC) +SFC = gfortran +RWORDSIZE = $(NATIVE_RWORDSIZE) +FCOPTIM = #-O3 -funroll-loops #-ftree-vectorize -ftree-loop-linear +FCDEBUG = #-g +FCBASEOPTS = -fno-second-underscore -w $(FCDEBUG) -ffree-form -ffree-line-length-none -fconvert=big-endian -frecord-marker=4 +FCFLAGS = $(FCOPTIM) $(FCBASEOPTS) +ARCHFLAGS = -DGRIB1 -DINTIO -DWRF_RSL_IO -DRSL -DDM_PARALLEL \ + -DIWORDSIZE=4 -DDWORDSIZE=8 -DRWORDSIZE=$(RWORDSIZE) -DLWORDSIZE=4 -DNETCDF \ + -DLIMIT_ARGS -DSTUBMPI +PROMOTION = +CFLAGS = -DLANDREADSTUB -DDM_PARALLEL -DWRF_RSL_IO \ + -DMAXDOM_MAKE=$(MAX_DOMAINS) -DMAXPROC_MAKE=$(MAX_PROC) -I../external/RSL/RSL $(ARCHFLAGS) +# machine-specific flags needed to link in ESMF library (C++ run-time-library, etc.) +ESMF_LIB_FLAGS = +ESMF_IO_LIB = -L../external/esmf_time_f90 -lesmf_time +ESMF_IO_LIB_EXT = -L../../external/esmf_time_f90 -lesmf_time +INCLUDE_MODULES = -J../main -I../main -I../external/io_netcdf -I../external/io_int -I../external/esmf_time_f90 \ + -I../external -I../frame -I../share -I../phys -I../chem -I../inc +PERL = perl +REGISTRY = Registry +LIB = -L../external/io_netcdf -lwrfio_nf -L/donnees/aslmd/MODELES/NETCDF/gfortran_64/netcdf-3.6.1//lib -lnetcdf -L../external/RSL/RSL -lrsl \ + -L../external/io_grib1 -lio_grib1 \ + -L../external/io_int -lwrfio_int \ + -L../external/io_grib_share -lio_grib_share \ + ../frame/module_internal_header_util.o ../frame/pack_utils.o -L../external/esmf_time_f90 -lesmf_time -L../mars_lmd/libo -llmd /donnees/aslmd/MODELES/NETCDF/gfortran_64/netcdf-3.6.1//lib/libnetcdf.a +LDFLAGS = $(FCFLAGS) -fconvert=big-endian +ENVCOMPDEFS = +WRF_CHEM = 0 +CPP = /lib/cpp -C -P -traditional +POUND_DEF = -DNO_RRTM_PHYSICS $(COREDEFS) -DNONSTANDARD_SYSTEM -DCONFIG_BUF_LEN=$(CONFIG_BUF_LEN) -DMAX_DOMAINS_F=$(MAX_DOMAINS) +CPPFLAGS = -I$(LIBINCLUDE) -C -P $(ARCHFLAGS) -I../external/RSL/RSL -C -P `cat ../inc/dm_comm_cpp_flags` $(ENVCOMPDEFS) $(POUND_DEF) +AR = ar ru +M4 = m4 -B 14000 +RANLIB = ranlib +NETCDFPATH = /donnees/aslmd/MODELES/NETCDF/gfortran_64/netcdf-3.6.1/ +CC_TOOLS = cc + +externals : wrf_ioapi_includes wrfio_nf wrfio_grib_share wrfio_grib1 ../external/RSL/RSL/librsl.a wrfio_int gen_comms.c module_dm.F esmf_time + +gen_comms.c : ../external/RSL/gen_comms.c + ( /bin/cp ../tools/gen_comms_warning ../tools/gen_comms.c ; \ + cat ../external/RSL/gen_comms.c >> ../tools/gen_comms.c ) + +module_dm.F : ../external/RSL/module_dm.F + ( /bin/cp module_dm_warning module_dm.F ; cat ../external/RSL/module_dm.F >> module_dm.F ) + +wrfio_nf : + ( cd ../external/io_netcdf ; \ + make NETCDFPATH=/donnees/aslmd/MODELES/NETCDF/gfortran_64/netcdf-3.6.1/ RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(FC) $(PROMOTION) $(FCFLAGS) -w" TRADFLAG="-traditional" ) + +wrfio_grib_share : + ( cd ../external/io_grib_share ; \ + make CC="$(SCC)" CFLAGS="$(CFLAGS)" RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) + +wrfio_grib1 : + ( cd ../external/io_grib1 ; \ + make CC="$(SCC)" CFLAGS="$(CFLAGS)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" archive) + +wrfio_grib2 : + ( cd ../external/io_grib2 ; \ + make CC="$(SCC)" CFLAGS="$(CFLAGS) " RM="$(RM)" RANLIB="$(RANLIB)" CPP="$(CPP)" FC="$(SFC) $(PROMOTION) -I. $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" FIXED="-Mfixed" archive) + +wrfio_int : + ( cd ../external/io_int ; \ + make CC="$(CC)" RANLIB="$(RANLIB)" CPP="$(CPP)" SFC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" FC="$(SFC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS) -w" TRADFLAG="-traditional" all ) + +esmf_time : + ( cd ../external/esmf_time_f90 ; \ + make FC="$(FC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS)" RANLIB="$(RANLIB)" CPP="$(CPP) -I../../inc -I. $(POUND_DEF)" ) + +../external/RSL/RSL/librsl.a : + ( cd ../external/RSL/RSL ; make CC="$(CC)" CFLAGS='-DSTUBS' FC="$(FC) $(PROMOTION)" FFLAGS='' MAX_DOMAINS=$(MAX_DOMAINS) MAX_PROC=$(MAX_PROC) LEARN_BCAST=-DLEARN_BCAST stub ) + +# compile these without high optimization to speed compile +solve_interface.o : solve_interface.F +shift_domain_em.o : shift_domain_em.F +module_io_mm5.o : module_io_mm5.F +module_si_io.o : module_si_io.F +module_io_wrf.o : module_io_wrf.F +module_domain.o : module_domain.F +module_start.o : module_start.F +module_initialize.o : module_initialize.F +module_initialize_b_wave.o : module_initialize_b_wave.F +module_initialize_hill2d_x.o : module_initialize_hill2d_x.F +module_initialize_quarter_ss.o : module_initialize_quarter_ss.F +module_initialize_squall2d_x.o : module_initialize_squall2d_x.F +module_initialize_squall2d_y.o : module_initialize_squall2d_y.F +module_initialize_real.o : module_initialize_real.F +module_dm.o : module_dm.F +start_domain.o : start_domain.F +mediation_integrate.o : mediation_integrate.F +module_configure.o : module_configure.F + +solve_interface.o \ +shift_domain_em.o \ +module_io_mm5.o module_si_io.o module_io_wrf.o module_domain.o \ +module_start.o module_initialize.o module_initialize_b_wave.o \ +module_initialize_hill2d_x.o module_initialize_quarter_ss.o \ +module_initialize_squall2d_x.o module_initialize_squall2d_y.o \ +module_configure.o \ +mediation_integrate.o \ +module_initialize_real.o module_dm.o start_domain.o : + $(RM) $@ + $(SED_FTN) $*.F > $*.b + $(CPP) -I../inc $(CPPFLAGS) $*.b > $*.f90 + $(RM) $*.b + $(FC) -c $(PROMOTION) $(FCBASEOPTS) $(MODULE_DIRS) $*.f90 + + +########################################################### + +# These sub-directory builds are identical across all architectures +wrfio_esmf : + ( cd ../external/io_esmf ; \ + make FC="$(FC) $(PROMOTION) $(FCDEBUG) $(FCBASEOPTS) $(ESMF_MOD_INC)" RANLIB="$(RANLIB)" CPP="$(CPP) $(POUND_DEF) " ) + + +# +# Macros, these should be generic for all machines + +LN = ln -sf +MAKE = make -i -r +RM = rm -f + +# There is probably no reason to modify these rules + +wrf_ioapi_includes : + ( cd ../external/ioapi_share ; \ + $(MAKE) NATIVE_RWORDSIZE="$(NATIVE_RWORDSIZE)" RWORDSIZE="$(RWORDSIZE)" ) + +.F.i: + $(RM) $@ + $(CPP) -I../inc $(CPPFLAGS) $*.F > $@ + mv $*.i $(DEVTOP)/pick/$*.f90 + cp $*.F $(DEVTOP)/pick + +.F.o: + $(RM) $@ + $(SED_FTN) $*.F > $*.b + $(CPP) -I../inc $(CPPFLAGS) $*.b > $*.f90 + $(RM) $*.b + $(FC) -c $(FCFLAGS) $(MODULE_DIRS) $(PROMOTION) $(FCSUFFIX) $*.f90 +#PGI if [ ! -e $@ ] ; then \ +#PGI sleep 10 ; $(FC) -c $(FCFLAGS) $(MODULE_DIRS) $(PROMOTION) $(FCSUFFIX) $*.f90 ; \ +#PGI fi +#PGI if [ ! -e $@ ] ; then \ +#PGI sleep 10 ; $(FC) -c $(FCFLAGS) $(MODULE_DIRS) $(PROMOTION) $(FCSUFFIX) $*.f90 ; \ +#PGI fi + +.F.f90: + $(RM) $@ + $(SED_FTN) $*.F > $*.b + $(CPP) -I../inc $(CPPFLAGS) $*.b > $@ + $(RM) $*.b + +.f90.o: + $(RM) $@ + $(FC) -c $(FCFLAGS) $(PROMOTION) $(FCSUFFIX) $*.f90 +#PGI if [ ! -e $@ ] ; then \ +#PGI sleep 10 ; $(FC) -c $(FCFLAGS) $(PROMOTION) $(FCSUFFIX) $*.f90 ; \ +#PGI fi +#PGI if [ ! -e $@ ] ; then \ +#PGI sleep 10 ; $(FC) -c $(FCFLAGS) $(PROMOTION) $(FCSUFFIX) $*.f90 ; \ +#PGI fi + +.c.o: + $(RM) $@ + $(CC) -c $(CFLAGS) $*.c + Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64_single.txt =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64_single.txt (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/gnu_64_single.txt (revision 1185) @@ -0,0 +1,53 @@ +gnu simple en 64 bits sur huchard, etc... +------------------- + +choix 99 de makemeso pour le moment + + +**** segmentation fault dans le cas de la montagne 2D **** + + + +autre possibilite: essayer les librairies NETCDF4.0 d'Ehouarn + + +essai sans O3 +NB: on a juste a modifier SRC/WRFV2/configure.defaults.gfortran64 +>>> NE REGLE PAS LE PB + + +IL Y A UN PROBLEME AVEC LES VALEURS CALCULEES... +ou: Floating exception sur auric + + +../mars_lmd/libo/liblmd.a(physdem1.o): In function `physdem1_': +physdem1.F:(.text+0x6bb): undefined reference to `nf_put_var_real_' +physdem1.F:(.text+0x74a): undefined reference to `nf_put_var_real_' +physdem1.F:(.text+0x7d9): undefined reference to `nf_put_var_real_' +physdem1.F:(.text+0x868): undefined reference to `nf_put_var_real_' +physdem1.F:(.text+0x8f7): undefined reference to `nf_put_var_real_' +../mars_lmd/libo/liblmd.a(physdem1.o):physdem1.F:(.text+0x986): more undefined +references to `nf_put_var_real_' follow +../mars_lmd/libo/liblmd.a(meso_readtesassim.o): In function +`meso_readtesassim_': +meso_readtesassim.F90:(.text+0x1224): undefined reference to +`nf_get_var_real_' +meso_readtesassim.F90:(.text+0x14cc): undefined reference to +`nf_get_var_real_' +meso_readtesassim.F90:(.text+0x179a): undefined reference to +`nf_get_var_real_' +meso_readtesassim.F90:(.text+0x17dd): undefined reference to +`nf_get_var_real_' +../mars_lmd/libo/liblmd.a(readtesassim.o): In function `readtesassim_': +readtesassim.F90:(.text+0x11ec): undefined reference to `nf_get_var_real_' +../mars_lmd/libo/liblmd.a(readtesassim.o):readtesassim.F90:(.text+0x121b): +more undefined references to `nf_get_var_real_' follow +collect2: ld returned 1 exit status +make[1]: [em_wrf] Error 1 (ignored) + + +>>>> corrige en ajoutant dans makemeso + + sed s+"-L../mars_lmd/libo -llmd"+"-L../mars_lmd/libo +-llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah +configure.wrf ;; Index: trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/meso_readtesassim.F90 =================================================================== --- trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/meso_readtesassim.F90 (revision 1185) +++ trunk/MESOSCALE/LMD_MM_MARS/SRC/DEV/WORK/meso_readtesassim.F90 (revision 1185) @@ -0,0 +1,297 @@ +subroutine readtesassim(ngrid,nlayer,zday,pplev,tauref) + +! Reading of the dust assimilation file + +implicit none + +#include "dimensions.h" +#include "dimphys.h" +#include "comgeomfi.h" +#include "netcdf.inc" +#include "datafile.h" +#include "callkeys.h" + +integer, intent(in) :: ngrid,nlayer +real, intent(in) :: zday ! date in martian days +real, dimension(ngrid,nlayer+1), intent(in) :: pplev +real, dimension(ngrid), intent(out) :: tauref + +! Local variables + +real :: realday +integer nid,nvarid,ierr +integer ltloop,lsloop,iloop,jloop,varloop,ig +real, dimension(2) :: taubuf +real tau1(4),tau +real alt(4) +integer latp(4),lonp(4) +integer yinf,ysup,xinf,xsup,tinf,tsup +real latinf,latsup,loninf,lonsup +real latintmp,lonintmp,latdeg,londeg +real colat,dlat,dlon,colattmp +logical, save :: firstcall=.true. +logical :: timeflag +real,save :: radeg,pi +integer :: timedim,londim,latdim +real, dimension(:), allocatable, save :: lat,lon,time +real, dimension(:,:,:), allocatable, save :: tautes +integer, save :: timelen,lonlen,latlen +character(len=33) :: filename + +realday=mod(zday,669.) + +if (firstcall) then + firstcall=.false. + + pi=acos(-1.) + radeg=180/pi + + ! assimilated dust file: (NB: iaervar is a common in "callkeys.h") + ! iaervar=4 means read dust_tes.nc file + ! iaervar=24 means read dust_tes_MY24.nc file + ! iaervar=25 means read dust_tes_MY25.nc file + ! iaervar=26 means read dust_tes_MY26.nc file + if (iaervar.eq.4) then + filename="dust_tes.nc" + else if (iaervar.eq.24) then + filename="dust_tes_MY24.nc" ! which currently is identical to dust_tes.nc + else if (iaervar.eq.25) then + filename="dust_tes_MY25.nc" + else if (iaervar.eq.26) then + filename="dust_tes_MY26.nc" + else if ((iaervar.eq.6).or.(iaervar.eq.7)) then ! modified TES scenario + filename="modifieddust.nc" + endif + + ! Note: datafile() is defined in "datafile.h" + ierr=NF_OPEN(trim(datafile)//"/"//trim(filename),NF_NOWRITE,nid) + IF (ierr.NE.NF_NOERR) THEN + write(*,*)'Problem opening ',trim(filename),' (in phymars/readtesassim.F90)' + write(*,*)'It should be in :',trim(datafile),'/' + write(*,*)'1) You can change this directory address in ' + write(*,*)' file phymars/datafile.h' + write(*,*)'2) If necessary, dust.nc (and other datafiles)' + write(*,*)' can be obtained online on:' + write(*,*)' http://www.lmd.jussieu.fr/~forget/datagcm/datafile' + CALL ABORT + ENDIF + + ierr=NF_INQ_DIMID(nid,"Time",timedim) + ierr=NF_INQ_DIMLEN(nid,timedim,timelen) + ierr=NF_INQ_DIMID(nid,"latitude",latdim) + ierr=NF_INQ_DIMLEN(nid,latdim,latlen) + ierr=NF_INQ_DIMID(nid,"longitude",londim) + ierr=NF_INQ_DIMLEN(nid,londim,lonlen) + + allocate(tautes(lonlen,latlen,timelen)) + allocate(lat(latlen), lon(lonlen), time(timelen)) + + ierr = NF_INQ_VARID (nid, "dustop",nvarid) +#ifdef NC_DOUBLE + ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tautes) +#else + ierr = NF_GET_VAR_REAL(nid, nvarid, tautes) +#endif + IF (ierr .NE. NF_NOERR) THEN + PRINT*, "Error: Readtesassim not found" + stop + ENDIF + ierr = NF_INQ_VARID (nid, "Time",nvarid) +#ifdef NC_DOUBLE + ierr = NF_GET_VAR_DOUBLE(nid, nvarid, time) +#else + ierr = NF_GET_VAR_REAL(nid, nvarid, time) +#endif + IF (ierr .NE. NF_NOERR) THEN + PRINT*, "Error: Readtesassim