Changeset 343
- Timestamp:
- Jun 24, 2009 12:59:09 PM (15 years ago)
- Location:
- palm/trunk
- Files:
-
- 12 edited
-
DOC/app/chapter_3.5.4.html (modified) (1 diff)
-
SCRIPTS/.mrun.config.imuk (modified) (7 diffs)
-
SCRIPTS/mrun (modified) (6 diffs)
-
SCRIPTS/subjob (modified) (5 diffs)
-
SOURCE/CURRENT_MODIFICATIONS (modified) (2 diffs)
-
SOURCE/check_parameters.f90 (modified) (1 diff)
-
SOURCE/data_output_2d.f90 (modified) (3 diffs)
-
SOURCE/flow_statistics.f90 (modified) (3 diffs)
-
SOURCE/modules.f90 (modified) (4 diffs)
-
SOURCE/time_integration.f90 (modified) (2 diffs)
-
SOURCE/timestep.f90 (modified) (2 diffs)
-
SOURCE/user_data_output_2d.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_3.5.4.html
r212 r343 275 275 ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;"> 276 276 grid = 'zu'</span><br> <br>The <span style="font-family: monospace;">ELSE</span> case is 277 only needed in case that output of time-averaged data is requested.<br><br><br></li><li>For 277 only needed in case that output of time-averaged data is requested.<br><br> 278 As a special case, xy cross section output can also be defined at one single level at height <span style="font-family: monospace;">k=nzb+1</span> on the u-grid. This features is useful for output of surface data (e.g. heat fluxes). In this case, the corresponding 2d data has to be resorted to the array <span style="font-family: monospace;">local_pf(i,j,nzb+1).</span> In addition to this, the grid in <a href="chapter_3.5.1.html#user_define_netcdf_grid"><span style="font-family: monospace;">user_define_netcdf_grid</span></a> as well as in <a href="chapter_3.5.1.html#user_data_output_2d"><span style="font-family: monospace;">user_data_output_2d</span></a> must be set to <span style="font-family: monospace;">grid = 'zu1'</span>.<br><br> 279 <span style="font-family: monospace;"> CASE 280 ( 'u2_xy' )</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 281 IF ( av == 0 ) THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 282 DO i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 283 DO j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 284 local_pf(i,j,nzb+1) = u2(j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 285 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 286 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 287 ELSE</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 288 DO i = nxl-1, nxr+1</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 289 DO j = nys-1, nyn+1</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 290 local_pf(i,j,nzb+1) = u2_av(j,i)</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 291 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 292 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 293 ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;"> 294 grid = 'zu1'</span><br><span style="font-family: monospace;"> 295 two_d = .TRUE.</span><br><br> 296 Note that <span style="font-family: monospace;">two_d = .TRUE.</span> is necessary for output of a 2d data slice.<br><br> 297 <br></li><li>For 278 298 output of 3d volume data, the gridpoint data of the quantity has to be 279 299 resorted to array <span style="font-family: monospace;">local_pf</span> -
palm/trunk/SCRIPTS/.mrun.config.imuk
r289 r343 21 21 %host_identifier b01* lcsgib 22 22 %host_identifier gallego lcmuk 23 %host_identifier cs* nech 24 %host_identifier hreg*-en0 ibmh 25 %host_identifier hanni*-en0 ibmh 26 %host_identifier breg*-en0 ibmb 27 %host_identifier berni*-en0 ibmb 28 %host_identifier tgg* lctit 29 # 30 # The following blocks contain all information for compiling the PALM code 31 # and for generating and running the PALM executable using MPI. Replace all 32 # required paths (given in <>) by the respective paths valid on your host. 33 # Please adjust compilernames, cpp-options, compiler-options etc. 34 # appropriately for your own system(s), and please delete or comment out 35 # blocks for unused systems. 36 # 37 # The following 'lcmuk parallel' block assumes PALM compilation and execution 38 # on host gallego. Please note that other hosts within the IMUK linux cluster, 39 # e.g. maestro or bora, may have require different settings. 40 # 41 #%remote_username <replace by your linux cluster username> lcmuk parallel 23 %host_identifier nid00* lcxt4 24 # 25 %remote_username maronga lcmuk parallel 42 26 %tmp_user_catalog /tmp lcmuk parallel 43 27 %compiler_name mpif90 lcmuk parallel … … 50 34 #%hostfile $base_directory/.hostfile lcmuk parallel 51 35 # 52 #%remote_username <replace by your linux cluster username>lcmuk36 #%remote_username maronga lcmuk 53 37 %tmp_user_catalog /tmp lcmuk 54 38 %compiler_name ifort lcmuk … … 59 43 %lopts -axW:-cpp:-r8:-nbs:-Vaxlib lcmuk 60 44 # 61 #%remote_username <replace by your linux cluster username>lcmuk trace45 %remote_username maronga lcmuk trace 62 46 %tmp_user_catalog /tmp lcmuk trace 63 47 %compiler_name ifort lcmuk trace … … 69 53 # 70 54 %sgi_feature ice lcsgih parallel 71 #%remote_username <replace by your HLRN username>lcsgih parallel55 %remote_username nikmaron lcsgih parallel 72 56 %compiler_name mpif90 lcsgih parallel 73 57 %compiler_name_ser ifort lcsgih parallel … … 78 62 %fopts -g:-w:-xT:-O3:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-convert:little_endian lcsgih parallel 79 63 %lopts -g:-w:-xT:-O3:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-Vaxlib lcsgih parallel 80 #%tmp_data_catalog /gfs1/work/ <replace by your HLRN username>/palm_restart_datalcsgih parallel81 #%tmp_user_catalog /gfs1/tmp/ <replace by your HLRN username>lcsgih parallel64 #%tmp_data_catalog /gfs1/work/nikmaron/palm_restart_data lcsgih parallel 65 #%tmp_user_catalog /gfs1/tmp/nikmaron lcsgih parallel 82 66 # 83 67 %sgi_feature ice lcsgih parallel debug 68 %remote_username nikmaron lcsgih parallel debug 84 69 %totalview true lcsgih parallel debug 85 #%remote_username <replace by your HLRN username> lcsgih parallel debug86 70 %compiler_name mpif90 lcsgih parallel debug 87 71 %compiler_name_ser ifort lcsgih parallel debug … … 92 76 %fopts -C:-fpe0:-debug:-traceback:-g:-w:-xT:-O0:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-convert:little_endian lcsgih parallel debug 93 77 %lopts -C:-fpe0:-debug:-traceback:-g:-w:-xT:-O0:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-Vaxlib lcsgih parallel debug 94 #%tmp_data_catalog /gfs1/work/<replace by your HLRN username>/palm_restart_data lcsgih parallel debug 95 #%tmp_user_catalog /gfs1/tmp/<replace by your HLRN username> lcsgih parallel debug 96 # 97 #%remote_username <replace by your tsubame username> lctit parallel 98 #%tmp_user_catalog /work/<replace by your tsubame username> lctit parallel 99 #%tmp_data_catalog /work/<replace by your tsubame username>/palm_restart_data lctit parallel 100 %cpp_options -Mpreprocess:-DMPI_REAL=MPI_DOUBLE_PRECISION:-DMPI_2REAL=MPI_2DOUBLE_PRECISION:-D__netcdf lctit parallel 101 %netcdf_inc -I:/home2/usr5/mkanda/netcdf-3.6.2/include lctit parallel 102 %netcdf_lib -L/home2/usr5/mkanda/netcdf-3.6.2/lib:-lnetcdf lctit parallel 103 %compiler_name_ser pgf95 lctit parallel 104 %compiler_name mpif90 lctit parallel 105 %fopts -r8:-fast:-fastsse lctit parallel 106 %lopts -r8:-fast:-fastsse lctit parallel 107 # 108 #%remote_username <replace by your HLRN username> ibmh parallel 109 %compiler_name mpxlf95_r ibmh parallel 110 %compiler_name_ser xlf95 ibmh parallel 111 %cpp_options -qsuffix=cpp=f90:-WF,-DMPI_REAL=MPI_DOUBLE_PRECISION,-D__netcdf=__netcdf,-D__netcdf_64bit=__netcdf_64bit ibmh parallel 112 %dvrp_inc -I/home/h/nikevita/lib_stable/optimize ibmh parallel 113 %dvrp_lib -L/home/h/nikevita/lib_stable/optimize:-lDVRP2:-lftp ibmh parallel 114 %netcdf_inc -I:/aws/dataformats/netcdf-3.6.0-p1/64-32/include-64 ibmh parallel 115 %netcdf_lib -L/aws/dataformats/netcdf-3.6.0-p1/64-32/lib:-lnetcdf ibmh parallel 116 %fopts -O3:-g:-qfloat=nomaf:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape ibmh parallel 117 %lopts -O3:-g:-qfloat=nomaf:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape:-lesslsmp ibmh parallel 118 %memory 1630 ibmh parallel 119 %cpumax 1000 ibmh parallel 120 %numprocs 4 ibmh parallel 121 #%tmp_data_catalog /fastfs/work/<replace by your HLRN username>/palm_restart_data ibmh parallel 122 #%email_notification <replace by your email addres> ibmh parallel 123 # 124 #%remote_username <replace by your HLRN username> ibmh parallel debug 125 %compiler_name mpxlf95_r ibmh parallel debug 126 %compiler_name_ser xlf95 ibmh parallel debug 127 %cpp_options -qsuffix=cpp=f90:-WF,-DMPI_REAL=MPI_DOUBLE_PRECISION,-D__netcdf=__netcdf,-D__netcdf_64bit=__netcdf_64bit ibmh parallel debug 128 %netcdf_inc -I:/aws/dataformats/netcdf-3.6.0-p1/64-32/include-64 ibmh parallel debug 129 %netcdf_lib -L/aws/dataformats/netcdf-3.6.0-p1/64-32/lib:-lnetcdf ibmh parallel debug 130 %fopts -qnoopt:-g:-C:-qinitauto=FFFFFFFF:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qnosave:-qnoescape:-qflttrap=overflow::zerodivide::invalid::enable:-qsigtrap ibmh parallel debug 131 %lopts -qnoopt:-g:-C:-qinitauto=FFFFFFFF:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qnosave:-qnoescape:-qflttrap=overflow::zerodivide::invalid::enable:-qsigtrap:-lesslsmp ibmh parallel debug 132 %memory 1000 ibmh parallel debug 133 %cpumax 1000 ibmh parallel debug 134 %numprocs 4 ibmh parallel debug 135 #%tmp_data_catalog /fastfs/work/<replace by your HLRN username>/palm_restart_data ibmh parallel debug 136 # 137 #%remote_username <replace by your HLRN username> ibmb parallel 138 %compiler_name mpxlf95_r ibmb parallel 139 %compiler_name_ser xlf95 ibmb parallel 140 %cpp_options -qsuffix=cpp=f90:-WF,-DMPI_REAL=MPI_DOUBLE_PRECISION,-D__netcdf=__netcdf,-D__netcdf_64bit=__netcdf_64bit ibmb parallel 141 %netcdf_inc -I:/aws/dataformats/netcdf-3.6.0-p1/64-32/include-64 ibmb parallel 142 %netcdf_lib -L/aws/dataformats/netcdf-3.6.0-p1/64-32/lib:-lnetcdf ibmb parallel 143 %fopts -O3:-g:-qfloat=nomaf:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape ibmb parallel 144 %lopts -O3:-g:-qfloat=nomaf:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape:-lesslsmp ibmb parallel 145 %memory 1630 ibmb parallel 146 %cpumax 1000 ibmb parallel 147 %numprocs 4 ibmb parallel 148 #%tmp_data_catalog /fastfs/work/<replace by your HLRN username>/palm_restart_data ibmb parallel 149 # 150 #%remote_username <replace by your DKRZ username> nech parallel 151 %compiler_name sxmpif90 nech parallel 152 %compiler_name_ser sxf90 nech parallel 153 %cpp_options -Ep:-DMPI_REAL=MPI_REAL8:-DMPI_2REAL=MPI_2REAL8:-DSCFFT=DZFFT:-DCSFFT=ZDFFT:-D__netcdf:-D__netcdf_64bit nech parallel 154 %netcdf_inc -I:/pool/SX-6/netcdf/netcdf-3.6.0-p1/include nech parallel 155 %netcdf_lib -L/pool/SX-6/netcdf/netcdf-3.6.0-p1/lib:-lnetcdf nech parallel 156 %fopts -C:hopt:-Wf:\'-init:stack=nan:-init:heap=nan:-A:idbl4:-pvctl:fullmsg\':-ftrace:-P:stack:-pi:auto:-pi:line=1000 nech parallel 157 %lopts -C:hopt:-Wf:\'-init:stack=nan:-init:heap=nan:-A:idbl4\':-ftrace:-P:stack:-pi:auto:-pi:line=1000:-L/SX/opt/MathKeisan/MK1_6/lib/:-l:fft nech parallel 158 %memory 4000 nech parallel 159 %cpumax 1000 nech parallel 160 %numprocs 4 nech parallel 161 %tmp_data_catalog $WRKSHR/palm_restart_data nech parallel 162 #%email_notification <replace by your email addres> nech parallel 163 # 164 #%remote_username <replace by your NEC-username> neck parallel 165 %compiler_name sxmpif90 neck parallel 166 %compiler_name_ser sxf90 neck parallel 167 %cpp_options -Ep:-DMPI_REAL=MPI_REAL8:-DMPI_2REAL=MPI_2REAL8:-DSCFFT=DZFFT:-DCSFFT=ZDFFT:-D__netcdf:-D__netcdf_64bit neck parallel 168 %netcdf_inc -I:/home/DSRC/NC/tatuyama/pub/netcdf-3.6.0-p1/include neck parallel 169 %netcdf_lib -L/home/DSRC/NC/tatuyama/pub/netcdf-3.6.0-p1/lib:-lnetcdf neck parallel 170 %dvrp_inc -I/home/DSRC/NC/tatuyama/pub/dvrp neck parallel 171 %dvrp_lib -L/home/DSRC/NC/tatuyama/pub/dvrp:-lDVRP2:-lftp neck parallel 172 %fopts -C:hopt:-Wf:\'-init:stack=nan:-init:heap=nan:-A:idbl4:-pvctl:fullmsg\':-ftrace:-P:stack:-pi:auto:-pi:line=1000 neck parallel 173 %lopts -C:hopt:-Wf:\'-init:stack=nan:-init:heap=nan:-A:idbl4\':-ftrace:-P:stack:-pi:auto:-pi:line=1000:-L/SX/opt/MathKeisan/MK1_6/lib/:-l:fft neck parallel 174 %memory 4000 neck parallel 175 %cpumax 1000 neck parallel 176 %numprocs 4 neck parallel 177 #%tmp_user_catalog <replace by your temporary directory> neck parallel 178 #%tmp_data_catalog <replace by your temp. dir.>/palm_restart_data neck parallel 179 # 180 %source_path $base_directory/SOURCE ibmy parallel 181 %compiler_name mpxlf95_r ibmy parallel 182 %compiler_name_ser xlf95 ibmy parallel 183 %fopts -O3:-g:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape ibmy parallel 184 %lopts -O3:-g:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qtune=pwr4:-qarch=pwr4:-qnosave:-qnoescape ibmy parallel 185 %memory 1000 ibmy parallel 186 %cpumax 1000 ibmy parallel 187 %numprocs 4 ibmy parallel 188 %remote_username raasch ibmy parallel 189 %tmp_data_catalog /usr1/users/raasch/work ibmy parallel 190 %tmp_user_catalog /usr1/users/raasch/work ibmy parallel 191 %email_notification raasch@muk.uni-hannover.de ibmy parallel 192 # 193 %compiler_name mpxlf95_r ibmy parallel trace 194 %compiler_name_ser xlf95 ibmy parallel trace 195 %cpp_options -Ep:-DMPI_REAL=MPI_REAL8:-DMPI_2REAL=MPI_2REAL8:-DSCFFT=DZFFT:-DCSFFT=ZDFFT:-D__netcdf:-D__netcdf_64bit nech parallel 196 %fopts -qnoopt:-g:-C:-qinitauto=FFFFFFFF:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qnosave:-qnoescape:-qflttrap=overflow::zerodivide::invalid::enable:-qsigtrap ibmy parallel trace 197 %lopts -qnoopt:-g:-C:-qinitauto=FFFFFFFF:-qrealsize=8:-Q:-q64:-qmaxmem=-1:-qnosave:-qnoescape:-qflttrap=overflow::zerodivide::invalid::enable:-qsigtrap ibmy parallel trace 198 %memory 1000 ibmy parallel trace 199 %cpumax 1000 ibmy parallel trace 200 %numprocs 4 ibmy parallel trace 201 %remote_username raasch ibmy parallel trace 202 %tmp_data_catalog /usr1/users/raasch/work ibmy parallel trace 203 %tmp_user_catalog /usr1/users/raasch/work ibmy parallel trace 204 %email_notification raasch@muk.uni-hannover.de ibmy parallel trace 78 #%tmp_data_catalog /gfs1/work/nikmaron/palm_restart_data lcsgih parallel debug 79 #%tmp_user_catalog /gfs1/tmp/nikmaron lcsgih parallel debug 80 # 81 %sgi_feature ice lcsgib parallel 82 %remote_username nikmaron lcsgib parallel 83 %compiler_name mpif90 lcsgib parallel 84 %compiler_name_ser ifort lcsgib parallel 85 %mpilib mvapich2/1.2.0-intel lcsgib parallel 86 %cpp_options -DMPI_REAL=MPI_DOUBLE_PRECISION:-DMPI_2REAL=MPI_2DOUBLE_PRECISION:-D__netcdf:-D__netcdf_64bit lcsgib parallel 87 %netcdf_inc -I:/sw/dataformats/netcdf/3.6.2/include lcsgib parallel 88 %netcdf_lib -L/sw/dataformats/netcdf/3.6.2/lib:-lnetcdf:-lnetcdff lcsgib parallel 89 90 %fopts -g:-w:-xT:-O3:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-convert:little_endian lcsgib parallel 91 %lopts -g:-w:-xT:-O3:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-Vaxlib lcsgib parallel 92 #%tmp_data_catalog /gfs1/work/nikmaron/palm_restart_data lcsgib parallel 93 #%tmp_user_catalog /gfs1/tmp/nikmaron lcsgib parallel 94 # 95 %sgi_feature ice lcsgib parallel debug 96 %remote_username nikmaron lcsgib parallel debug 97 %compiler_name mpif90 lcsgib parallel debug 98 %compiler_name_ser ifort lcsgib parallel debug 99 %mpilib mvapich2/1.2.0-intel lcsgib parallel debug 100 %cpp_options -DMPI_REAL=MPI_DOUBLE_PRECISION:-DMPI_2REAL=MPI_2DOUBLE_PRECISION:-D__netcdf:-D__netcdf_64bit lcsgib parallel debug 101 %netcdf_inc -I:/sw/dataformats/netcdf/3.6.2/include lcsgib parallel debug 102 %netcdf_lib -L/sw/dataformats/netcdf/3.6.2/lib:-lnetcdf:-lnetcdff lcsgib parallel debug 103 %fopts -C:-fpe0:-debug:-traceback:-g:-w:-xT:-O0:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-convert:little_endian lcsgib parallel debug 104 %lopts -C:-fpe0:-debug:-traceback:-g:-w:-xT:-O0:-cpp:-r8:-ftz:-fno-alias:-no-prec-div:-no-prec-sqrt:-ip:-nbs:-Vaxlib lcsgib parallel debug 105 #%tmp_data_catalog /gfs1/work/nikmaron/palm_restart_data lcsgib parallel debug 106 #%tmp_user_catalog /gfs1/tmp/nikmaron lcsgib parallel debug 107 108 # 109 %remote_username bjornm lcxt4 parallel 110 %tmp_user_catalog /work/bjornm lcxt4 parallel 111 %tmp_data_catalog /work/bjornm/PALM_RESTART_DATA lcxt4 parallel 112 #%job_catalog /work/bjornm/job_queue lcxt4 parallel 113 %compiler_name ftn lcxt4 parallel 114 %compiler_name_ser ftn lcxt4 parallel 115 %cpp_options -DMPI_REAL=MPI_DOUBLE_PRECISION:-DMPI_2REAL=MPI_2DOUBLE_PRECISION:-D__netcdf:-D__intel_openmp_bug lcxt4 parallel 116 %netcdf_inc -I/work/apps/netcdf-3.6.2-cnl-pgi/include lcxt4 parallel 117 %netcdf_lib -L/work/apps/netcdf-3.6.2-cnl-pgi/lib:-lnetcdf lcxt4 parallel 118 %fopts -fastsse:-Mpreprocess:-r8:-Mnobackslash lcxt4 parallel 119 %lopts -fastsse:-Mpreprocess:-r8:-Mnobackslash lcxt4 parallel 120 #%fopts -Mbounds:-fastsse:-Mpreprocess:-r8:-Mnobackslash lcxt4 parallel 121 #%lopts -Mbounds:-fastsse:-Mpreprocess:-r8:-Mnobackslash lcxt4 parallel 122 123 %email_notification maronga@muk.uni-hannover.de lcxt4 parallel 205 124 # 206 125 %write_binary true restart … … 264 183 # 265 184 DATA_1D_PR_NETCDF out:loc:tr pr#:prf $base_data/$fname/OUTPUT _pr nc 185 #DATA_1D_PR_NETCDF out:loc pr#:prf $base_data/$fname/OUTPUT _pr nc 266 186 DATA_1D_SP_NETCDF out:loc:tr sp#:spf $base_data/$fname/OUTPUT _sp nc 267 187 DATA_1D_TS_NETCDF out:loc:tr ts#:tsf $base_data/$fname/OUTPUT _ts nc -
palm/trunk/SCRIPTS/mrun
r331 r343 181 181 # 16/04/09 - Siggi - new option -y for precursor (uncoupled) ocean runs 182 182 # to be followed by a coupled atmosphere-ocean run 183 # 21/04/09 - Siggi - adjustments for new IBM at DKRZ, which is now ibmh 183 # 21/04/09 - Siggi - adjustments for new IBM at DKRZ, which is now ibmh 184 # 24/06/09 - BjornM - adjustments for coupled/precursor runs on ibmy 184 185 185 186 … … 2956 2957 if [[ -n $numprocs ]] 2957 2958 then 2958 2959 2959 # RUNNING THE PROGRAM ON PARALLEL MACHINES 2960 2960 if [[ $(echo $host | cut -c1-3) = ibm ]] … … 2976 2976 [[ $node_usage = not_shared ]] && export MP_SINGLE_THREAD=yes 2977 2977 fi 2978 2979 2978 if [[ "$LOADLBATCH" = yes ]] 2980 2979 then … … 2984 2983 /opt/optibm/HPM_2_4_1/bin/hpmcount a.out $ROPTS 2985 2984 else 2986 poe ./a.out $ROPTS 2985 if [[ $run_coupled_model = false ]] 2986 then 2987 if [[ "$ocean_file_appendix" = true ]] 2988 then 2989 echo "precursor_ocean" > runfile_atmos 2990 else 2991 echo "precursor_atmos" > runfile_atmos 2992 fi 2993 else 2994 (( iia = $numprocs_atmos / $threads_per_task )) 2995 (( iio = $numprocs_ocean / $threads_per_task )) 2996 printf "\n coupled run ($iia atmosphere, $iio ocean)" 2997 printf "\n using $coupled_mode coupling" 2998 printf "\n\n" 2999 echo "coupled_run $iia $iio" > runfile_atmos 3000 fi 3001 poe ./a.out $ROPTS < runfile_atmos 2987 3002 fi 2988 3003 else … … 3004 3019 fi 3005 3020 export MP_HOSTFILE=hostfile 3006 3007 3021 if [[ $run_coupled_model = false ]] 3008 3022 then … … 3021 3035 echo "coupled_run $iia $iio" > runfile_atmos 3022 3036 fi 3023 3024 3037 ./a.out -procs $tasks_per_node $ROPTS < runfile_atmos 3025 3038 -
palm/trunk/SCRIPTS/subjob
r326 r343 99 99 # 21/04/09 - Siggi - adjustments for new IBM at DKRZ, which is now ibmh 100 100 # 18/05/09 - Siggi - Settings for serial jobs on lcsgi changed 101 101 # 24/06/09 - BjornM- adjustments for lcxt4 (loading modules manually) 102 102 103 103 … … 676 676 #PBS -l pmem=${memory}mb 677 677 #PBS -m abe 678 #PBS -M bjorn.maronga@student.uib.no678 #PBS -M maronga@muk.uni-hannover.de 679 679 #PBS -o $remote_dayfile 680 680 #PBS -j oe … … 692 692 #PBS -l pmem=${memory}mb 693 693 #PBS -m abe 694 #PBS -M bjorn.maronga@student.uib.no694 #PBS -M maronga@muk.uni-hannover.de 695 695 #PBS -o $remote_dayfile 696 696 #PBS -j oe … … 761 761 #PBS -l mppnppn=${tasks_per_node} 762 762 #PBS -m abe 763 #PBS -M bjorn.maronga@student.uib.no763 #PBS -M maronga@muk.uni-hannover.de 764 764 #PBS -o $remote_dayfile 765 766 . /opt/modules/default/init/ksh 767 module load modules/3.1.6 768 module load PrgEnv-pgi/2.1.56HD 769 module load Base-opts/2.1.56HD 770 module load moab/5.2.5-s13441 771 module load torque/2.3.4-snap.200809221601 765 772 766 773 %%END%% … … 777 784 #PBS -l pmem=${memory}mb 778 785 #PBS -m abe 779 #PBS -M bjorn.maronga@student.uib.no786 #PBS -M maronga@muk.uni-hannover.de 780 787 #PBS -o $remote_dayfile 788 789 . /opt/modules/default/init/ksh 790 module load modules/3.1.6 791 module load PrgEnv-pgi/2.1.56HD 792 module load Base-opts/2.1.56HD 793 module load moab/5.2.5-s13441 794 module load torque/2.3.4-snap.200809221601 781 795 782 796 %%END%% -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r336 r343 37 37 write_var_list) 38 38 39 check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, modules, mrun, package_parin, parin, prognostic_equations, read_var_list, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_dvrp, user_header, user_init_grid, write_var_list 39 Option for user defined 2D data output in xy cross sections at z=nzb+1 40 (data_output_2d, user_data_output_2d) 41 42 check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, modules, mrun, package_parin, parin, prognostic_equations, read_var_list, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_2d, user_data_output_dvrp, user_header, user_init_grid, write_var_list 40 43 41 44 New: set_particle_attributes, set_slicer_attributes_dvrp … … 91 94 Default value of grid_matching changed to strict. 92 95 93 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list 96 Adjustments for runs on lcxt4 (necessary due to an software update on CRAY) and 97 for coupled runs on ibmy (mrun, subjob) 98 99 advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, mrun, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, subjob, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list 94 100 95 101 -
palm/trunk/SOURCE/check_parameters.f90
r332 r343 391 391 ELSEIF ( TRIM( coupling_mode ) == 'uncoupled' .AND. & 392 392 TRIM( coupling_char ) == '_O' ) THEN 393 393 394 ! 394 395 !-- Check whether an (uncoupled) atmospheric run has been declared as an 395 396 !-- ocean run (this setting is done via mrun-option -y) 397 396 398 message_string = 'ocean = .F. does not allow coupling_char = "' // & 397 399 TRIM( coupling_char ) // '" set by mrun-option "-y"' -
palm/trunk/SOURCE/data_output_2d.f90
r336 r343 8 8 ! Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0) 9 9 ! Output of messages replaced by message handling routine. 10 ! Output of user defined 2D (XY) arrays at z=nzb+1 is now possible 10 11 ! 11 12 ! … … 512 513 !-- User defined quantity 513 514 CALL user_data_output_2d( av, do2d(av,if), found, grid, & 514 local_pf )515 local_pf, two_d ) 515 516 resorted = .TRUE. 516 517 … … 519 520 ELSEIF ( grid == 'zw' ) THEN 520 521 IF ( mode == 'xy' ) level_z = zw 522 ELSEIF ( grid == 'zu1' ) THEN 523 IF ( mode == 'xy' ) level_z(nzb+1) = zu(nzb+1) 521 524 ENDIF 522 525 -
palm/trunk/SOURCE/flow_statistics.f90
r305 r343 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! Added missing timeseries calculation of w"q"(0), moved timeseries q* to the 7 ! end. 6 8 ! Temperature gradient criterion for estimating the boundary layer height 7 9 ! replaced by the gradient criterion of Sullivan et al. (1998). … … 1025 1027 ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr) ! u'w' at k=0 1026 1028 ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr) ! v'w' at k=0 1027 ts_value(21,sr) = hom(nzb +12,1,pr_palm,sr) ! q*1029 ts_value(21,sr) = hom(nzb,1,48,sr) ! w"q" at k=0 1028 1030 1029 1031 IF ( ts_value(5,sr) /= 0.0 ) THEN … … 1034 1036 ENDIF 1035 1037 1038 ts_value(23,sr) = hom(nzb+12,1,pr_palm,sr) ! q* 1036 1039 ! 1037 1040 !-- Calculate additional statistics provided by the user interface -
palm/trunk/SOURCE/modules.f90
r336 r343 3 3 4 4 !------------------------------------------------------------------------------! 5 ! Actualrevisions:5 ! Current revisions: 6 6 ! ----------------- 7 ! translation error of actual -> current revisions fixed 8 ! +q* in dots_label, dots_unit. increased dots_num respectively 7 9 ! typographical error in dots_unit fixed 8 10 ! +clip_dvrp_*, cluster_size, color_interval, dvrpsize_interval, dvrp_overlap, … … 824 826 replace_num = 6 825 827 826 INTEGER :: dots_num = 2 2828 INTEGER :: dots_num = 23 827 829 828 830 CHARACTER, DIMENSION( replace_num ) :: & … … 856 858 'vmax ', 'wmax ', 'div_new', 'div_old', 'z_i_wpt', 'z_i_pt ', & 857 859 'w* ', 'w"pt"0 ', 'w"pt" ', 'wpt ', 'pt(0) ', 'pt(zp) ', & 858 'w"u"0 ', 'w"v"0 ', 'w"q"0 ', 'mo_L ', 859 ( 'unknown', i9 = 1, 7 8) /)860 'w"u"0 ', 'w"v"0 ', 'w"q"0 ', 'mo_L ', 'q* ', & 861 ( 'unknown', i9 = 1, 77) /) 860 862 861 863 CHARACTER (LEN=7), DIMENSION(dots_max) :: dots_unit = & … … 863 865 'm/s ', 'm/s ', 's-1 ', 's-1 ', 'm ', 'm ', & 864 866 'm/s ', 'K m/s ', 'K m/s ', 'K m/s ', 'K ', 'K ', & 865 'm2/s2 ', 'm2/s2 ', 'kg m/s ', 'm ', 866 ( 'unknown', i9 = 1, 7 8) /)867 'm2/s2 ', 'm2/s2 ', 'kg m/s ', 'm ', 'kg/kg ', & 868 ( 'unknown', i9 = 1, 77 ) /) 867 869 868 870 CHARACTER (LEN=7), DIMENSION(300) :: dopr_unit = 'unknown' -
palm/trunk/SOURCE/time_integration.f90
r291 r343 116 116 !-- Determine size of next time step 117 117 IF ( simulated_time /= 0.0 ) CALL timestep 118 119 118 ! 120 119 !-- Execute the user-defined actions … … 334 333 IF ( coupling_mode /= 'uncoupled' .AND. run_coupled ) THEN 335 334 time_coupling = time_coupling + dt_3d 335 336 336 ! 337 337 !-- In case of model termination initiated by the local model -
palm/trunk/SOURCE/timestep.f90
r320 r343 55 55 REAL, DIMENSION(2) :: uv_gtrans, uv_gtrans_l 56 56 REAL, DIMENSION(nzb+1:nzt) :: dxyz2_min 57 58 57 59 58 CALL cpu_log( log_point(12), 'calculate_timestep', 'start' ) … … 385 384 386 385 ENDIF 387 388 386 CALL cpu_log( log_point(12), 'calculate_timestep', 'stop' ) 389 387 390 391 388 END SUBROUTINE timestep -
palm/trunk/SOURCE/user_data_output_2d.f90
r226 r343 1 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf )1 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf, two_d ) 2 2 3 3 !------------------------------------------------------------------------------! … … 29 29 INTEGER :: av, i, j, k 30 30 31 LOGICAL :: found 31 LOGICAL :: found, two_d 32 32 33 33 REAL, DIMENSION(nxl-1:nxr+1,nys-1:nyn+1,nzb:nzt+1) :: local_pf … … 64 64 ! grid = 'zu' 65 65 66 66 67 CASE DEFAULT 67 68 found = .FALSE.
Note: See TracChangeset
for help on using the changeset viewer.
