Ignore:
Timestamp:
2014-07-10T17:17:46+02:00 (7 years ago)
Author:
rfurner
Message:

changes to remove tracer modules and produce barotropic model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r4624 r4708  
    303303      ! ocean albedo assumed to be constant + modify now Qsr to include the diurnal cycle                    ! Short Wave 
    304304      zztmp = 1. - albo 
     305#if defined key_surge 
     306      qsr(:,:)=0._wp 
     307      zqlw(:,:) = 0._wp 
     308#else 
    305309      IF( ln_dm2dc ) THEN   ;   qsr(:,:) = zztmp * sbc_dcy( sf(jp_qsr)%fnow(:,:,1) ) * tmask(:,:,1) 
    306310      ELSE                  ;   qsr(:,:) = zztmp *          sf(jp_qsr)%fnow(:,:,1)   * tmask(:,:,1) 
     
    308312!CDIR COLLAPSE 
    309313      zqlw(:,:) = (  sf(jp_qlw)%fnow(:,:,1) - Stef * zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:)  ) * tmask(:,:,1)   ! Long  Wave 
     314#endif 
    310315      ! ----------------------------------------------------------------------------- ! 
    311316      !     II    Turbulent FLUXES                                                    ! 
     
    382387      !  Turbulent fluxes over ocean 
    383388      ! ----------------------------- 
     389#if ! defined key_surge 
    384390      IF( ln_2m .OR. ( ln_bulk2z .AND. rn_zqt /= rn_zu ) ) THEN 
    385391         ! Values of temp. and hum. adjusted to height of wind must be used 
     
    394400!CDIR COLLAPSE 
    395401      zqla (:,:) = Lv * zevap(:,:)                                                              ! Latent Heat 
     402#endif 
    396403 
    397404      IF(ln_ctl) THEN 
     
    410417      ! ----------------------------------------------------------------------------- ! 
    411418      
     419#if defined key_surge 
     420      emp (:,:) = 0._wp 
     421      qns(:,:)  = 0._wp 
     422#else 
    412423!CDIR COLLAPSE 
    413424      emp (:,:) = (  zevap(:,:)                                          &   ! mass flux (evap. - precip.) 
     
    421432         &     + sf(jp_snow)%fnow(:,:,1) * rn_pfac                                &   ! add solid  precip heat content at min(Tair,Tsnow) 
    422433         &     * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic  
     434#endif 
    423435      ! 
    424436      CALL iom_put( "qlw_oce",   zqlw )                 ! output downward longwave heat over the ocean 
Note: See TracChangeset for help on using the changeset viewer.