Changeset 102 for codes/icosagcm/trunk


Ignore:
Timestamp:
08/06/12 03:13:10 (12 years ago)
Author:
ymipsl
Message:

implement virtual temperature for moist physics from dcmip

YM

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r99 r102  
    1 MODULE physics_dcmip_mod 
     1        MODULE physics_dcmip_mod 
    22  USE ICOSA 
    33  PRIVATE 
    … …  
    133133    vfi=v     
    134134     
     135    DO l=1,llm 
     136      DO j=jj_begin,jj_end 
     137        DO i=ii_begin,ii_end 
     138          ij=(j-1)*iim+i 
     139          T(ij,l)=T(ij,l)/(1+0.608*q(ij,l)) 
     140        ENDDO 
     141      ENDDO 
     142    ENDDO        
     143     
     144     
    135145    CALL simple_physics(iim*jjm, llm, dt, lat, t, q , ufi, vfi, pmid, p, pdel, 1/pdel, ps, precl, testcase)  
    136146     
     147    DO l=1,llm 
     148      DO j=jj_begin,jj_end 
     149        DO i=ii_begin,ii_end 
     150          ij=(j-1)*iim+i 
     151          T(ij,l)=T(ij,l)*(1+0.608*q(ij,l)) 
     152        ENDDO 
     153      ENDDO 
     154    ENDDO        
    137155         
    138156    ufi=ufi-u 
Note: See TracChangeset for help on using the changeset viewer.