Changeset 7031


Ignore:
Timestamp:
2016-10-14T19:06:24+02:00 (4 years ago)
Author:
kuniko
Message:

test mode

File:
1 edited

Legend:

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

    r6913 r7031  
    3232   PUBLIC   sbc_flx_adj_init   ! routine called in sbcmod 
    3333 
     34   !ky 06/09/2016 uncommented below two lines for FA test 
    3435   !ky 11/12/2015 recommented below two lines 
    3536   !!ky 3/12/2015 uncommented below two lines for FA test! 
    36    !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s] 
    37    !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2] 
     37   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   erp   !: evaporation damping   [kg/m2/s] 
     38   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   qrp   !: heat flux damping        [w/m2] 
    3839 
    3940   !                                   !!* Namelist namsbc_flx_adj * 
     
    102103                        &        * sf_qrp(1)%fnow(ji,jj,1) 
    103104                     !zqrp = sf_qrp(1)%fnow(ji,jj,1) 
     105                     !ky 07/09/2016 copied 1 line below for FA test 
    104106                     !ky 11/12/2015 commented out 1 line below 
    105107                     !!ky 3/12/2015 1 line below for FA test! 
    106108                     !qrp(ji,jj) = sf_qrp(1)%fnow(ji,jj,1) 
     109                     !ky 07/09/2016 copied 1 line above and modified as below for FA test 
     110                     qrp(ji,jj) = zqrp 
    107111                     qns(ji,jj) = qns(ji,jj) + zqrp 
    108112                  END DO 
    109113               END DO 
     114               !ky 06/09/2016 uncommented below 1 line below for FA test! 
    110115               !ky 11/12/2015 recommented out below 1 line below 
    111116               !!ky 3/12/2015 uncommented below 1 line below for FA test! 
    112                !CALL iom_put( "qrp", qrp )                             ! heat flux damping 
     117               CALL iom_put( "qrp", qrp )                             ! heat flux damping 
    113118            ENDIF 
    114119            ! 
     
    127132                     !!ky 3/12/2015 1 line below for FA test! 
    128133                     !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 
     134                     !ky 06/09/2016 copied line above and changed to below two lines for FA test! 
     135                     erp(ji,jj) = ( 1. - 2.*rnfmsk_arcmsk(ji,jj) )  &        ! No damping in vicinity of river mouths 
     136                        &        * sf_erp(1)%fnow(ji,jj,1)  
    129137                     sfx(ji,jj) = sfx(ji,jj) + zerp                 ! salt flux 
    130138                  END DO 
    131139               END DO 
     140               !ky 06/09/2016 uncommented below one line for FA test! 
    132141               !ky 11/12/2015 recommented out below one line 
    133142               !!ky 3/12/2015 uncommented below one line for FA test! 
    134                !CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     143               CALL iom_put( "erp", erp )                             ! freshwater flux damping 
    135144               ! 
    136145            ELSEIF( nn_flxadjfw == 2 ) THEN                               !* Salinity damping term (volume flux (emp) and associated heat flux (qns) 
     
    145154                     !!ky 3/12/2015 1 line below for FA test! 
    146155                     !erp(ji,jj) = sf_erp(1)%fnow(ji,jj,1) 
     156                     !ky 06/09/2016 copied line above and changed to below line for FA test! 
     157                     erp(ji,jj) = zerp 
    147158                     emp(ji,jj) = emp (ji,jj) + zerp 
    148159                     qns(ji,jj) = qns(ji,jj) - zerp * rcp * sst_m(ji,jj) 
    149160                  END DO 
    150161               END DO 
     162               !ky 06/09/2016 uncommented below 1 line for FA test! 
    151163               !ky 11/12/2015 recommented out below 1 line 
    152164               !!ky 3/12/2015 uncommented below 1 line for FA test! 
    153                !CALL iom_put( "erp", erp )                             ! freshwater flux damping 
     165               CALL iom_put( "erp", erp )                             ! freshwater flux damping 
    154166            ENDIF 
    155167            ! 
     
    210222      ! 
    211223      !                            !* Allocate erp and qrp array 
     224      !ky 06/09/2016 uncommented below two lines for FA test! 
    212225      !ky 11/12/2015 recommented out below two lines 
    213226      !!ky 3/12/2015 uncommented below two lines for FA test! 
    214       !ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror ) 
    215       !IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' ) 
     227      ALLOCATE( qrp(jpi,jpj), erp(jpi,jpj), STAT=ierror ) 
     228      IF( ierror > 0 )   CALL ctl_stop( 'STOP', 'sbc_ssr: unable to allocate erp and qrp array' ) 
    216229      ! 
    217230      IF( nn_flxadjht == 1 ) THEN      !* set sf_qrp structure & allocate arrays 
     
    244257      ! 
    245258      !ky!!                            !* Initialize qrp and erp if no restoring  
     259      !ky 06/09/2016 below two lines for FA test! 
    246260      !ky 11/12/2015 commented out below two lines 
    247261      !!ky 3/12/2015 below two lines for FA test! 
    248       !IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp 
    249       !IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp 
     262      IF( nn_flxadjht /= 1                   )   qrp(:,:) = 0._wp 
     263      IF( nn_flxadjfw /= 1 .OR. nn_flxadjfw /= 2 )   erp(:,:) = 0._wp 
    250264      !!ky!!IF( nn_sstr /= 1                   )   qrp(:,:) = 0._wp 
    251265      !!ky!!IF( nn_sssr /= 1 .OR. nn_sssr /= 2 )   erp(:,:) = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.