Changeset 13241 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

Timestamp:

2020-07-03T14:42:49+02:00 (4 years ago)

Author:

dford

Message:

Update NEMO-FABM coupler for compatability with FABM v1.0.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90 (modified) (7 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_utils366_fabmv1/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

@@ -45 +45 @@
    LOGICAL , PUBLIC ::   ln_qsr_bio   !: bio-model      light absorption flag
    LOGICAL , PUBLIC ::   ln_qsr_ice   !: light penetration for ice-model LIM3 (clem)
+   LOGICAL , PUBLIC ::   ln_qsr_spec  !: spectral model heating from ERSEM
    INTEGER , PUBLIC ::   nn_chldta    !: use Chlorophyll data (=1) or not (=0)
    INTEGER , PUBLIC ::   nn_kd490dta  !: use kd490dta data (=1) or not (=0)
@@ -150 +151 @@
       
       !                                           ! ============================================== !
-      IF( lk_qsr_bio .AND. ln_qsr_bio ) THEN      !  bio-model fluxes  : all vertical coordinates  !
+      IF( ln_qsr_spec ) THEN                      !  ERSEM spectral heating                        !
+         !                                        ! ============================================== !
+         DO jk = 1, jpkm1
+           qsr_hc(:,:,jk) = r1_rau0_rcp * ( etot3(:,:,jk) )
+         END DO
+         !                                        Add to the general trend
+         DO jk = 1, jpkm1
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  z1_e3t = zfact / fse3t(ji,jj,jk)
+                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) + ( qsr_hc_b(ji,jj,jk) + qsr_hc(ji,jj,jk) ) * z1_e3t
+               END DO
+            END DO
+         END DO
+         CALL iom_put( 'qsr3d', etot3 )   ! Shortwave Radiation 3D distribution
+         IF ( ln_qsr_ice ) THEN
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  IF ( qsr(ji,jj) /= 0._wp ) THEN
+                     fraqsr_1lev(ji,jj) = ( qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) ) )
+                  ELSE
+                     fraqsr_1lev(ji,jj) = 1.
+                  ENDIF
+               END DO
+            END DO
+         ENDIF
+         ! 
+
+         !                                        ! ============================================== !
+      ELSEIF( lk_qsr_bio .AND. ln_qsr_bio ) THEN      !  bio-model fluxes  : all vertical coordinates  !
          !                                        ! ============================================== !
          DO jk = 1, jpkm1
@@ -423 +453 @@
       !!
       NAMELIST/namtra_qsr/  sn_chl, sn_kd490, cn_dir, ln_traqsr, ln_qsr_rgb, ln_qsr_2bd, ln_qsr_bio, ln_qsr_ice,  &
-         &                  nn_chldta, rn_abs, rn_si0, rn_si1, nn_kd490dta
+         &                  ln_qsr_spec, nn_chldta, rn_abs, rn_si0, rn_si1, nn_kd490dta
       !!----------------------------------------------------------------------
 
@@ -451 +481 @@
          WRITE(numout,*) '      2 band               light penetration   ln_qsr_2bd = ', ln_qsr_2bd
          WRITE(numout,*) '      bio-model            light penetration   ln_qsr_bio = ', ln_qsr_bio
+         WRITE(numout,*) '      ERSEM spectral heating model             ln_qsr_spec= ', ln_qsr_spec
          WRITE(numout,*) '      light penetration for ice-model LIM3     ln_qsr_ice = ', ln_qsr_ice
          WRITE(numout,*) '      RGB : Chl data (=1) or cst value (=0)    nn_chldta  = ', nn_chldta
@@ -470 +501 @@
          IF( ln_qsr_2bd  )   ioptio = ioptio + 1
          IF( ln_qsr_bio  )   ioptio = ioptio + 1
+         IF( ln_qsr_spec )   ioptio = ioptio + 1
          IF( nn_kd490dta == 1 )   ioptio = ioptio + 1
          !
@@ -481 +513 @@
          IF( ln_qsr_bio                      )   nqsr =  4
          IF( nn_kd490dta == 1                )   nqsr =  5
+         IF( ln_qsr_spec                     )   nqsr =  6
          !
          IF(lwp) THEN                   ! Print the choice
@@ -489 +522 @@
             IF( nqsr ==  4 )   WRITE(numout,*) '         bio-model light penetration'
             IF( nqsr ==  5 )   WRITE(numout,*) '         KD490 light penetration'
-         ENDIF
+            IF( nqsr ==  6 )   WRITE(numout,*) '         ERSEM spectral light penetration'
+         ENDIF
+#if ! defined key_fabm
+         !
+         IF( nqsr ==  6 ) THEN
+            CALL ctl_stop( 'ln_qsr_spec=.true. so trying to use ERSEM spectral light penetration', &
+               &           'but not running with ERSEM' )
+         ENDIF
+#endif
          !
       ENDIF