Changeset 228 for LMDZ.3.3/branches


Ignore:
Timestamp:
Jun 20, 2001, 12:33:59 PM (23 years ago)
Author:
lmdzadmin
Message:

detail dans l'initialisation
LF

Location:
LMDZ.3.3/branches/rel-LF/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/gcm.F

    r200 r228  
    213213
    214214c     OPEN( 99,file ='run.def',status='old',form='formatted')
    215       CALL defrun_new( 99, .TRUE. , clesphy0 )
     215c      CALL defrun_new( 99, .TRUE. , clesphy0 )
     216      CALL conf_gcm( 99, .TRUE. , clesphy0 )
    216217c     CLOSE(99)
    217218
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/grid_noro.F

    r177 r228  
    5252      IMPLICIT REAL(X,Z)
    5353     
    54           parameter(iusn=2160,jusn=1080,iext=216)
     54          parameter(iusn=2160,jusn=1080,iext=216, epsfra = 1.e-5)
    5555#include "dimensions.h"
    5656          REAL xusn(iusn+2*iext),yusn(jusn+2)   
     
    7878      REAL zpic(imar+1,jmar),zval(imar+1,jmar)
    7979c$$$ PB     integer mask(imar+1,jmar)
    80       real mask(imar+1,jmar)
     80      real mask(imar+1,jmar), mask_tmp(imar+1,jmar)
    8181      real num_tot(2200,1100),num_lan(2200,1100)
    8282c
     
    300300       CALL MVA9(zxtzy,iim+1,jjm+1)
    301301c$$$       CALL MVA9(zytzy,iim+1,jjm+1)
     302C$$$   Masque prenant en compte maximum de terre
     303       mask_tmp= 0.0
     304       WHERE(mask .GE. epsfra) mask_tmp = 1.
    302305
    303306       DO ii = 1, imar
     
    314317           if(xq.le.xw) xq=xw
    315318           if(abs(xm).le.xw) xm=xw*sign(1.,xm)
    316 c$$$c slope:
     319c slope:
    317320c$$$           zsig(ii,jj)=sqrt(xq)*mask(ii,jj)
    318321c$$$c isotropy:
     
    327330C$$* PB modif pour maque de terre fractionnaire
    328331c slope:
    329            zsig(ii,jj)=sqrt(xq)
     332           zsig(ii,jj)=sqrt(xq)*mask_tmp(ii,jj)
    330333c isotropy:
    331            zgam(ii,jj)=xp/xq
     334           zgam(ii,jj)=xp/xq*mask_tmp(ii,jj)
    332335c angle theta:
    333            zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.
    334            zphi(ii,jj)=zmea(ii,jj)
    335 c$$$           zmea(ii,jj)=zmea(ii,jj)
    336 c$$$           zpic(ii,jj)=zpic(ii,jj)
    337 c$$$           zval(ii,jj)=zval(ii,jj)
    338 c$$$           zstd(ii,jj)=zstd(ii,jj)
     336           zthe(ii,jj)=57.29577951*atan2(xm,xl)/2.*mask_tmp(ii,jj)
     337           zphi(ii,jj)=zmea(ii,jj)*mask_tmp(ii,jj)
     338           zmea(ii,jj)=zmea(ii,jj)*mask_tmp(ii,jj)
     339           zpic(ii,jj)=zpic(ii,jj)*mask_tmp(ii,jj)
     340           zval(ii,jj)=zval(ii,jj)*mask_tmp(ii,jj)
     341           zstd(ii,jj)=zstd(ii,jj)*mask_tmp(ii,jj)
    339342c          print 101,ii,jj,
    340343c    *           zmea(ii,jj),zstd(ii,jj),zsig(ii,jj),zgam(ii,jj),
Note: See TracChangeset for help on using the changeset viewer.