Changeset 7034

Timestamp:

2016-10-17T18:02:09+02:00 (8 years ago)

Author:

kuniko

Message:

Initial commit

File:

• branches/UKMO/dev_r5518_fa_am_dt_deltadelta_toa/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90 (modified) (10 diffs)

branches/UKMO/dev_r5518_fa_am_dt_deltadelta_toa/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90

@@ -32 +32 @@
    PUBLIC   sbc_flx_adj_init   ! routine called in sbcmod
 
+   !ky 16/10/2016 commented out below two lines
    !ky 06/09/2016 uncommented below two lines for FA test
    !ky 11/12/2015 recommented below two lines
    !!ky 3/12/2015 uncommented below two lines for FA test!
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s]
-   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2]
+   !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s]
+   !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2]
 
    !                                   !!* Namelist namsbc_flx_adj *
@@ -43 +44 @@
    LOGICAL         ::   ln_sssr_bnd     ! flag to bound erp term 
    REAL(wp)        ::   rn_sssr_bnd     ! ABS(Max./Min.) value of erp term [mm/day]
+   REAL(wp)        ::   delta_toa       ! global mean annual mean TOA for FA to perturbed ensemble [W/m2]
+   REAL(wp)        ::   delta_delta_toa ! perturbation to delta_toa [W/m2]
 
    REAL(wp) , ALLOCATABLE, DIMENSION(:) ::   buffer   ! Temporary buffer for exchange
@@ -83 +86 @@
       !!----------------------------------------------------------------------
       !
-      !!write(numout,*) '*** In sbcflx_adj ***'
+      write(numout,*) '*** In sbcflx_adj *** delta_toa, delta_delta_toa=', &
+         & delta_toa, delta_delta_toa
 
       IF( nn_timing == 1 )  CALL timing_start('sbc_flx_adj')
@@ -101 +105 @@
                      ! K.Y. 16/03/2016 Apply Arctic mask to heat flux adjustment
                      zqrp = ( 1. - 2.*only_arcmsk(ji,jj) )  &
-                        &        * sf_qrp(1)%fnow(ji,jj,1)
+                        &        * sf_qrp(1)%fnow(ji,jj,1)  &
+                        &        + delta_toa + delta_delta_toa
+                     ! K.Y. 16/10/2016 added "&" to 2 lines above and 
+                     !   delta_toa, delta_delta_toa to above line
                      !zqrp = sf_qrp(1)%fnow(ji,jj,1)
                      !ky 07/09/2016 copied 1 line below for FA test
@@ -108 +115 @@
                      !qrp(ji,jj) = sf_qrp(1)%fnow(ji,jj,1)
                      !ky 07/09/2016 copied 1 line above and modified as below for FA test
-                     qrp(ji,jj) = zqrp
+                     !ky 16/10/2016 commented out 1 line below
+                     !qrp(ji,jj) = zqrp
                      qns(ji,jj) = qns(ji,jj) + zqrp
                   END DO
                END DO
+               !ky 16/10/2016 commented out below 1 line below
                !ky 06/09/2016 uncommented below 1 line below for FA test!
                !ky 11/12/2015 recommented out below 1 line below
                !!ky 3/12/2015 uncommented below 1 line below for FA test!
-               CALL iom_put( "qrp", qrp )                             ! heat flux damping
+               !CALL iom_put( "qrp", qrp )                             ! heat flux damping
             ENDIF
             !
@@ -132 +141 @@
                      !!ky 3/12/2015 1 line below for FA test!
                      !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1)
+                     !ky 16/10/2016 commented out below two lines
                      !ky 06/09/2016 copied line above and changed to below two lines for FA test!
-                     erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) )  &        ! No damping in vicinity of river mouths
-                        &        * sf_erp(1)%fnow(ji,jj,1) 
+                     !erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) )  &        ! No damping in vicinity of river mouths
+                     !   &        * sf_erp(1)%fnow(ji,jj,1) 
                      sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux
                   END DO
                END DO
+               !ky 16/10/2016 commented out below one line
                !ky 06/09/2016 uncommented below one line for FA test!
                !ky 11/12/2015 recommented out below one line
                !!ky 3/12/2015 uncommented below one line for FA test!
-               CALL iom_put( "erp", erp )                             ! freshwater flux damping
+               !CALL iom_put( "erp", erp )                             ! freshwater flux damping
                !
             ELSEIF( nn_flxadjfw == 2 ) THEN                               !* Salinity damping term (volume flux (emp) and associated heat flux (qns)
@@ -155 +166 @@
                      !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1)
                      !ky 06/09/2016 copied line above and changed to below line for FA test!
-                     erp(ji,jj) = zerp
+                     !ky 16/10/2016 commented out below line
+                     !erp(ji,jj) = zerp
                      emp(ji,jj) = emp (ji,jj) + zerp
                      qns(ji,jj) = qns(ji,jj) - zerp * rcp * sst_m(ji,jj)
                   END DO
                END DO
+               !ky 16/10/2016 commented out below 1 line
                !ky 06/09/2016 uncommented below 1 line for FA test!
                !ky 11/12/2015 recommented out below 1 line
                !!ky 3/12/2015 uncommented below 1 line for FA test!
-               CALL iom_put( "erp", erp )                             ! freshwater flux damping
+               !CALL iom_put( "erp", erp )                             ! freshwater flux damping
             ENDIF
             !
@@ -193 +206 @@
       TYPE(FLD_N) ::   sn_qrp, sn_erp        ! informations about the fields to be read
       NAMELIST/namsbc_flx_adj/ cn_dir, nn_flxadjht, nn_flxadjfw, sn_qrp, sn_erp,   &
-         &                     ln_sssr_bnd, rn_sssr_bnd
+         &                     ln_sssr_bnd, rn_sssr_bnd, delta_toa, delta_delta_toa
       INTEGER     ::  ios
       !!----------------------------------------------------------------------
@@ -219 +232 @@
          WRITE(numout,*) '      flag to bound erp term                 ln_sssr_bnd = ', ln_sssr_bnd
          WRITE(numout,*) '      ABS(Max./Min.) erp threshold           rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day'
+         WRITE(numout,*) '      global mean, annual mean TOA to supplement FA delta_toa = ', delta_toa, ' W/m2'
+         WRITE(numout,*) '      perturbation to delta_toa = ', delta_delta_toa, ' W/m2'
       ENDIF
       !
       !                            !* Allocate erp and qrp array
+      !ky 16/10/2016 commented out below two lines
       !ky 06/09/2016 uncommented below two lines for FA test!
       !ky 11/12/2015 recommented out below two lines
       !!ky 3/12/2015 uncommented below two lines for FA test!
-      ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror )
-      IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' )
+      !ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror )
+      !IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' )
       !
       IF( nn_flxadjht == 1 ) THEN      !* set sf_qrp structure & allocate arrays
@@ -257 +273 @@
       !
       !ky!!                            !* Initialize qrp and erp if no restoring 
+      !ky 16/10/2016 commented out below two lines
       !ky 06/09/2016 below two lines for FA test!
       !ky 11/12/2015 commented out below two lines
       !!ky 3/12/2015 below two lines for FA test!
-      IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp
-      IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp
+      !IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp
+      !IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp
       !!ky!!IF( nn_sstr /= 1                   )   qrp(:,:) = 0._wp
       !!ky!!IF( nn_sssr /= 1 .OR. nn_sssr /= 2 )   erp(:,:) = 0._wp