Changeset 550


Ignore:
Timestamp:
09/18/17 13:58:05 (7 years ago)
Author:
dubos
Message:

devel : backported commits 541 and 545 from trunk

Location:
codes/icosagcm/devel/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/devel/src/dissip/dissip_gcm.f90

    r533 r550  
    8787   
    8888             
    89   INTEGER :: i,j,ij,ind,it,iter 
     89  INTEGER :: i,j,ij,ind,it,iter,M 
    9090 
    9191   rayleigh_friction_key='none' 
     
    159159      u=f_u(ind) 
    160160 
     161! set random seed to get reproductibility when using a different number of process 
     162      CALL RANDOM_SEED(size=M) 
     163      CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 
     164 
    161165      DO j=jj_begin,jj_end 
    162166        DO i=ii_begin,ii_end 
     
    247251      u=f_u(ind) 
    248252 
    249       DO j=jj_begin,jj_end 
     253! set random seed to get reproductibility when using a different number of process 
     254      CALL RANDOM_SEED(size=M) 
     255      CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 
     256  
     257       DO j=jj_begin,jj_end 
    250258        DO i=ii_begin,ii_end 
    251259          ij=(j-1)*iim+i    
     
    335343      CALL swap_geometry(ind) 
    336344      theta=f_theta(ind) 
    337   
     345 
     346! set random seed to get reproductibility when using a different number of process 
     347      CALL RANDOM_SEED(size=M) 
     348      CALL RANDOM_SEED(put=(/(i,i=1,M)/)) 
     349 
    338350      DO j=jj_begin,jj_end 
    339351        DO i=ii_begin,ii_end 
  • codes/icosagcm/devel/src/vertical/disvert_strato.f90

    r531 r550  
    1212  SUBROUTINE init_disvert 
    1313  USE icosa 
     14  USE getin_mod 
    1415  IMPLICIT NONE 
    15    
     16  REAL(rstd) :: dsigmin 
    1617    ALLOCATE(ap(llm+1)) 
    1718    ALLOCATE(bp(llm+1)) 
    1819    ALLOCATE(presnivs(llm)) 
    1920 
    20     CALL disvert(ap,bp,presnivs)     
     21    dsigmin=1.0  ! Should be 0.3 for CMIP5 
     22    CALL getin('disvert_dsigmin', dsigmin)  
     23    CALL disvert(dsigmin,ap,bp,presnivs) 
    2124 
    2225  END SUBROUTINE init_disvert   
    2326 
    2427 
    25   SUBROUTINE disvert(ap,bp,presnivs) 
     28  SUBROUTINE disvert(dsigmin,ap,bp,presnivs) 
    2629  USE icosa 
    2730  USE mpipara 
    2831  USE earth_const 
    2932  IMPLICIT NONE 
     33  REAL(rstd) :: dsigmin 
    3034  REAL(rstd),INTENT(OUT) :: ap(:) 
    3135  REAL(rstd),INTENT(OUT) :: bp(:) 
     
    4145    DO l = 1, llm 
    4246      x = Pi*(l-0.5)/(llm+1) 
    43       dsig(l) = (1.0 + 7.0 * SIN(x)**2) * (0.5* (1.-tanh( (x-Pi/2)/(Pi/2) ) ))**2 
     47      dsig(l) = (dsigmin + 7.0 * SIN(x)**2) * (0.5* (1.-tanh( (x-Pi/2)/(Pi/2) ) ))**2 
    4448      snorm = snorm + dsig(l) 
    4549    ENDDO     
Note: See TracChangeset for help on using the changeset viewer.