Changeset 1171


Ignore:
Timestamp:
2008-09-10T17:30:28+02:00 (12 years ago)
Author:
ctlod
Message:

now input air temperature data for CORE & CLIO bulks must be in Kelvin, see ticket: #242

Location:
trunk/NEMO
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limthd.F90

    r1156 r1171  
    123123      !------------------------------------------------------------------------------! 
    124124      zeps = 1.0e-10 
    125       tatm_ice(:,:) = tatm_ice(:,:) + 273.15 ! convert C to K 
    126125 
    127126      !-------------------- 
  • trunk/NEMO/LIM_SRC_3/limthd_dh.F90

    r1156 r1171  
    275275         qt_s_in(ji,jl) =  qt_s_in(ji,jl) + zqprec(ji) * zdh_s_pre(ji) 
    276276         zqt_s(ji)      =  zqt_s(ji) + zqprec(ji) * zdh_s_pre(ji) 
    277       END DO 
    278  
    279       ! Update total snow heat content 
    280       zqt_s(ji)         =  MAX ( zqt_s(ji) - zqfont_su(ji) , 0.0 )  
    281       IF( lk_mpp ) CALL mpp_max(zqt_s(ji), kcom = ncomm_ice )  
     277         zqt_s(ji)      =  MAX ( zqt_s(ji) - zqfont_su(ji) , 0.0 )  
     278      END DO 
     279 
    282280 
    283281      ! Snow melt due to surface heat imbalance 
  • trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

    r1152 r1171  
    295295         DO ji = 1, jpi 
    296296            ! 
    297             zsst  = pst(ji,jj)              + rt0           ! converte Celcius to Kelvin the SST and Tair 
    298             ztatm = sf(jp_tair)%fnow(ji,jj) + rt0           ! and set minimum value far above 0 K (=rt0 over land) 
     297            zsst  = pst(ji,jj)              + rt0           ! converte Celcius to Kelvin the SST 
     298            ztatm = sf(jp_tair)%fnow(ji,jj)                 ! and set minimum value far above 0 K (=rt0 over land) 
    299299            zcco1 = 1.0 - sf(jp_ccov)%fnow(ji,jj)           ! fraction of clear sky ( 1 - cloud cover) 
    300300            zrhoa = zpatm / ( 287.04 * ztatm )              ! air density (equation of state for dry air)  
     
    495495!CDIR NOVERRCHK 
    496496         DO ji = 1, jpi 
    497             ztatm (ji,jj) = sf(jp_tair)%fnow(ji,jj) + rt0            ! air temperature in Kelvins  
     497            ztatm (ji,jj) = sf(jp_tair)%fnow(ji,jj)                  ! air temperature in Kelvins  
    498498       
    499499            zrhoa(ji,jj) = zpatm / ( 287.04 * ztatm(ji,jj) )         ! air density (equation of state for dry air)  
     
    703703!CDIR NOVERRCHK 
    704704         DO ji = 1, jpi 
    705             ztamr = sf(jp_tair)%fnow(ji,jj) + rt0 - rtt 
     705            ztamr = sf(jp_tair)%fnow(ji,jj) - rtt 
    706706            zmt1  = SIGN( 17.269,  ztamr ) 
    707707            zmt2  = SIGN( 21.875,  ztamr ) 
    708708            zmt3  = SIGN( 28.200, -ztamr ) 
    709709            zes = 611.0 * EXP(  ABS( ztamr ) * MIN ( zmt1, zmt2 )   &              ! Saturation water vapour 
    710                &                     / ( sf(jp_tair)%fnow(ji,jj) + rt0 - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
     710               &                     / ( sf(jp_tair)%fnow(ji,jj) - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
    711711            zev(ji,jj) = sf(jp_humi)%fnow(ji,jj) * zes * 1.0e-05                   ! vapour pressure   
    712712         END DO 
     
    833833!CDIR NOVERRCHK 
    834834         DO ji = 1, jpi            
    835             ztamr = sf(jp_tair)%fnow(ji,jj) + rt0 - rtt            
     835            ztamr = sf(jp_tair)%fnow(ji,jj) - rtt            
    836836            zmt1  = SIGN( 17.269,  ztamr ) 
    837837            zmt2  = SIGN( 21.875,  ztamr ) 
    838838            zmt3  = SIGN( 28.200, -ztamr ) 
    839839            zes = 611.0 * EXP(  ABS( ztamr ) * MIN ( zmt1, zmt2 )   &              ! Saturation water vapour 
    840                &                     / ( sf(jp_tair)%fnow(ji,jj) + rt0 - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
     840               &                     / ( sf(jp_tair)%fnow(ji,jj) - 35.86  + MAX( 0.e0, zmt3 ) )  ) 
    841841            zev(ji,jj) = sf(jp_humi)%fnow(ji,jj) * zes * 1.0e-05                   ! vapour pressure   
    842842         END DO 
Note: See TracChangeset for help on using the changeset viewer.