Changeset 417


Ignore:
Timestamp:
06/10/16 22:28:02 (8 years ago)
Author:
dubos
Message:

Fixed computation of omega

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

Legend:

Unmodified
Added
Removed

  • codes/icosagcm/trunk/src/observable.f90

    r413 r417  
    5959!    IF (is_master) PRINT *,'CALL write_output_fields_basic' 
    6060 
     61    CALL transfert_request(f_ps,req_i1) 
     62     
    6163    IF(init) THEN 
    6264       scalar(1)=dt 

  • codes/icosagcm/trunk/src/omega.f90

    r295 r417  
    146146                  + le(ij+u_rdown)*u(ij+u_rdown,l)*( ne(ij,rdown)*ps(ij) + ne(ij+t_rdown,lup)*ps(ij+t_rdown) ) 
    147147             ugradps = .5*(bp(l)+bp(l+1)) *ugradps/(-4*Ai(ij))  ! sign convention as in Ringler et al. 2010, Eq. 22 p.3072 
    148              w( ij, l) =  ugradps - .5*(convm( ij,l+1)+convm(ij,l))  
     148             w( ij, l) =  ugradps - g*.5*(convm( ij,l+1)+convm(ij,l))  
    149149          ENDDO 
    150150        ENDDO 
Note: See TracChangeset for help on using the changeset viewer.