Changeset 15433

Timestamp:

2021-10-22T11:16:15+02:00 (3 years ago)

Author:

dancopsey

Message:

Add evap only option

File:

• NEMO/branches/UKMO/NEMO_4.0.1_NGMS_couple_dan/src/OCE/SBC/sbccpl.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_NGMS_couple_dan/src/OCE/SBC/sbccpl.F90

@@ -463 +463 @@
          IF ( k_ice <= 1 )  srcv(jpr_ievp)%laction = .FALSE.
       CASE( 'ice only'  )       ;   srcv(jpr_ievp)%laction = .TRUE. 
+      CASE( 'evap only'  )      ;   srcv( (/jpr_ievp, jpr_tevp/) )%laction = .TRUE.
       CASE( 'oce and ice'   )   ;   srcv( (/jpr_ievp, jpr_sbpr, jpr_semp, jpr_oemp/) )%laction = .TRUE.
       CASE default              ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_rcv_emp%cldes' )
@@ -1499 +1500 @@
             CASE( 'oce only', 'oce and ice' )
                zemp(:,:) = frcv(jpr_oemp)%z3(:,:,1)
+            CASE( 'evap only' )
+               zemp(:,:) = frcv(jpr_tevp)%z3(:,:,1)
             CASE default
                CALL ctl_stop( 'sbc_cpl_rcv: wrong definition of sn_rcv_emp%cldes' )
@@ -1775 +1778 @@
          zemp_tot(:,:) =   frcv(jpr_tevp)%z3(:,:,1) - ztprecip(:,:)
          zemp_ice(:,:) = ( frcv(jpr_ievp)%z3(:,:,1) - frcv(jpr_snow)%z3(:,:,1) ) * picefr(:,:)
+      CASE( 'evap only' )   ! received fields: jpr_ievp, jpr_tevp
+         zsprecip(:,:) =   0._wp                  ! May need to ensure positive here
+         ztprecip(:,:) =   0._wp  ! May need to ensure positive here
+         zemp_tot(:,:) =   frcv(jpr_tevp)%z3(:,:,1)
+         zemp_ice(:,:) =   frcv(jpr_ievp)%z3(:,:,1)
       CASE( 'oce and ice'   )   ! received fields: jpr_sbpr, jpr_semp, jpr_oemp, jpr_ievp
          zemp_tot(:,:) = ziceld(:,:) * frcv(jpr_oemp)%z3(:,:,1) + picefr(:,:) * frcv(jpr_sbpr)%z3(:,:,1)