Changeset 11997


Ignore:
Timestamp:
2019-11-28T12:49:14+01:00 (11 months ago)
Author:
laurent
Message:

Better namelist sanity test in sbcblk.F90 w.r.t algo choice.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90

    r11996 r11997  
    207207 
    208208 
    209       !LB: 
     209 
    210210      !                             !** initialization of the cool-skin / warm-layer parametrization 
    211       IF( ln_NCAR .AND. (ln_skin_cs .OR. ln_skin_wl) ) & 
    212          & CALL ctl_stop( 'sbc_blk_init: Cool-skin/warm-layer param. not compatible with NCAR algorithm!' ) 
    213       ! 
     211      IF( ln_skin_cs .OR. ln_skin_wl ) THEN 
     212         !! Some namelist sanity tests: 
     213         IF( ln_NCAR )      & 
     214            & CALL ctl_stop( 'sbc_blk_init: Cool-skin/warm-layer param. not compatible with NCAR algorithm' )          
     215         IF( nn_fsbc /= 1 ) & 
     216            & CALL ctl_stop( 'sbc_blk_init: Please set "nn_fsbc" to 1 when using cool-skin/warm-layer param.')                   
     217      END IF 
     218       
    214219      ioptio = 0 
    215220      IF( ln_humi_sph ) THEN 
     
    532537         CALL ctl_stop( 'STOP', 'sbc_oce: non-existing bulk formula selected' ) 
    533538      END SELECT 
    534  
     539       
    535540      IF( ln_skin_cs .OR. ln_skin_wl ) THEN 
    536541         !! In the presence of sea-ice we forget about the cool-skin/warm-layer update of zst and zsq: 
     
    546551         tsk(:,:) = zst(:,:) !#LB: Update of tsk, the "official" array for skin temperature 
    547552      END IF 
    548        
    549       !ELSE !IF( ln_skin_cs .OR. ln_skin_wl ) 
    550       ! 
    551       !   SELECT CASE( nblk )        !==  transfer coefficients  ==!   Cd, Ch, Ce at T-point 
    552       ! 
    553       !CASE( np_COARE_3p0 ) 
    554       !      CALL turb_coare3p0 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 
    555       !         &                Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 
    556       ! 
    557       !   CASE( np_COARE_3p6 ) 
    558       !CALL turb_coare3p6 ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl, & 
    559       !         &                Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 
    560       ! 
    561       !   CASE( np_ECMWF     ) 
    562       !      CALL turb_ecmwf   ( kt, rn_zqt, rn_zu, zst, ztpot, zsq, zqair, wndm, ln_skin_cs, ln_skin_wl,  & 
    563       !         &                Cd_atm, Ch_atm, Ce_atm, t_zu, q_zu, zU_zu, cdn_oce, chn_oce, cen_oce ) 
    564       !      !  
    565       !   CASE DEFAULT 
    566       !      CALL ctl_stop( 'STOP', 'sbc_oce: non-existing bulk formula selected' ) 
    567       !   END SELECT! 
    568       ! 
    569       !ENDIF ! IF( ln_skin_cs .OR. ln_skin_wl ) 
    570553 
    571554      !!      CALL iom_put( "Cd_oce", Cd_atm)  ! output value of pure ocean-atm. transfer coef. 
Note: See TracChangeset for help on using the changeset viewer.