Changeset 1944 for LMDZ5/trunk/libf
- Timestamp:
- Jan 22, 2014, 6:39:15 PM (11 years ago)
- Location:
- LMDZ5/trunk/libf
- Files:
-
- 1 added
- 14 deleted
- 7 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d_common/grid_noro.F
r1938 r1944 83 83 REAL zsig(imar+1,jmar),zgam(imar+1,jmar),zthe(imar+1,jmar) 84 84 REAL zpic(imar+1,jmar),zval(imar+1,jmar) 85 cx $$PB integer mask(imar+1,jmar)85 cxxx PB integer mask(imar+1,jmar) 86 86 real mask(imar+1,jmar), mask_tmp(imar+1,jmar) 87 87 real num_tot(2200,1100),num_lan(2200,1100) … … 274 274 IF (weight(ii,jj) .NE. 0.0) THEN 275 275 c Mask 276 c x$$if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then277 c x$$mask(ii,jj)=1278 c x$$else279 c x$$mask(ii,jj)=0280 c x$$ENDIF276 cXXX if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then 277 cXXX mask(ii,jj)=1 278 cXXX else 279 cXXX mask(ii,jj)=0 280 cXXX ENDIF 281 281 if (.not. masque_lu) then 282 282 mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj) … … 319 319 CALL MVA9(zxtzy,iim+1,jjm+1) 320 320 CALL MVA9(zytzy,iim+1,jjm+1) 321 C x$$Masque prenant en compte maximum de terre322 C x$$On seuil a 10% de terre de terre car en dessous les parametres de surface n'on323 C x$$pas de sens (PB)321 CXXX Masque prenant en compte maximum de terre 322 CXXX On seuil a 10% de terre de terre car en dessous les parametres de surface n'on 323 CXXX pas de sens (PB) 324 324 mask_tmp= 0.0 325 325 WHERE(mask .GE. 0.1) mask_tmp = 1. … … 339 339 if(abs(xm).le.xw) xm=xw*sign(1.,xm) 340 340 c slope: 341 c x$$zsig(ii,jj)=sqrt(xq)*mask(ii,jj)342 c x$$c isotropy:343 c x$$zgam(ii,jj)=xp/xq*mask(ii,jj)344 c x$$c angle theta:345 c x$$zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj)346 c x$$zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)347 c x$$zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)348 c x$$zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)349 c x$$zval(ii,jj)=zval(ii,jj)*mask(ii,jj)350 c x$$zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)351 C x$* PB modif pour maque de terre fractionnaire341 cXXX zsig(ii,jj)=sqrt(xq)*mask(ii,jj) 342 cXXXc isotropy: 343 cXXX zgam(ii,jj)=xp/xq*mask(ii,jj) 344 cXXXc angle theta: 345 cXXX zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj) 346 cXXX zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj) 347 cXXX zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj) 348 cXXX zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj) 349 cXXX zval(ii,jj)=zval(ii,jj)*mask(ii,jj) 350 cXXX zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj) 351 CXX* PB modif pour maque de terre fractionnaire 352 352 c slope: 353 353 zsig(ii,jj)=sqrt(xq)*mask_tmp(ii,jj)
Note: See TracChangeset
for help on using the changeset viewer.