Changeset 734 for LMDZ4/branches/V3_test/libf/dyn3dpar/grid_noro.F

Timestamp:

Nov 10, 2006, 3:42:49 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Nettoyage version parallele
YM/LF

File:

• LMDZ4/branches/V3_test/libf/dyn3dpar/grid_noro.F (modified) (4 diffs)

LMDZ4/branches/V3_test/libf/dyn3dpar/grid_noro.F

@@ -83 +83 @@
       REAL zsig(imar+1,jmar),zgam(imar+1,jmar),zthe(imar+1,jmar)
       REAL zpic(imar+1,jmar),zval(imar+1,jmar)
-c$$$ PB     integer mask(imar+1,jmar)
+cxxx PB     integer mask(imar+1,jmar)
       real mask(imar+1,jmar), mask_tmp(imar+1,jmar)
       real num_tot(2200,1100),num_lan(2200,1100)
@@ -274 +274 @@
          IF (weight(ii,jj) .NE. 0.0) THEN
 c  Mask
-c$$$           if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then
-c$$$             mask(ii,jj)=1
-c$$$           else
-c$$$             mask(ii,jj)=0
-c$$$           ENDIF
+cXXX           if(num_lan(ii,jj)/num_tot(ii,jj).ge.0.5)then
+cXXX             mask(ii,jj)=1
+cXXX           else
+cXXX             mask(ii,jj)=0
+cXXX           ENDIF
              if (.not. masque_lu) then
                mask(ii,jj) = num_lan(ii,jj)/num_tot(ii,jj)
@@ -319 +319 @@
        CALL MVA9(zxtzy,iim+1,jjm+1) 
        CALL MVA9(zytzy,iim+1,jjm+1)
-C$$$   Masque prenant en compte maximum de terre
-C$$$  On seuil a 10% de terre de terre car en dessous les parametres de surface n'on
-C$$$ pas de sens (PB)
+CXXX   Masque prenant en compte maximum de terre
+CXXX  On seuil a 10% de terre de terre car en dessous les parametres de surface n'on
+CXXX pas de sens (PB)
        mask_tmp= 0.0
        WHERE(mask .GE. 0.1) mask_tmp = 1.
@@ -339 +339 @@
            if(abs(xm).le.xw) xm=xw*sign(1.,xm)
 c slope: 
-c$$$           zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
-c$$$c isotropy:
-c$$$           zgam(ii,jj)=xp/xq*mask(ii,jj)
-c$$$c angle theta:
-c$$$           zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj)
-c$$$           zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)
-c$$$           zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)
-c$$$           zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)
-c$$$           zval(ii,jj)=zval(ii,jj)*mask(ii,jj)
-c$$$           zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)
-C$$* PB modif pour maque de terre fractionnaire
+cXXX           zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
+cXXXc isotropy:
+cXXX           zgam(ii,jj)=xp/xq*mask(ii,jj)
+cXXXc angle theta:
+cXXX           zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask(ii,jj)
+cXXX           zphi(ii,jj)=zmea(ii,jj)*mask(ii,jj)
+cXXX           zmea(ii,jj)=zmea(ii,jj)*mask(ii,jj)
+cXXX           zpic(ii,jj)=zpic(ii,jj)*mask(ii,jj)
+cXXX           zval(ii,jj)=zval(ii,jj)*mask(ii,jj)
+cXXX           zstd(ii,jj)=zstd(ii,jj)*mask(ii,jj)
+CXX* PB modif pour maque de terre fractionnaire
 c slope: 
            zsig(ii,jj)=sqrt(xq)*mask_tmp(ii,jj)