- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/satellite_out_spla.F90
r2630 r5082 170 170 lon_cen=(ATAN(TAN(tempo)*COS(incli))-duree_orb*tempo)*RADEG 171 171 demi_larg=demi_larg_eq/COS(rlat(i)*DTOR) 172 IF (ABS(SIN(rlat(i)*DTOR)/SIN(incli)) .GE.1.0) demi_larg=200.0173 IF (rlat(i) .GE.lat_fin.AND.rlat(i).LE.lat_debut) THEN174 IF (demi_larg .GE. 180.) THEN172 IF (ABS(SIN(rlat(i)*DTOR)/SIN(incli))>=1.0) demi_larg=200.0 173 IF (rlat(i)>=lat_fin.AND.rlat(i)<=lat_debut) THEN 174 IF (demi_larg>= 180.) THEN 175 175 masque(i)=1 176 176 ELSE … … 178 178 lon_east = MOD(lon0+lon_cen+demi_larg, 360.) 179 179 zlon = MOD(rlon(i)+360., 360.) 180 IF (lon_west .LE.lon_east) THEN181 IF (zlon .GE.lon_west.AND.zlon.LE.lon_east) masque(i)=1180 IF (lon_west<=lon_east) THEN 181 IF (zlon>=lon_west.AND.zlon<=lon_east) masque(i)=1 182 182 ELSE 183 IF (zlon .GE.lon_west.OR.zlon.LE.lon_east) masque(i)=1183 IF (zlon>=lon_west.OR.zlon<=lon_east) masque(i)=1 184 184 ENDIF 185 185 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.