- Timestamp:
- May 30, 2019, 4:13:37 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/StratAer/traccoag_mod.F90
r3527 r3528 235 235 ! SAI continuous emission o 236 236 dlat_loc=180./RPI/2.*(boundslat(i,1)-boundslat(i,3)) ! dlat = half difference of boundary latitudes 237 WRITE(lunout,*) "IN traccoag, dlon=",dlon238 WRITE(lunout,*) "IN traccoag, dlat=",dlat_loc239 237 IF ( xlat(i).GE.xlat_sai-dlat_loc .AND. xlat(i).LT.xlat_sai+dlat_loc .AND. & 240 238 & xlon(i).GE.xlon_sai-dlon .AND. xlon(i).LT.xlon_sai+dlon ) THEN 241 239 ! 242 WRITE(lunout,*) 'coordinates of SAI point=',xlat(i), xlon(i), day_cur, mth_cur, year_cur243 240 ! compute altLMDz 244 241 altLMDz(:)=0.0 … … 303 300 ! lat_min and lat_max 304 301 305 WRITE(lunout,*) 'IN traccoag, dlon=',dlon306 302 DO i=1,klon 307 303 ! SAI scenario with continuous emission
Note: See TracChangeset
for help on using the changeset viewer.