Changeset 7169


Ignore:
Timestamp:
2016-11-02T19:56:47+01:00 (4 years ago)
Author:
kuniko
Message:

Corrected to only apply delta_toa & delta_delta_toa to heat flux adjusted ocean points

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_fa_am_dt_deltadelta_toa/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx_adj.F90

    r7057 r7169  
    3232   PUBLIC   sbc_flx_adj_init   ! routine called in sbcmod 
    3333 
     34   !ky 01/11/2016 uncommented below two lines for test 
    3435   !ky 16/10/2016 commented out below two lines 
    3536   !ky 06/09/2016 uncommented below two lines for FA test 
    3637   !ky 11/12/2015 recommented below two lines 
    3738   !!ky 3/12/2015 uncommented below two lines for FA test! 
    38    !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s] 
    39    !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2] 
     39   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s] 
     40   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2] 
    4041 
    4142   !                                   !!* Namelist namsbc_flx_adj * 
     
    106107                     ! K.Y. 16/03/2016 Apply Arctic mask to heat flux adjustment 
    107108                     zqrp = ( 1. - 2.*only_arcmsk(ji,jj) )  & 
    108                         &        * sf_qrp(1)%fnow(ji,jj,1)  & 
    109                         &        + delta_toa + delta_delta_toa 
     109                        &        * (sf_qrp(1)%fnow(ji,jj,1)  & 
     110                        &        + delta_toa + delta_delta_toa) 
     111                     !   &        * sf_qrp(1)%fnow(ji,jj,1)  & 
     112                     !   &        + delta_toa + delta_delta_toa 
     113                     ! K.Y. 02/11/2016 recoded two lines above to only add delta_toas to FA regions 
    110114                     ! K.Y. 16/10/2016 added "&" to 2 lines above and  
    111115                     !   delta_toa, delta_delta_toa to above line 
     
    117121                     !ky 07/09/2016 copied 1 line above and modified as below for FA test 
    118122                     !ky 16/10/2016 commented out 1 line below 
    119                      !qrp(ji,jj) = zqrp 
     123                     !ky 01/11/2016 uncommented 1 line below for test 
     124                     qrp(ji,jj) = zqrp 
    120125                     qns(ji,jj) = qns(ji,jj) + zqrp 
    121126                  END DO 
    122127               END DO 
     128               !ky 01/11/2016 uncommented 1 line below for test 
    123129               !ky 16/10/2016 commented out below 1 line below 
    124130               !ky 06/09/2016 uncommented below 1 line below for FA test! 
    125131               !ky 11/12/2015 recommented out below 1 line below 
    126132               !!ky 3/12/2015 uncommented below 1 line below for FA test! 
    127                !CALL iom_put( "qrp", qrp )                             ! heat flux damping 
     133               CALL iom_put( "qrp", qrp )                             ! heat flux damping 
     134               !ky 01/11/2016 added 1 line below for test 
     135               CALL iom_put( "only_arcmsk", only_arcmsk )                             ! arctic mask 
    128136            ENDIF 
    129137            ! 
     
    142150                     !!ky 3/12/2015 1 line below for FA test! 
    143151                     !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 
     152                     !ky 01/11/2016 uncommented two lines below for test. erp is freshwater flux. 
    144153                     !ky 16/10/2016 commented out below two lines 
    145154                     !ky 06/09/2016 copied line above and changed to below two lines for FA test! 
    146                      !erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) )  &        ! No damping in vicinity of river mouths 
    147                      !   &        * sf_erp(1)%fnow(ji,jj,1)  
     155                     erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) )  &        ! No damping in vicinity of river mouths 
     156                        &        * sf_erp(1)%fnow(ji,jj,1)  
    148157                     sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux 
    149158                  END DO 
    150159               END DO 
     160               !ky 01/11/2016 uncommented below one line for test 
    151161               !ky 16/10/2016 commented out below one line 
    152162               !ky 06/09/2016 uncommented below one line for FA test! 
    153163               !ky 11/12/2015 recommented out below one line 
    154164               !!ky 3/12/2015 uncommented below one line for FA test! 
    155                !CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     165               CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     166               !ky 01/11/2016 added 1 line below for test 
     167               CALL iom_put( "rnfmsk_arcmsk", rnfmsk_arcmsk )                         ! river mough & arctic  
    156168               ! 
    157169            ELSEIF( nn_flxadjfw == 2 ) THEN                               !* Salinity damping term (volume flux (emp) and associated heat flux (qns) 
     
    167179                     !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 
    168180                     !ky 06/09/2016 copied line above and changed to below line for FA test! 
     181                     !ky 01/11/2016 uncommented below line for test 
    169182                     !ky 16/10/2016 commented out below line 
    170                      !erp(ji,jj) = zerp 
     183                     erp(ji,jj) = zerp 
    171184                     emp(ji,jj) = emp (ji,jj) + zerp 
    172185                     qns(ji,jj) = qns(ji,jj) - zerp * rcp * sst_m(ji,jj) 
    173186                  END DO 
    174187               END DO 
     188               !ky 01/11/2016 uncommented below 1 line for test 
    175189               !ky 16/10/2016 commented out below 1 line 
    176190               !ky 06/09/2016 uncommented below 1 line for FA test! 
    177191               !ky 11/12/2015 recommented out below 1 line 
    178192               !!ky 3/12/2015 uncommented below 1 line for FA test! 
    179                !CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     193               CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     194               !ky 01/11/2016 added 1 line below for test 
     195               CALL iom_put( "rnfmsk_arcmsk", rnfmsk_arcmsk )                         ! river mough & arctic  
    180196            ENDIF 
    181197            ! 
     
    238254      ! 
    239255      !                            !* Allocate erp and qrp array 
     256      !ky 01/11/2016 uncommented out below two lines for test 
    240257      !ky 16/10/2016 commented out below two lines 
    241258      !ky 06/09/2016 uncommented below two lines for FA test! 
    242259      !ky 11/12/2015 recommented out below two lines 
    243260      !!ky 3/12/2015 uncommented below two lines for FA test! 
    244       !ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror ) 
    245       !IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' ) 
     261      ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror ) 
     262      IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' ) 
    246263      ! 
    247264      IF( nn_flxadjht == 1 ) THEN      !* set sf_qrp structure & allocate arrays 
     
    274291      ! 
    275292      !ky!!                            !* Initialize qrp and erp if no restoring  
     293      !ky 01/11/2016 uncommented below two lines for test 
    276294      !ky 16/10/2016 commented out below two lines 
    277295      !ky 06/09/2016 below two lines for FA test! 
    278296      !ky 11/12/2015 commented out below two lines 
    279297      !!ky 3/12/2015 below two lines for FA test! 
    280       !IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp 
    281       !IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp 
     298      IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp 
     299      IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp 
    282300      !!ky!!IF( nn_sstr /= 1                   )   qrp(:,:) = 0._wp 
    283301      !!ky!!IF( nn_sssr /= 1 .OR. nn_sssr /= 2 )   erp(:,:) = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.