Changeset 49 for codes/icosagcm/trunk
- Timestamp:
- 07/30/12 08:13:16 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/wind.f90
r38 r49 1 1 MODULE wind_mod 2 2 3 4 5 3 CONTAINS 4 6 5 7 6 … … 19 18 ij=(j-1)*iim+i 20 19 ucenter(ij,:,l)=1/Ai(ij)* & 21 ( ne(ij,right)*ue(ij+u_right,l)*le(ij+u_right)*((xyz_v(ij+z_rdown,:)+xyz_ i(ij+z_rup,:))/2 -centroid(ij,:))&22 + ne(ij,rup)*ue(ij+u_rup,l)*le(ij+u_rup)*((xyz_v(ij+z_rup,:)+xyz_ i(ij+z_up,:))/2-centroid(ij,:)) &23 + ne(ij,lup)*ue(ij+u_lup,l)*le(ij+u_lup)*((xyz_v(ij+z_up,:)+xyz_ i(ij+z_lup,:))/2-centroid(ij,:)) &24 + ne(ij,left)*ue(ij+u_left,l)*le(ij+u_left)*((xyz_ i(ij+z_lup,:)+xyz_i(ij+z_ldown,:))/2-centroid(ij,:)) &25 + ne(ij,ldown)*ue(ij+u_ldown,l)*le(ij+u_ldown)*((xyz_ i(ij+z_ldown,:)+xyz_i(ij+z_down,:))/2-centroid(ij,:))&26 + ne(ij,rdown)*ue(ij+u_rdown,l)*le(ij+u_rdown)*((xyz_ i(ij+z_down,:)+xyz_i(ij+z_rdown,:))/2-centroid(ij,:)))20 ( ne(ij,right)*ue(ij+u_right,l)*le(ij+u_right)*((xyz_v(ij+z_rdown,:)+xyz_v(ij+z_rup,:))/2-centroid(ij,:)) & 21 + ne(ij,rup)*ue(ij+u_rup,l)*le(ij+u_rup)*((xyz_v(ij+z_rup,:)+xyz_v(ij+z_up,:))/2-centroid(ij,:)) & 22 + ne(ij,lup)*ue(ij+u_lup,l)*le(ij+u_lup)*((xyz_v(ij+z_up,:)+xyz_v(ij+z_lup,:))/2-centroid(ij,:)) & 23 + ne(ij,left)*ue(ij+u_left,l)*le(ij+u_left)*((xyz_v(ij+z_lup,:)+xyz_v(ij+z_ldown,:))/2-centroid(ij,:)) & 24 + ne(ij,ldown)*ue(ij+u_ldown,l)*le(ij+u_ldown)*((xyz_v(ij+z_ldown,:)+xyz_v(ij+z_down,:))/2-centroid(ij,:))& 25 + ne(ij,rdown)*ue(ij+u_rdown,l)*le(ij+u_rdown)*((xyz_v(ij+z_down,:)+xyz_v(ij+z_rdown,:))/2-centroid(ij,:))) 27 26 ENDDO 28 27 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.