Changeset 4708 for branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2014-07-10T17:17:46+02:00 (7 years ago)

Author:

rfurner

Message:

changes to remove tracer modules and produce barotropic model

File:

• branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (6 diffs)

branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -303 +303 @@
       ! ocean albedo assumed to be constant + modify now Qsr to include the diurnal cycle                    ! Short Wave
       zztmp = 1. - albo
+#if defined key_surge
+      qsr(:,:)=0._wp
+      zqlw(:,:) = 0._wp
+#else
       IF( ln_dm2dc ) THEN   ;   qsr(:,:) = zztmp * sbc_dcy( sf(jp_qsr)%fnow(:,:,1) ) * tmask(:,:,1)
       ELSE                  ;   qsr(:,:) = zztmp *          sf(jp_qsr)%fnow(:,:,1)   * tmask(:,:,1)
@@ -308 +312 @@
 !CDIR COLLAPSE
       zqlw(:,:) = (  sf(jp_qlw)%fnow(:,:,1) - Stef * zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:)  ) * tmask(:,:,1)   ! Long  Wave
+#endif
       ! ----------------------------------------------------------------------------- !
       !     II    Turbulent FLUXES                                                    !
@@ -382 +387 @@
       !  Turbulent fluxes over ocean
       ! -----------------------------
+#if ! defined key_surge
       IF( ln_2m .OR. ( ln_bulk2z .AND. rn_zqt /= rn_zu ) ) THEN
          ! Values of temp. and hum. adjusted to height of wind must be used
@@ -394 +400 @@
 !CDIR COLLAPSE
       zqla (:,:) = Lv * zevap(:,:)                                                              ! Latent Heat
+#endif
 
       IF(ln_ctl) THEN
@@ -410 +417 @@
       ! ----------------------------------------------------------------------------- !
      
+#if defined key_surge
+      emp (:,:) = 0._wp
+      qns(:,:)  = 0._wp
+#else
 !CDIR COLLAPSE
       emp (:,:) = (  zevap(:,:)                                          &   ! mass flux (evap. - precip.)
@@ -421 +432 @@
          &     + sf(jp_snow)%fnow(:,:,1) * rn_pfac                                &   ! add solid  precip heat content at min(Tair,Tsnow)
          &     * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic 
+#endif
       !
       CALL iom_put( "qlw_oce",   zqlw )                 ! output downward longwave heat over the ocean