Changeset 286 for codes/icosagcm/trunk/src/etat0_dcmip2.f90
- Timestamp:
- 10/20/14 23:42:26 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_dcmip2.f90
r186 r286 86 86 DO i=ii_begin,ii_end 87 87 ij=(j-1)*iim+i 88 CALL comp_all( xyz_i(ij,:), ps(ij),phis(ij),temp(ij,l), dummy1,dummy2)88 CALL comp_all(lon_i(ij),lat_i(ij), ps(ij),phis(ij),temp(ij,l), dummy1,dummy2) 89 89 END DO 90 90 END DO … … 100 100 DO i=ii_begin,ii_end 101 101 ij=(j-1)*iim+i 102 CALL comp_all(xyz_e(ij+u_right,:), dummy1,dummy2,dummy3, ulon(ij+u_right,l),ulat(ij+u_right,l)) 103 CALL comp_all(xyz_e(ij+u_lup,:), dummy1,dummy2,dummy3, ulon(ij+u_lup,l),ulat(ij+u_lup,l)) 104 CALL comp_all(xyz_e(ij+u_ldown,:), dummy1,dummy2,dummy3, ulon(ij+u_ldown,l),ulat(ij+u_ldown,l)) 102 CALL comp_all(lon_e(ij+u_right), lat_e(ij+u_right), & 103 dummy1,dummy2,dummy3, ulon(ij+u_right,l),ulat(ij+u_right,l)) 104 CALL comp_all(lon_e(ij+u_lup), lat_e(ij+u_lup), & 105 dummy1,dummy2,dummy3, ulon(ij+u_lup,l),ulat(ij+u_lup,l)) 106 CALL comp_all(lon_e(ij+u_ldown), lat_e(ij+u_ldown), & 107 dummy1,dummy2,dummy3, ulon(ij+u_ldown,l),ulat(ij+u_ldown,l)) 105 108 END DO 106 109 END DO … … 112 115 CONTAINS 113 116 114 SUBROUTINE comp_all( xyz, psj,phisj,tempj, ulonj,ulatj)117 SUBROUTINE comp_all(lon,lat, psj,phisj,tempj, ulonj,ulatj) 115 118 USE dcmip_initial_conditions_test_1_2_3 116 REAL(rstd), INTENT(IN) :: xyz(3)119 REAL(rstd), INTENT(IN) :: lon, lat 117 120 REAL(rstd), INTENT(OUT) :: psj,phisj,tempj,ulonj,ulatj 118 REAL :: lon,lat,dummy121 REAL :: dummy 119 122 dummy=0. 120 CALL xyz2lonlat(xyz,lon,lat)121 123 SELECT CASE (icase) 122 124 CASE(0)
Note: See TracChangeset
for help on using the changeset viewer.