Changeset 388 for palm/trunk/SOURCE
- Timestamp:
- Sep 23, 2009 9:40:33 AM (15 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 10 edited
-
CURRENT_MODIFICATIONS (modified) (2 diffs)
-
check_parameters.f90 (modified) (7 diffs)
-
eqn_state_seawater.f90 (modified) (6 diffs)
-
flow_statistics.f90 (modified) (7 diffs)
-
init_3d_model.f90 (modified) (3 diffs)
-
init_ocean.f90 (modified) (6 diffs)
-
modules.f90 (modified) (2 diffs)
-
production_e.f90 (modified) (4 diffs)
-
prognostic_equations.f90 (modified) (9 diffs)
-
time_integration.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r377 r388 125 125 Errors: 126 126 ------ 127 Bugfix: Initial hydrostatic pressure profile in case of ocean runs is now 128 calculated in 5 iteration steps. (init_ocean) 129 130 Bugfix: wrong sign in buoyancy production of ocean part in case of not using 131 the reference density (only in 3D routine production_e) (production_e) 132 127 133 Bugfix: output of averaged 2d/3d quantities requires that an avaraging 128 134 interval has been set, respective error message is included (check_parameters) … … 183 189 Bugfix: initial setting of time_coupling in coupled restart runs (time_integration) 184 190 185 advec_particles, check_parameters, cpu_log, data_output_2d, data_output_3d, header, init_3d_model, init_particles, modules, netcdf, prandtl_fluxes, production_e, read_var_list, time_integration, user_last_actions, write_var_list191 advec_particles, check_parameters, cpu_log, data_output_2d, data_output_3d, header, init_3d_model, init_particles, init_ocean, modules, netcdf, prandtl_fluxes, production_e, read_var_list, time_integration, user_last_actions, write_var_list -
palm/trunk/SOURCE/check_parameters.f90
r376 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Check profiles fpr prho and hyp. 6 7 ! Bugfix: output of averaged 2d/3d quantities requires that an avaraging 7 8 ! interval has been set, respective error message is included … … 933 934 934 935 ! 935 !-- If required compute the profile of leaf area density used in the plant canopy model 936 !-- If required compute the profile of leaf area density used in the plant 937 !-- canopy model 936 938 IF ( plant_canopy ) THEN 937 939 … … 967 969 968 970 ! 969 !-- In case of no given leaf area density gradients, choose a vanishing gradient 971 !-- In case of no given leaf area density gradients, choose a vanishing 972 !-- gradient 970 973 IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN 971 974 lad_vertical_gradient_level(1) = 0.0 … … 2089 2092 2090 2093 CASE ( 'rho' ) 2091 dopr_index(i) = 64 2092 dopr_unit(i) = 'kg/m3' 2093 hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 ) 2094 IF ( .NOT. ocean ) THEN 2095 message_string = 'data_output_pr = ' // & 2096 TRIM( data_output_pr(i) ) // ' is not imp' // & 2097 'lemented for ocean = .FALSE.' 2098 CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 ) 2099 ELSE 2100 dopr_index(i) = 64 2101 dopr_unit(i) = 'kg/m3' 2102 hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 ) 2103 ENDIF 2094 2104 2095 2105 CASE ( 'w"sa"' ) … … 2150 2160 hom(:,2,70,:) = SPREAD( zu, 2, statistic_regions+1 ) 2151 2161 ENDIF 2162 2163 CASE ( 'prho' ) 2164 IF ( .NOT. ocean ) THEN 2165 message_string = 'data_output_pr = ' // & 2166 TRIM( data_output_pr(i) ) // ' is not imp' // & 2167 'lemented for ocean = .FALSE.' 2168 CALL message( 'check_parameters', 'PA0091', 1, 2, 0, 6, 0 ) 2169 ELSE 2170 dopr_index(i) = 71 2171 dopr_unit(i) = 'kg/m3' 2172 hom(:,2,71,:) = SPREAD( zu, 2, statistic_regions+1 ) 2173 ENDIF 2174 2175 CASE ( 'hyp' ) 2176 dopr_index(i) = 72 2177 dopr_unit(i) = 'kPa' 2178 hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 ) 2152 2179 2153 2180 CASE DEFAULT … … 2929 2956 !-- Check pressure gradient conditions 2930 2957 IF ( dp_external .AND. conserve_volume_flow ) THEN 2931 WRITE( message_string, * ) 'Both dp_external and conserve_volume_flo w', &2932 ' are .TRUE. but one of them must be .FALSE.'2958 WRITE( message_string, * ) 'Both dp_external and conserve_volume_flo', & 2959 'w are .TRUE. but one of them must be .FALSE.' 2933 2960 CALL message( 'check_parameters', 'PA0150', 1, 2, 0, 6, 0 ) 2934 2961 ENDIF … … 2940 2967 ENDIF 2941 2968 IF ( .NOT. ANY( dpdxy /= 0.0 ) ) THEN 2942 WRITE( message_string, * ) 'dp_external is .TRUE. but dpdxy is ze ro',&2943 ' , i.e. the external pressure gradient & will not be applied'2969 WRITE( message_string, * ) 'dp_external is .TRUE. but dpdxy is ze', & 2970 'ro, i.e. the external pressure gradient & will not be applied' 2944 2971 CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 ) 2945 2972 ENDIF -
palm/trunk/SOURCE/eqn_state_seawater.f90
r336 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! First constant in array den also defined as type double 6 ! Potential density is additionally calculated in eqn_state_seawater, 7 ! first constant in array den also defined as type double 7 8 ! 8 9 ! Former revisions: … … 69 70 INTEGER :: i, j, k 70 71 71 REAL :: p 1, p2, p3, pt1, pt2, pt3, pt4, sa1, sa15, sa272 REAL :: pden, pnom, p1, p2, p3, pt1, pt2, pt3, pt4, sa1, sa15, sa2 72 73 73 74 DO i = nxl, nxr … … 76 77 ! 77 78 !-- Pressure is needed in dbar 78 ! p1 = hyp(0) * 1E-479 ! p1 = 0.080 79 p1 = hyp(k) * 1E-4 81 80 p2 = p1 * p1 … … 93 92 sa2 = sa1 * sa1 94 93 95 rho(k,j,i) = & 96 ( nom(1) + nom(2)*pt1 + nom(3)*pt2 + & 97 nom(4)*pt3 + nom(5)*sa1 + nom(6)*sa1*pt1 + & 98 nom(7)*sa2 + nom(8)*p1 + nom(9)*p1*pt2 + & 99 nom(10)*p1*sa1 + nom(11)*p2 + nom(12)*p2*pt2 & 100 ) / & 101 ( den(1) + den(2)*pt1 + den(3)*pt2 + & 102 den(4)*pt3 + den(5)*pt4 + den(6)*sa1 + & 103 den(7)*sa1*pt1 + den(8)*sa1*pt3 + den(9)*sa15 + & 104 den(10)*sa15*pt2 + den(11)*p1 + den(12)*p2*pt3 + & 105 den(13)*p3*pt1 & 106 ) 94 pnom = nom(1) + nom(2)*pt1 + nom(3)*pt2 + & 95 nom(4)*pt3 + nom(5)*sa1 + nom(6)*sa1*pt1 + & 96 nom(7)*sa2 97 98 pden = den(1) + den(2)*pt1 + den(3)*pt2 + & 99 den(4)*pt3 + den(5)*pt4 + den(6)*sa1 + & 100 den(7)*sa1*pt1 + den(8)*sa1*pt3 + den(9)*sa15 + & 101 den(10)*sa15*pt2 102 103 ! 104 !-- Potential density (without pressure terms) 105 prho(k,j,i) = pnom / pden 106 107 pnom = pnom + nom(8)*p1 + nom(9)*p1*pt2 + & 108 nom(10)*p1*sa1 + nom(11)*p2 + nom(12)*p2*pt2 109 110 pden = pden + den(11)*p1 + den(12)*p2*pt3 + & 111 den(13)*p3*pt1 112 113 ! 114 !-- In-situ density 115 rho(k,j,i) = pnom / pden 107 116 108 117 ENDDO 109 118 ! 110 119 !-- Neumann conditions are assumed at bottom and top boundary 111 rho(nzt+1,j,i) = rho(nzt,j,i) 112 rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i) 120 prho(nzt+1,j,i) = prho(nzt,j,i) 121 prho(nzb_s_inner(j,i),j,i) = prho(nzb_s_inner(j,i)+1,j,i) 122 rho(nzt+1,j,i) = rho(nzt,j,i) 123 rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i) 124 113 125 ENDDO 114 126 ENDDO … … 129 141 INTEGER :: i, j, k 130 142 131 REAL :: p 1, p2, p3, pt1, pt2, pt3, pt4, sa1, sa15, sa2143 REAL :: pden, pnom, p1, p2, p3, pt1, pt2, pt3, pt4, sa1, sa15, sa2 132 144 133 145 DO k = nzb_s_inner(j,i)+1, nzt 134 146 ! 135 147 !-- Pressure is needed in dbar 136 ! p1 = hyp(0) * 1E-4137 ! p1 = 0.0138 148 p1 = hyp(k) * 1E-4 139 149 p2 = p1 * p1 … … 151 161 sa2 = sa1 * sa1 152 162 153 rho(k,j,i) = ( nom(1) + nom(2)*pt1 + nom(3)*pt2 + & 154 nom(4)*pt3 + nom(5)*sa1 + nom(6)*sa1*pt1 + & 155 nom(7)*sa2 + nom(8)*p1 + nom(9)*p1*pt2 + & 156 nom(10)*p1*sa1 + nom(11)*p2 + nom(12)*p2*pt2 & 157 ) / & 158 ( den(1) + den(2)*pt1 + den(3)*pt2 + & 159 den(4)*pt3 + den(5)*pt4 + den(6)*sa1 + & 160 den(7)*sa1*pt1 + den(8)*sa1*pt3 + den(9)*sa15 + & 161 den(10)*sa15*pt2 + den(11)*p1 + den(12)*p2*pt3 + & 162 den(13)*p3*pt1 & 163 ) 163 pnom = nom(1) + nom(2)*pt1 + nom(3)*pt2 + & 164 nom(4)*pt3 + nom(5)*sa1 + nom(6)*sa1*pt1 + & 165 nom(7)*sa2 166 167 pden = den(1) + den(2)*pt1 + den(3)*pt2 + & 168 den(4)*pt3 + den(5)*pt4 + den(6)*sa1 + & 169 den(7)*sa1*pt1 + den(8)*sa1*pt3 + den(9)*sa15 + & 170 den(10)*sa15*pt2 171 172 ! 173 !-- Potential density (without pressure terms) 174 prho(k,j,i) = pnom / pden 175 176 pnom = pnom + nom(8)*p1 + nom(9)*p1*pt2 + & 177 nom(10)*p1*sa1 + nom(11)*p2 + nom(12)*p2*pt2 178 pden = pden + den(11)*p1 + den(12)*p2*pt3 + & 179 den(13)*p3*pt1 180 181 ! 182 !-- In-situ density 183 rho(k,j,i) = pnom / pden 184 185 164 186 ENDDO 187 165 188 ! 166 189 !-- Neumann conditions are assumed at bottom and top boundary 167 rho(nzt+1,j,i) = rho(nzt,j,i) 168 rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i) 190 prho(nzt+1,j,i) = prho(nzt,j,i) 191 prho(nzb_s_inner(j,i),j,i) = prho(nzb_s_inner(j,i)+1,j,i) 192 rho(nzt+1,j,i) = rho(nzt,j,i) 193 rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i) 169 194 170 195 END SUBROUTINE eqn_state_seawater_ij -
palm/trunk/SOURCE/flow_statistics.f90
r343 r388 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! Vertical profiles of potential density and hydrostatic pressure are 7 ! calculated. 6 8 ! Added missing timeseries calculation of w"q"(0), moved timeseries q* to the 7 9 ! end. … … 531 533 IF ( humidity ) THEN 532 534 sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + & 533 qswst(j,i) * rmask(j,i,sr) 535 qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv") 534 536 IF ( cloud_physics ) THEN 535 537 sums_l(nzt,45,tn) = sums_l(nzt,45,tn) + ( & … … 546 548 IF ( passive_scalar ) THEN 547 549 sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + & 548 qswst(j,i) * rmask(j,i,sr) 550 qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv") 549 551 ENDIF 550 552 ENDIF … … 597 599 rmask(j,i,sr) 598 600 sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) * & 601 rmask(j,i,sr) 602 sums_l(k,71,tn) = sums_l(k,71,tn) + prho(k,j,i) * & 599 603 rmask(j,i,sr) 600 604 ENDIF … … 644 648 ! 645 649 !-- Density at top follows Neumann condition 646 IF ( ocean ) sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn) 650 IF ( ocean ) THEN 651 sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn) 652 sums_l(nzt+1,71,tn) = sums_l(nzt,71,tn) 653 ENDIF 647 654 648 655 ! … … 897 904 hom(:,1,69,sr) = sums(:,69) ! w"e + w"p"/rho 898 905 hom(:,1,70,sr) = sums(:,70) ! q*2 906 hom(:,1,71,sr) = sums(:,71) ! prho 907 hom(:,1,72,sr) = hyp * 1E-4 ! hyp in kPa 899 908 900 909 hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1) … … 916 925 !-- is less than 1.5 times the height where the heat flux becomes negative 917 926 !-- (positive) for the first time. 918 !-- NOTE: This criterion is still capable of improving!919 927 z_i(1) = 0.0 920 928 first = .TRUE. -
palm/trunk/SOURCE/init_3d_model.f90
r359 r388 7 7 ! Current revisions: 8 8 ! ----------------- 9 ! Initialization of prho added. 9 10 ! bugfix: correction of initial volume flow for non-flat topography 10 11 ! bugfix: zero initialization of arrays within buildings for 'cyclic_fill' … … 234 235 ALLOCATE( saswsb_1(nys-1:nyn+1,nxl-1:nxr+1), & 235 236 saswst_1(nys-1:nyn+1,nxl-1:nxr+1) ) 236 ALLOCATE( rho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 237 sa_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 238 sa_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 237 ALLOCATE( prho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 238 rho_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 239 sa_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 240 sa_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), & 239 241 sa_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) 240 rho => rho_1 ! routine calc_mean_profile requires density to be a 241 ! pointer 242 prho => prho_1 243 rho => rho_1 ! routines calc_mean_profile and diffusion_e require 244 ! density to be apointer 242 245 IF ( humidity_remote ) THEN 243 246 ALLOCATE( qswst_remote(nys-1:nyn+1,nxl-1:nxr+1) ) … … 1328 1331 !-- Initialize quantities needed for the ocean model 1329 1332 CALL init_ocean 1333 1330 1334 ELSE 1331 1335 ! -
palm/trunk/SOURCE/init_ocean.f90
r366 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Bugfix: First calculation of hyp(0) changed 6 ! Bugfix: Initial profiles of hydrostatic pressure and density are calculated 7 ! iteratively. First calculation of hyp(0) changed. 7 8 ! 8 9 ! Former revisions: … … 32 33 INTEGER :: k, n 33 34 34 REAL :: sa_l, pt_l , rho_l35 REAL :: sa_l, pt_l 35 36 36 37 REAL, DIMENSION(nzb:nzt+1) :: rho_init … … 45 46 !-- Calculate initial vertical profile of hydrostatic pressure (in Pa) 46 47 !-- and the reference density (used later in buoyancy term) 48 !-- First step: Calculate pressure using reference density 47 49 hyp(nzt+1) = surface_pressure * 100.0 48 50 … … 55 57 hyp(0) = hyp(1) + rho_surface * g * dzu(1) 56 58 57 IF ( myid == 0 ) THEN 58 print*,'hydro pres using rho_surface' 59 DO k = nzt+1, 0, -1 60 print*, 'k = ', k, ' hyp = ', hyp(k) 61 ENDDO 62 print*, ' ' 63 ENDIF 64 59 ! 60 !-- Second step: Iteratively calculate in situ density (based on presssure) 61 !-- and pressure (based on in situ density) 65 62 DO n = 1, 5 66 63 … … 85 82 ENDDO 86 83 87 IF ( myid == 0 ) THEN88 print*,'hydro pres / rho n = ', n89 DO k = nzt+1, 0, -190 print*, 'k = ', k, ' hyp = ', hyp(k), ' rho = ', rho_init(k)91 ENDDO92 print*, ' '93 ENDIF94 95 84 ENDDO 96 85 … … 111 100 112 101 ! 113 !-- Calculate the initial potential density, based on the initial114 !-- temperature and salinity profile102 !-- Calculate the 3d array of initial in situ and potential density, 103 !-- based on the initial temperature and salinity profile 115 104 CALL eqn_state_seawater 116 105 -
palm/trunk/SOURCE/modules.f90
r367 r388 5 5 ! Current revisions: 6 6 ! ----------------- 7 ! +prho, prho_1 7 8 ! +bc_lr_cyc, bc_ns_cyc 8 9 ! +output_for_t0 … … 169 170 170 171 REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET :: & 171 e_1, e_2, e_3, kh_1, kh_2, km_1, km_2, p, p t_1, pt_2, pt_3, q_1,&172 q_ 2, q_3, ql_1, ql_2, rho_1, sa_1, sa_2, sa_3, u_1, u_2, u_3, v_1, &173 v_ 2, v_3, vpt_1, vpt_2, w_1, w_2, w_3172 e_1, e_2, e_3, kh_1, kh_2, km_1, km_2, p, prho_1, pt_1, pt_2, pt_3, & 173 q_1, q_2, q_3, ql_1, ql_2, rho_1, sa_1, sa_2, sa_3, u_1, u_2, u_3, & 174 v_1, v_2, v_3, vpt_1, vpt_2, w_1, w_2, w_3 174 175 175 176 REAL, DIMENSION(:,:,:), POINTER :: & 176 e, e_m, e_p, kh, kh_m, km, km_m, p t, pt_m, pt_p, q, q_m, q_p, ql,&177 ql _c, rho, sa, sa_p, te_m, tpt_m, tq_m, tsa_m, tu_m, tv_m, tw_m, u,&178 u _m, u_p, v, v_m, v_p, vpt, vpt_m, w, w_m, w_p177 e, e_m, e_p, kh, kh_m, km, km_m, prho, pt, pt_m, pt_p, q, q_m, q_p, & 178 ql, ql_c, rho, sa, sa_p, te_m, tpt_m, tq_m, tsa_m, tu_m, tv_m, tw_m, & 179 u, u_m, u_p, v, v_m, v_p, vpt, vpt_m, w, w_m, w_p 179 180 180 181 REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: rif_wall -
palm/trunk/SOURCE/production_e.f90
r364 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Bugfix: wrong sign in buoyancy production of ocean part in case of not using 7 ! the reference density (only in 3D routine production_e) 6 8 ! Bugfix to avoid zero division by km_neutral 7 9 ! … … 450 452 DO j = nys, nyn 451 453 DO k = nzb_s_inner(j,i)+1, nzt 452 tend(k,j,i) = tend(k,j,i) + 453 kh(k,j,i) * g / prho_reference * &454 tend(k,j,i) = tend(k,j,i) + & 455 kh(k,j,i) * g / rho_reference * & 454 456 ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k) 455 457 ENDDO … … 487 489 DO j = nys, nyn 488 490 DO k = nzb_s_inner(j,i)+1, nzt 489 tend(k,j,i) = tend(k,j,i) -&491 tend(k,j,i) = tend(k,j,i) + & 490 492 kh(k,j,i) * g / rho(k,j,i) * & 491 493 ( rho(k+1,j,i)-rho(k-1,j,i) ) * dd2zu(k) … … 946 948 !-- bottom and top surface layer 947 949 DO k = nzb_s_inner(j,i)+1, nzt 948 tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &950 tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / rho_reference * & 949 951 ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k) 950 952 ENDDO -
palm/trunk/SOURCE/prognostic_equations.f90
r240 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! prho is used instead of rho in diffusion_e, 6 7 ! external pressure gradient 7 8 ! … … 330 331 CALL coriolis( i, j, 3 ) 331 332 IF ( ocean ) THEN 332 CALL buoyancy( i, j, rho, prho_reference, 3, 64 )333 CALL buoyancy( i, j, rho, rho_reference, 3, 64 ) 333 334 ELSE 334 335 IF ( .NOT. humidity ) THEN … … 725 726 IF ( .NOT. humidity ) THEN 726 727 IF ( ocean ) THEN 727 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, 728 l_grid, rho, prho_reference, rif, tend,&729 zu, zw )728 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, & 729 l_grid, prho, prho_reference, rif, & 730 tend, zu, zw ) 730 731 ELSE 731 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, &732 l_grid, pt, pt_reference, rif, tend, &732 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, & 733 l_grid, pt, pt_reference, rif, tend, & 733 734 zu, zw ) 734 735 ENDIF … … 769 770 IF ( ocean ) THEN 770 771 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, & 771 km, l_grid, rho, prho_reference,&772 km, l_grid, prho, prho_reference, & 772 773 rif, tend, zu, zw ) 773 774 ELSE … … 1013 1014 CALL coriolis( i, j, 3 ) 1014 1015 IF ( ocean ) THEN 1015 CALL buoyancy( i, j, rho, prho_reference, 3, 64 )1016 CALL buoyancy( i, j, rho, rho_reference, 3, 64 ) 1016 1017 ELSE 1017 1018 IF ( .NOT. humidity ) THEN … … 1263 1264 IF ( .NOT. humidity ) THEN 1264 1265 IF ( ocean ) THEN 1265 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, &1266 km, l_grid, rho, prho_reference, &1266 CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, & 1267 km, l_grid, prho, prho_reference, & 1267 1268 rif, tend, zu, zw ) 1268 1269 ELSE … … 1552 1553 CALL coriolis( 3 ) 1553 1554 IF ( ocean ) THEN 1554 CALL buoyancy( rho, prho_reference, 3, 64 )1555 CALL buoyancy( rho, rho_reference, 3, 64 ) 1555 1556 ELSE 1556 1557 IF ( .NOT. humidity ) THEN … … 1967 1968 IF ( .NOT. humidity ) THEN 1968 1969 IF ( ocean ) THEN 1969 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, rho,&1970 prho _reference, rif, tend, zu, zw )1970 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, & 1971 prho, prho_reference, rif, tend, zu, zw ) 1971 1972 ELSE 1972 1973 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, pt, & … … 2004 2005 IF ( ocean ) THEN 2005 2006 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, & 2006 rho, prho_reference, rif, tend, zu, zw )2007 prho, prho_reference, rif, tend, zu, zw ) 2007 2008 ELSE 2008 2009 CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, & -
palm/trunk/SOURCE/time_integration.f90
r348 r388 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Using prho instead of rho in diffusvities. 6 7 ! Coupling with independent precursor runs. 7 8 ! Bugfix: output of particle time series only if particle advection is switched … … 284 285 IF ( .NOT. humidity ) THEN 285 286 IF ( ocean ) THEN 286 CALL diffusivities( rho, prho_reference )287 CALL diffusivities( prho, prho_reference ) 287 288 ELSE 288 289 CALL diffusivities( pt, pt_reference )
Note: See TracChangeset
for help on using the changeset viewer.
