New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15473 for branches/UKMO/dev_r5518_GO6_package_FOAMv14_dsst_CO2/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2021-11-05T12:40:18+01:00 (3 years ago)
Author:
dford
Message:

Use cool skin SST in air-sea CO2 flux calculations.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_dsst_CO2/NEMOGCM/NEMO/OPA_SRC/step.F90

    r13355 r15473  
    253253      ENDIF 
    254254 
     255      ! Cool skin 
     256      IF ( ln_diurnal ) THEN  
     257         IF ( ln_blk_core ) THEN 
     258            CALL diurnal_sst_coolskin_step(qns(:,:)+(rn_abs*qsr(:,:)), taum, rhop(:,:,1), rdt) 
     259         ELSE 
     260            CALL diurnal_sst_coolskin_step(qns, taum, rhop(:,:,1), rdt) 
     261         ENDIF 
     262      ENDIF 
     263 
    255264      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    256265      ! diagnostics and outputs             (ua, va, tsa used as workspace) 
     
    266275      ! 
    267276      IF( ln_crs     )      CALL crs_fld( kstp )         ! ocean model: online field coarsening & output 
     277 
     278      !Diurnal warm layer model        
     279      IF ( ln_diurnal ) THEN 
     280         IF ( ln_blk_core ) THEN 
     281            IF( kstp == nit000 )THEN  
     282               CALL diurnal_sst_takaya_step( &  
     283                  &    qsr(:,:)-(rn_abs*qsr(:,:)), qns(:,:)+(rn_abs*qsr(:,:)), & 
     284                  &    taum, rhop(:,:,1), & 
     285                  &    rdt, ld_calcfrac = .TRUE.)  
     286            ELSE  
     287               CALL diurnal_sst_takaya_step( &  
     288                  &    qsr(:,:)-(rn_abs*qsr(:,:)), qns(:,:)+(rn_abs*qsr(:,:)), & 
     289                  &    taum, rhop(:,:,1), rdt )  
     290            ENDIF 
     291         ELSE 
     292            IF( kstp == nit000 )THEN  
     293               CALL diurnal_sst_takaya_step( &  
     294                  &    qsr, qns, taum, rhop(:,:,1), & 
     295                  &    rdt, ld_calcfrac = .TRUE.)  
     296            ELSE  
     297               CALL diurnal_sst_takaya_step( &  
     298                  &    qsr, qns, taum, rhop(:,:,1), rdt )  
     299            ENDIF 
     300         ENDIF 
     301      ENDIF 
    268302 
    269303#if defined key_top 
Note: See TracChangeset for help on using the changeset viewer.