Ignore:
Timestamp:
Jan 22, 2014, 6:39:15 PM (10 years ago)
Author:
lguez
Message:

Created directory dyn3d_common. Began moving to dyn3d_common files
that were duplicated in dyn3d, dyn3dpar and dyn3dmem. Adapted
compilation files for makegcm, makelmdz and makelmdz_fcm. Adapted
install_1d_src.sh.

Location:
LMDZ5/trunk/libf/dyn3d_common
Files:
1 added
7 moved

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3d_common/grid_noro.F

    r1938 r1944  
    8383      REAL zsig(imar+1,jmar),zgam(imar+1,jmar),zthe(imar+1,jmar)
    8484      REAL zpic(imar+1,jmar),zval(imar+1,jmar)
    85 cx$$ PB     integer mask(imar+1,jmar)
     85cxxx PB     integer mask(imar+1,jmar)
    8686      real mask(imar+1,jmar), mask_tmp(imar+1,jmar)
    8787      real num_tot(2200,1100),num_lan(2200,1100)
     
    274274         IF (weight(ii,jj) .NE. 0.0) THEN
    275275c  Mask
    276 cx$$           if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then
    277 cx$$             mask(ii,jj)=1
    278 cx$$           else
    279 cx$$             mask(ii,jj)=0
    280 cx$$           ENDIF
     276cXXX           if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then
     277cXXX             mask(ii,jj)=1
     278cXXX           else
     279cXXX             mask(ii,jj)=0
     280cXXX           ENDIF
    281281             if (.not. masque_lu) then
    282282               mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj)
     
    319319       CALL MVA9(zxtzy,iim+1,jjm+1)
    320320       CALL MVA9(zytzy,iim+1,jjm+1)
    321 Cx$$   Masque prenant en compte maximum de terre
    322 Cx$$  On seuil a 10% de terre de terre car en dessous les parametres de surface n'on
    323 Cx$$ pas de sens (PB)
     321CXXX   Masque prenant en compte maximum de terre
     322CXXX  On seuil a 10% de terre de terre car en dessous les parametres de surface n'on
     323CXXX pas de sens (PB)
    324324       mask_tmp= 0.0
    325325       WHERE(mask .GE. 0.1) mask_tmp = 1.
     
    339339           if(abs(xm).le.xw) xm=xw*sign(1.,xm)
    340340c slope:
    341 cx$$           zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
    342 cx$$c isotropy:
    343 cx$$           zgam(ii,jj)=xp/xq*mask(ii,jj)
    344 cx$$c angle theta:
    345 cx$$           zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj)
    346 cx$$           zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)
    347 cx$$           zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)
    348 cx$$           zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)
    349 cx$$           zval(ii,jj)=zval(ii,jj)*mask(ii,jj)
    350 cx$$           zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)
    351 Cx$* PB modif pour maque de terre fractionnaire
     341cXXX           zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
     342cXXXc isotropy:
     343cXXX           zgam(ii,jj)=xp/xq*mask(ii,jj)
     344cXXXc angle theta:
     345cXXX           zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj)
     346cXXX           zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)
     347cXXX           zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)
     348cXXX           zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)
     349cXXX           zval(ii,jj)=zval(ii,jj)*mask(ii,jj)
     350cXXX           zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)
     351CXX* PB modif pour maque de terre fractionnaire
    352352c slope:
    353353           zsig(ii,jj)=sqrt(xq)*mask_tmp(ii,jj)
Note: See TracChangeset for help on using the changeset viewer.