Changeset 407 for palm/trunk


Ignore:
Timestamp:
Dec 1, 2009 3:01:15 PM (15 years ago)
Author:
maronga
Message:

humidity for non-flat topography implemented, re-adjustments for lcxt4

Location:
palm/trunk
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_4.1.html

    r354 r407  
    1438614386Non-flat <span style="font-weight: bold;">topography</span>
    1438714387modes may assign a</font>
    14388 kinematic sensible<font color="#000000"> <a href="chapter_4.1.html#wall_heatflux">wall_heatflux</a> at the five topography faces.</font><br>
     14388kinematic sensible<font color="#000000"> <a href="chapter_4.1.html#wall_heatflux">wall_heatflux</a> and a kinematic <a href="chapter_4.1.html#wall_humidityflux">wall_humidityflux</a> (requires <a href="chapter_4.1.html#humidity">humidity</a> = .T.) or a <a href="chapter_4.1.html#wall_scalarflux">wall_scalarflux</a> (requires <a href="chapter_4.1.html#passive_scalar">passive_scalar</a> = .T.) at the five topography faces.</font><br>
    1438914389
    1439014390      <font color="#000000">
     
    1440814408      <i>&nbsp;</i><a href="#alpha_surface">alpha_surface</a>
    1440914409= 0.0,<span style="font-style: italic;"></span>&nbsp;<a style="" href="#galilei_transformation">galilei_transformation</a>
    14410 = <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp;&nbsp;<a href="#humidity">humidity</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
     14410= <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
    1441114411
    1441214412
     
    1702817028
    1702917029
    17030 
    17031 
    17032  
    17033  
    17034  
    17035  
    17036  
     17030  <tr>
     17031
     17032
     17033
     17034
     17035
     17036
     17037 <td style="vertical-align: top;"><span style="font-weight: bold;"><a name="wall_humidityflux"></a>wall_humidityflux</span></td>
     17038
     17039
     17040
     17041
     17042
     17043
     17044
     17045      <td style="vertical-align: top;">R(5)</td>
     17046
     17047
     17048
     17049
     17050
     17051
     17052 <td style="vertical-align: top;"><span style="font-style: italic;">5 * 0.0</span></td>
     17053
     17054
     17055
     17056
     17057
     17058
     17059 <td>Prescribed
     17060kinematic humidity flux in m/s
     17061at the five topography faces:<br>
     17062
     17063
     17064
     17065
     17066
     17067
     17068 <br>
     17069
     17070
     17071
     17072
     17073
     17074
     17075 
     17076     
     17077     
     17078     
     17079     
     17080     
     17081     
     17082      <div style="margin-left: 40px;"><span style="font-weight: bold;">wall_humidityflux(0)&nbsp;&nbsp;
     17083&nbsp;</span>top face<br>
     17084
     17085
     17086
     17087
     17088
     17089
     17090 <span style="font-weight: bold;">wall_humidityflux(1)&nbsp;&nbsp;&nbsp;
     17091      </span>left face<br>
     17092
     17093
     17094
     17095
     17096
     17097
     17098 <span style="font-weight: bold;">wall_humidityflux(2)&nbsp;&nbsp;&nbsp;
     17099      </span>right face<br>
     17100
     17101
     17102
     17103
     17104
     17105
     17106 <span style="font-weight: bold;">wall_humidityflux(3)&nbsp;&nbsp;&nbsp;
     17107      </span>south face<br>
     17108
     17109
     17110
     17111
     17112
     17113
     17114 <span style="font-weight: bold;">wall_humidityflux(4)&nbsp;&nbsp;&nbsp;
     17115      </span>north face</div>
     17116
     17117
     17118
     17119
     17120
     17121
     17122 <br>
     17123
     17124
     17125
     17126
     17127
     17128
     17129
     17130This parameter applies only in case of a non-flat <a href="#topography">topography</a> and <a href="#humidity">humidity</a> = .T.
     17131</td>
     17132
     17133
     17134
     17135
     17136
     17137
     17138 </tr>
     17139 
     17140 
     17141 <tr>
     17142
     17143
     17144
     17145
     17146
     17147
     17148 <td style="vertical-align: top;"><span style="font-weight: bold;"><a name="wall_scalarflux"></a>wall_scalarflux</span></td>
     17149
     17150
     17151
     17152
     17153
     17154
     17155
     17156      <td style="vertical-align: top;">R(5)</td>
     17157
     17158
     17159
     17160
     17161
     17162
     17163 <td style="vertical-align: top;"><span style="font-style: italic;">5 * 0.0</span></td>
     17164
     17165
     17166
     17167
     17168
     17169
     17170 <td>Prescribed scalar flux in kg/(m<sup>2</sup> s)
     17171at the five topography faces:<br>
     17172
     17173
     17174
     17175
     17176
     17177
     17178 <br>
     17179
     17180
     17181
     17182
     17183
     17184
     17185 
     17186     
     17187     
     17188     
     17189     
     17190     
     17191     
     17192      <div style="margin-left: 40px;"><span style="font-weight: bold;">wall_scalarflux(0)&nbsp;&nbsp;
     17193&nbsp;</span>top face<br>
     17194
     17195
     17196
     17197
     17198
     17199
     17200 <span style="font-weight: bold;">wall_scalarflux(1)&nbsp;&nbsp;&nbsp;
     17201      </span>left face<br>
     17202
     17203
     17204
     17205
     17206
     17207
     17208 <span style="font-weight: bold;">wall_scalarflux(2)&nbsp;&nbsp;&nbsp;
     17209      </span>right face<br>
     17210
     17211
     17212
     17213
     17214
     17215
     17216 <span style="font-weight: bold;">wall_scalarflux(3)&nbsp;&nbsp;&nbsp;
     17217      </span>south face<br>
     17218
     17219
     17220
     17221
     17222
     17223
     17224 <span style="font-weight: bold;">wall_scalarflux(4)&nbsp;&nbsp;&nbsp;
     17225      </span>north face</div>
     17226
     17227
     17228
     17229
     17230
     17231
     17232 <br>
     17233
     17234
     17235
     17236
     17237
     17238
     17239
     17240This parameter applies only in case of a non-flat <a href="#topography">topography</a> and <a href="#passive_scalar">passive_scalar</a> = .T.
     17241</td>
     17242
     17243
     17244
     17245
     17246
     17247
     17248 </tr>
     17249
     17250
     17251
     17252 
    1703717253 
    1703817254 
  • palm/trunk/DOC/app/chapter_4.6.html

    r350 r407  
    1509215092
    1509315093
    15094 
     15094 <tr>
     15095
     15096
     15097
     15098 <td><a href="chapter_4.1.html#wall_humidityflux"><span style="font-weight: bold;">wall_humidityflux</span></a></td>
     15099
     15100
     15101
     15102
     15103      <td>I</td>
     15104
     15105
     15106
     15107 <td>R(5)</td>
     15108
     15109
     15110
     15111 <td><span style="font-style: italic;">5 * 0.0</span></td>
     15112
     15113
     15114
     15115 <td>Prescribed
     15116kinematic humidtiy flux in m/s&nbsp;at
     15117the five topography faces.</td>
     15118
     15119
     15120
     15121 </tr>
     15122 
     15123  <tr>
     15124
     15125
     15126
     15127 <td><a href="chapter_4.1.html#wall_scalarflux"><span style="font-weight: bold;">wall_scalarflux</span></a></td>
     15128
     15129
     15130
     15131
     15132      <td>I</td>
     15133
     15134
     15135
     15136 <td>R(5)</td>
     15137
     15138
     15139
     15140 <td><span style="font-style: italic;">5 * 0.0</span></td>
     15141
     15142
     15143
     15144 <td>Prescribed
     15145scalar flux in kg/(m<sup>2</sup> s) at
     15146the five topography faces.</td>
     15147
     15148
     15149
     15150 </tr>
     15151 
    1509515152 <tr>
    1509615153
  • palm/trunk/SCRIPTS/mbuild

    r367 r407  
    107107     # 21/04/09 - Siggi - adjustments for new IBM at DKRZ, which is now ibmh
    108108     # 25/08/09 - BjornM- adapted for lck
     109     # 01/12/09 - BjornM- re-adjustments for lcxt4
    109110
    110111
     
    11101111             then
    11111112
    1112                 print "cd ${remote_md}; echo $make_call_string > LAST_MAKE_CALL; chmod u+x LAST_MAKE_CALL; $make_call_string; [[ \$? != 0 ]] && echo MAKE_ERROR" | ssh  ${remote_username}@${remote_addres} 2>&1 | tee ${remote_host}_last_make_protokoll
     1113                print ". /opt/modules/default/init/ksh; module load Base-opts/2.1.56HDA; module load modules/3.1.6; module load pgi/9.0.4; module load PrgEnv-pgi/2.1.56HDA; cd ${remote_md}; echo $make_call_string > LAST_MAKE_CALL; chmod u+x LAST_MAKE_CALL; $make_call_string; [[ \$? != 0 ]] && echo MAKE_ERROR" | ssh  ${remote_username}@${remote_addres} 2>&1 | tee ${remote_host}_last_make_protokoll
    11131114
    11141115             else
     
    12441245                done
    12451246
     1247             elif [[ $remote_host = lcxt4 ]]
     1248             then
     1249
     1250                print ". /opt/modules/default/init/ksh; module load Base-opts/2.1.56HDA; module load modules/3.1.6; module load pgi/9.0.4; module load PrgEnv-pgi/2.1.56HDA; cd ${remote_ud}; $make_call_string; [[ \$? != 0 ]] && echo MAKE_ERROR" | ssh  ${remote_username}@${remote_addres} 2>&1 | tee ${remote_host}_last_make_protokoll
     1251
    12461252             else
    12471253
  • palm/trunk/SCRIPTS/subjob

    r400 r407  
    107107     # 03/09/09 - Siggi - PBS sgi feature directive only used if explicitly
    108108     #                    set in the config file by the user
    109      # 16/10/09 - Carolin - adjustments for archiving on SGI-ICE of binary files;
     109     # 16/10/09 - Carolin adjustments for archiving on SGI-ICE of binary files;
    110110     #                      adding special1q
     111     # 01/12/09 - BjornM- re-adjustments for lcxt4, added sno (130.75.105.113)
    111112
    112113    # VARIABLENVEREINBARUNGEN + DEFAULTWERTE
     
    203204     (blizzard1)             local_addres=136.172.40.15;  local_host=ibmh;;
    204205     (irifi)                 local_addres=130.75.105.104; local_host=lcmuk;;
     206     (sno)                   local_addres=130.75.105.113; local_host=lcmuk;;
    205207     (levanto)               local_addres=130.75.105.45;  local_host=lcmuk;;
    206208     (maestro)               local_addres=130.75.105.2;   local_host=lcmuk;;
     
    342344        (ibmy)    queue=parallel; remote_addres=165.132.26.58; submcom=/usr/lpp/LoadL/full/bin/llsubmit;;
    343345        (lcfimm)  remote_addres=172.20.4.2; submcom=/opt/torque/bin/qsub;;
    344         (lcsgib)  queue=smallq; remote_addres=130.73.232.103; submcom=/opt/moab/bin/msub;;
    345         (lcsgih)  queue=smallq; remote_addres=130.75.4.102; submcom=/opt/moab/bin/msub;;
     346        (lcsgib)  queue=smallq; remote_addres=130.73.232.102; submcom=/opt/moab/bin/msub;;
     347        (lcsgih)  queue=smallq; remote_addres=130.75.4.103; submcom=/opt/moab/bin/msub;;
    346348        (lck)     remote_addres=165.132.26.61; submcom=/usr/torque/bin/qsub;;
    347349        (lctit)   queue=lctit; remote_addres=172.17.75.161; submcom=/n1ge/TITECH_GRID/tools/bin/n1ge;;
     
    844846. /opt/modules/default/init/ksh
    845847module load modules/3.1.6
    846 module load PrgEnv-pgi/2.1.56HD
    847 module load Base-opts/2.1.56HD
    848 module load moab/5.2.5-s13441
     848module load PrgEnv-pgi/2.1.56HDA
     849module load Base-opts/2.1.56HDA
     850module load moab/5.2.5-s13570
    849851module load torque/2.3.4-snap.200809221601
     852
    850853
    851854%%END%%
     
    12741277          echo "$submcom  $job_on_remhost"
    12751278          chmod  u+x  $job_on_remhost
     1279
    12761280       elif [[ $(echo $local_host | cut -c1-5) = lcsgi ]]
    12771281       then
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r405 r407  
    11New:
    22---
     3humidity=.T. is now usable for runs with topography. wall_humidityflux and
     4wall_scalarflux are the corresponding new parin arrays.
     5(check_parameters, init_3d_model, parin)
    36
     7check parameters, init_3d_model, parin
    48
    59Changed:
     
    1923
    2024Bugfix: initial data assignments to some dvrp arrays changed due to error
    21 messages from gfortran compiler (modules=
     25messages from gfortran compiler (modules)
    2226
    2327calc_precipitation, modules, sum_up_3d_data, time_integration
  • palm/trunk/SOURCE/check_parameters.f90

    r392 r407  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Enabled passive scalar/humidity wall fluxes for non-flat topography
    77!
    88! Former revisions:
     
    331331       IF ( cloud_droplets )  THEN
    332332          WRITE( action, '(A)' )  'cloud_droplets = .TRUE.'
    333        ENDIF
    334        IF ( humidity )  THEN
    335           WRITE( action, '(A)' )  'humidity = .TRUE.'
    336333       ENDIF
    337334       IF ( .NOT. prandtl_layer )  THEN
     
    858855          q_vertical_gradient_level = s_vertical_gradient_level
    859856          surface_waterflux         = surface_scalarflux
     857          wall_humidityflux         = wall_scalarflux
    860858       ENDIF
    861859
  • palm/trunk/SOURCE/init_3d_model.f90

    r394 r407  
    88! -----------------
    99! var_ts is replaced by dots_max
    10 !
     10! Enabled passive scalar/humidity wall fluxes for non-flat topography
    1111!
    1212! Former revisions:
     
    735735             IF ( constant_waterflux )  THEN
    736736                qsws   = surface_waterflux
     737!
     738!--             Over topography surface_waterflux is replaced by
     739!--             wall_humidityflux(0)
     740                IF ( TRIM( topography ) /= 'flat' )  THEN
     741                   wall_qflux = wall_humidityflux
     742                   DO  i = nxl-1, nxr+1
     743                      DO  j = nys-1, nyn+1
     744                         IF ( nzb_s_inner(j,i) /= 0 )  THEN
     745                            qsws(j,i) = wall_qflux(0)
     746                         ENDIF
     747                      ENDDO
     748                   ENDDO
     749                ENDIF
    737750                IF ( ASSOCIATED( qsws_m ) )  qsws_m = qsws
    738751             ENDIF
  • palm/trunk/SOURCE/parin.f90

    r392 r407  
    44! Current revisions:
    55! -----------------
    6 !
     6! +wall_humidityflux, wall_scalarflux
    77!
    88! Former revisions:
     
    131131             use_surface_fluxes, use_top_fluxes, use_ug_for_galilei_tr, &
    132132             use_upstream_for_tke, vg_surface, vg_vertical_gradient, &
    133              vg_vertical_gradient_level, v_bulk, wall_adjustment, wall_heatflux
     133             vg_vertical_gradient_level, v_bulk, wall_adjustment, &
     134             wall_heatflux, wall_humidityflux, wall_scalarflux
    134135
    135136
  • palm/trunk/SOURCE/prognostic_equations.f90

    r392 r407  
    405405       ENDIF
    406406    ENDIF
    407          
     407
    408408!
    409409!-- pt-tendency terms with no communication
Note: See TracChangeset for help on using the changeset viewer.