Ignore:
Timestamp:
May 11, 2014, 3:56:22 PM (11 years ago)
Author:
aslmd
Message:

MESOSCALE. Adapted the mesoscale model to previous commit to LMD GCM physics where number of scatterers is no longer set at compiling time.

Location:
trunk/MESOSCALE/LMD_MM_MARS
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars0

    r1187 r1247  
    148148photochem  = .false.
    149149
     150# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     151naerkind = 2
     152
     153
    150154## Thermospheric options (relevant if tracer=T) :
    151155##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars11_radact

    r1187 r1247  
    142142photochem  = .false.
    143143
     144# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     145naerkind = 2
     146
     147
    144148## Thermospheric options (relevant if tracer=T) :
    145149##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars12_radact_scav

    r1115 r1247  
    140140photochem  = .false.
    141141
     142# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     143naerkind = 2
     144
     145
    142146## Thermospheric options (relevant if tracer=T) :
    143147##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars1_radact

    r1187 r1247  
    145145photochem  = .false.
    146146
     147# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     148naerkind = 2
     149
     150
    147151## Thermospheric options (relevant if tracer=T) :
    148152##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars2

    r1187 r1247  
    147147photochem  = .false.
    148148
     149# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     150naerkind = 2
     151
     152
    149153## Thermospheric options (relevant if tracer=T) :
    150154##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/DEF/newphys_THARSIS_WATER/callphys.def.mars3

    r1187 r1247  
    144144photochem  = .false.
    145145
     146# SCATTERERS: set number of scatterers. must be compliant with preceding options.
     147naerkind = 2
     148
     149
    146150## Thermospheric options (relevant if tracer=T) :
    147151##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  • trunk/MESOSCALE/LMD_MM_MARS/SIMU/runmeso

    r1045 r1247  
    191191echo '** Domain(s)             >>> '${nest}
    192192echo '** Tracer(s)             >>> '${tracers}
    193 
    194    ############ SCATTERERS SCATTERERS ######
    195    if [[ ${newphys} -eq 1 ]]
    196    then
    197      #### CHANGE THIS IF YOU ADD NEW SCENARIOS
    198      case ${mars_scenario} in
    199        0) scat=1 ;;
    200        1) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
    201        2) scat=1 ;;
    202        3) scat=1 ;;
    203       11) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
    204       12) echo "!! How many scatterers are you using : 1,2,..." ; read scat ;;
    205       20) scat=1 ;;
    206        *) scat=1 ;;
    207      esac
    208      echo '** Scatterers            >>> '${scat}
    209    else
    210      #### USELESS here BUT COMPATIBILITY
    211      scat=1
    212    fi
    213    #########################################
    214 
    215193echo '** Grid points           >>> '${nx}" x "${ny}" x "${nz}
    216194echo '** Horizontal resolution >>> '${dx}
    … …  
    379357 cd ${simu_folder}
    380358     ## corresponding executables
    381     wrfexe=${work_folder}/wrf_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
    382     realexe=${work_folder}/real_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}_s${scat}.exe
     359    wrfexe=${work_folder}/wrf_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}.exe
     360    realexe=${work_folder}/real_x${nx}_y${ny}_z${nz}_d${nest}_t${tracers}_p${numproc}.exe
    383361 if [[ ${fresh} -eq 0 && -f ${wrfexe} ]]
    384362 then
    … …  
    404382   touch ${simu_folder}/makemeso_input
    405383   echo ${compilo} >> ${simu_folder}/makemeso_input
    406    if [[ "${compilo}" = "3" || "${compilo}" = "5" ]]
     384   if [[ ${newphys} -eq 0 ]]
     385   then
     386     if [[ "${compilo}" = "3" || "${compilo}" = "5" ]]
    407387     then
    408388         echo ${numproc} >> ${simu_folder}/makemeso_input
    409    fi
    410    echo ${nx} >> ${simu_folder}/makemeso_input
    411    echo ${ny} >> ${simu_folder}/makemeso_input
     389     fi
     390     echo ${nx} >> ${simu_folder}/makemeso_input
     391     echo ${ny} >> ${simu_folder}/makemeso_input
     392   fi
    412393   echo ${nz} >> ${simu_folder}/makemeso_input
    413394   echo ${nest} >> ${simu_folder}/makemeso_input
    414    echo ${tracers} >> ${simu_folder}/makemeso_input
    415    if [[ ${newphys} -eq 1 ]]
    416    then
    417      echo ${scat} >> ${simu_folder}/makemeso_input
    418      #echo 'n' >> ${simu_folder}/makemeso_input
     395   if [[ ${newphys} -eq 0 ]]
     396   then
     397     echo ${tracers} >> ${simu_folder}/makemeso_input
    419398   fi
    420399   #

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/SCRIPTS/NETCDF/gfortran_netcdf4

    r331 r1247  
    11#! /bin/bash
    22
     3ze_netcdf=netcdf-4.0.1
     4#ze_netcdf=netcdf-3.6.1
    35
    4 ##wget http://www.lmd.jussieu.fr/~hourdin/PATCH_LMDZ4OR_G95/netcdf-3.6.1.tar.gz
    5 ##gunzip netcdf-3.6.1.tar.gz
    6 #tar xvf netcdf-3.6.1.tar
    7 #cd netcdf-3.6.1
    8 
    9 cp /donnees/emlmd/Downloads/Software/netcdf-4.0.1.tar.gz .
    10 gunzip netcdf-4.0.1.tar.gz
    11 tar xvf netcdf-4.0.1.tar
    12 cd netcdf-4.0.1
    13 
    14 #export CC=gcc
    15 #export CXX=g++
    16 #export CPP='gcc -E'
    17 #export FC=gfortran
    18 #export F77=gfortran
    19 #export F90=gfortran
    20 #export CFLAGS='-O2 -DpgiFortran'
    21 #export CXXFLAGS='-O2 -DpgiFortran'
    22 #export CPPFLAGS='-DpgiFortran'
    23 #export FFLAGS='-O2 -fno-second-underscore'
    24 #export FCFLAGS='-O2 -fno-second-underscore'
    25 #export F90FLAGS='-O2 -fno-second-underscore'
     6wget ftp://ftp.unidata.ucar.edu/pub/netcdf/$ze_netcdf.tar.gz
     7tar xzvf $ze_netcdf.tar.gz
     8\rm $ze_netcdf.tar.gz
    269
    2710export FC=gfortran
    … …  
    3518export CXXFLAGS="-O2"
    3619
    37 mkdir /donnees/aslmd/MODELES/MESOSCALE_DEV/NETCDF/netcdf64-4.0.1_gfortran
    38 ./configure --prefix=/donnees/aslmd/MODELES/MESOSCALE_DEV/NETCDF/netcdf64-4.0.1_gfortran #> configure.log 2>&1
    39 make #> make.log 2>&1
    40 make test #> make_test.log 2>&1
    41 make install #> make_install.log 2>&1
     20mv $ze_netcdf gfortran_$ze_netcdf
     21cd gfortran_$ze_netcdf
     22PREFIX=$PWD
     23if [[ "$ze_netcdf" == "netcdf-3.6.1" ]]
     24then
     25      cd src
     26fi
    4227
    43 #cd src
    44 ## Compilation
    45 #./configure
    46 #make check
    47 #make install
    48 
     28./configure --prefix=${PREFIX} #--disable-cxx
     29make > make.log 2>&1
     30make test > make.log 2>&1
     31make install > make.log 2>&1

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

    r1236 r1247  
    584584
    585585  !! Read callphys.def !!
    586   call conf_phys(nq)
     586  call conf_phys(ngrid,nlayer,nq)
    587587
    588588#endif

  • trunk/MESOSCALE/LMD_MM_MARS/makemeso

    r1243 r1247  
    207207                tra=1   
    208208              fi
    209             fi
    210             # scatterers: a specific stuff for newphys
    211             if [[ "${phys}" == "newphys_" ]]
    212             then
    213               echo "How many types of scatterer are you using ? 1 [dust] or 2 [dust + water ice]" ; read scat
    214             else
    215               scat=1 ## dummy stuff
    216209            fi
    217210         fi
    … …  
    367360  if [[ "${phys}" == "newphys_" ]]
    368361  then
    369    echo you have ${scat} scatterers
    370362   echo number of processors can be changed without recompiling
    371363   echo x points can be changed without recompiling
    … …  
    398390   echo ${tra} >> last
    399391  fi
    400   echo ${scat} >> last
    401   #echo ${answer} >> last
    402392
    403393#------------
    … …  
    713703         ln -sf ../dyn3d/gr_fi_dyn.F .  ## dommage, a corriger
    714704     ###############################
    715 \rm scatterers.h
    716 cat << EOF > scatterers.h
    717 !-----------------------------------------------------------------------
    718 !   INCLUDE 'scatterers.h'
    719 !
    720 ! Number of kind of tracer radiative properties
    721 ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    722 ! (ex: naerkind=1 if you use one dust mode without ice ...)
    723 ! (ex: naerkind=2 if you use one dust mode and active ice ...)
    724 !-----------------------------------------------------------------------
    725 
    726       integer, parameter :: naerkind=${scat}
    727 
    728 !-----------------------------------------------------------------------
    729 EOF
    730705     ###############################
    731706     ###############################
    … …  
    789764  # save a copy
    790765  #--------------
    791   cp -f libo/liblmd.a libo/liblmd.a_${lon}_${lat}_${level}_${dom}_${tra}_${scat}
     766  cp -f libo/liblmd.a libo/liblmd.a_${lon}_${lat}_${level}_${dom}_${tra}
    792767  echo '****************************************'
    793768
    … …  
    849824                        echo 'Looks good ! wrf.exe is here...'
    850825                      fi
    851                     cp -f real.exe ../../real_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
    852                     cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
     826                    cp -f real.exe ../../real_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
     827                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
    853828                    cd ..  ;;
    854829    'ideal')        #mkdir 'test/em_quarter_ss' 2> /dev/null
    … …  
    865840                        echo 'Looks good ! wrf.exe is here...'
    866841                      fi
    867                     cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
    868                     cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
     842                    cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
     843                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
    869844                    cd ..  ;;
    870845    'les')          echo '>>> YOUR CONFIG IS : '${config}
    … …  
    881856                        echo 'Looks good ! wrf.exe is here...'
    882857                      fi
    883                     cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
    884                     cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.exe
     858                    cp -f ideal.exe ../../ideal_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
     859                    cp -f wrf.exe ../../wrf_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.exe
    885860                    cd ..  ;;
    886861    *)              echo not supported... please use ; echo ideal les ; exit ;;
    … …  
    894869  echo 'done.'
    895870  echo '****************************************'
    896   cp last ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}
     871  cp last ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}
    897872  mv last ../../
    898   svn info ../../ > ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}_s${scat}.svn.info
     873  svn info ../../ > ../makemeso_x${lon}_y${lat}_z${level}_d${dom}_t${tra}_p${numproc}.svn.info
    899874
    900875  # add here specific messages
Note: See TracChangeset for help on using the changeset viewer.