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 15313 – NEMO

Changeset 15313


Ignore:
Timestamp:
2021-10-01T13:57:20+02:00 (3 years ago)
Author:
hadjt
Message:

POLCOMS style radiatve transfer for radiative fluxes

TRA/tradwl.F90 added. step.F90 and step_oce.F90 edited.
Namelist added: namtra_dwl, with two keys: ln_tradwl, ln_vary_lambda

This version is still in progress, as the variable hbatt from CO6 is no in CO8, and I need to confirm my approach is
valid.

This was crashing due to rlambda2 being infinity, when depth was 0. I have set a max value of 0.25, which is
equivalent to about 10cm of water.

The version is committed in restrospect, so may not include all the changes to work.

Location:
NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/step.F90

    r14078 r15313  
    228228                         CALL tra_sbc       ( kstp )  ! surface boundary condition 
    229229      IF( ln_traqsr  )   CALL tra_qsr       ( kstp )  ! penetrative solar radiation qsr 
     230      !JT 
     231      IF( ln_tradwl  )   CALL tra_dwl       ( kstp )  ! Polcoms Style Short Wave Radiation  
     232      !JT 
     233 
    230234      IF( ln_trabbc  )   CALL tra_bbc       ( kstp )  ! bottom heat flux 
    231235      IF( ln_trabbl  )   CALL tra_bbl       ( kstp )  ! advective (and/or diffusive) bottom boundary layer scheme 
  • NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/step_oce.F90

    r14078 r15313  
    2323 
    2424   USE traqsr          ! solar radiation penetration      (tra_qsr routine) 
     25   !JT 
     26   USE tradwl          ! POLCOMS style solar radiation    (tra_dwl routine)  
     27   !JT 
    2528   USE trasbc          ! surface boundary condition       (tra_sbc routine) 
    2629   USE trabbc          ! bottom boundary condition        (tra_bbc routine) 
Note: See TracChangeset for help on using the changeset viewer.