Changeset 3825 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/init_be.F90
- Timestamp:
- May 6, 2015, 12:14:12 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/init_be.F90
r3814 r3825 5 5 6 6 USE dimphy 7 USE comgeomphy8 7 USE infotrac_phy, ONLY : nbtr 9 8 USE indice_sol_mod 9 USE geometry_mod, ONLY : longitude, latitude 10 10 11 11 IMPLICIT NONE … … 114 114 115 115 DO i = 1,klon 116 qcos=sin(glt)*sin( rlatd(i))116 qcos=sin(glt)*sin(latitude(i)) 117 117 !!jyg 118 118 !! qcos=qcos+cos(glt)*cos(rlatd(i))*cos(rlond(i)+glg) 119 qcos=qcos+cos(glt)*cos( rlatd(i))*cos(rlond(i)-glg)119 qcos=qcos+cos(glt)*cos(latitude(i))*cos(longitude(i)-glg) 120 120 !!jyg end 121 121 IF ( qcos .LT. -1.) qcos = -1.
Note: See TracChangeset
for help on using the changeset viewer.