Changeset 521 for codes/icosagcm/trunk/src/caldyn_kernels.f90
- Timestamp:
- 01/17/17 16:27:34 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/caldyn_kernels.f90
r519 r521 60 60 !DIR$ SIMD 61 61 DO ij=ij_begin_ext,ij_end_ext 62 IF(DEC) THEN ! ps is actually Ms 63 m = mass_dak(l)+(ps(ij)*g+ptop)*mass_dbk(l) 64 ELSE 65 m = mass_dak(l)+ps(ij)*mass_dbk(l) 66 END IF 62 m = mass_dak(l)+(ps(ij)*g+ptop)*mass_dbk(l) ! ps is actually Ms 67 63 rhodz(ij,l) = m/g 68 64 theta(ij,l) = theta_rhodz(ij,l)/rhodz(ij,l) … … 85 81 !DIR$ SIMD 86 82 DO ij=ij_begin_ext,ij_end_ext 87 IF(DEC) THEN 88 etav= 1./Av(ij+z_up)*( ne_rup * u(ij+u_rup,l) & 89 + ne_left * u(ij+t_rup+u_left,l) & 90 - ne_lup * u(ij+u_lup,l) ) 91 92 hv = Riv2(ij,vup) * rhodz(ij,l) & 93 + Riv2(ij+t_rup,vldown) * rhodz(ij+t_rup,l) & 94 + Riv2(ij+t_lup,vrdown) * rhodz(ij+t_lup,l) 95 qv(ij+z_up,l) = ( etav+fv(ij+z_up) )/hv 96 97 etav = 1./Av(ij+z_down)*( ne_ldown * u(ij+u_ldown,l) & 98 + ne_right * u(ij+t_ldown+u_right,l) & 99 - ne_rdown * u(ij+u_rdown,l) ) 100 hv = Riv2(ij,vdown) * rhodz(ij,l) & 101 + Riv2(ij+t_ldown,vrup) * rhodz(ij+t_ldown,l) & 102 + Riv2(ij+t_rdown,vlup) * rhodz(ij+t_rdown,l) 103 qv(ij+z_down,l) =( etav+fv(ij+z_down) )/hv 104 ELSE 105 etav= 1./Av(ij+z_up)*( ne_rup * u(ij+u_rup,l) * de(ij+u_rup) & 106 + ne_left * u(ij+t_rup+u_left,l) * de(ij+t_rup+u_left) & 107 - ne_lup * u(ij+u_lup,l) * de(ij+u_lup) ) 108 109 hv = Riv2(ij,vup) * rhodz(ij,l) & 110 + Riv2(ij+t_rup,vldown) * rhodz(ij+t_rup,l) & 111 + Riv2(ij+t_lup,vrdown) * rhodz(ij+t_lup,l) 112 qv(ij+z_up,l) = ( etav+fv(ij+z_up) )/hv 113 114 etav = 1./Av(ij+z_down)*( ne_ldown * u(ij+u_ldown,l) * de(ij+u_ldown) & 115 + ne_right * u(ij+t_ldown+u_right,l) * de(ij+t_ldown+u_right) & 116 - ne_rdown * u(ij+u_rdown,l) * de(ij+u_rdown) ) 117 hv = Riv2(ij,vdown) * rhodz(ij,l) & 118 + Riv2(ij+t_ldown,vrup) * rhodz(ij+t_ldown,l) & 119 + Riv2(ij+t_rdown,vlup) * rhodz(ij+t_rdown,l) 120 qv(ij+z_down,l) =( etav+fv(ij+z_down) )/hv 121 END IF 83 etav= 1./Av(ij+z_up)*( ne_rup * u(ij+u_rup,l) & 84 + ne_left * u(ij+t_rup+u_left,l) & 85 - ne_lup * u(ij+u_lup,l) ) 86 hv = Riv2(ij,vup) * rhodz(ij,l) & 87 + Riv2(ij+t_rup,vldown) * rhodz(ij+t_rup,l) & 88 + Riv2(ij+t_lup,vrdown) * rhodz(ij+t_lup,l) 89 qv(ij+z_up,l) = ( etav+fv(ij+z_up) )/hv 90 91 etav = 1./Av(ij+z_down)*( ne_ldown * u(ij+u_ldown,l) & 92 + ne_right * u(ij+t_ldown+u_right,l) & 93 - ne_rdown * u(ij+u_rdown,l) ) 94 hv = Riv2(ij,vdown) * rhodz(ij,l) & 95 + Riv2(ij+t_ldown,vrup) * rhodz(ij+t_ldown,l) & 96 + Riv2(ij+t_rdown,vlup) * rhodz(ij+t_rdown,l) 97 qv(ij+z_down,l) =( etav+fv(ij+z_down) )/hv 122 98 ENDDO 123 99
Note: See TracChangeset
for help on using the changeset viewer.