Changeset 8408


Ignore:
Timestamp:
2017-08-07T16:52:13+02:00 (3 years ago)
Author:
jcastill
Message:

Remove the code that does not belong to this branch, and that is already included in the amm15 configuration branch AMM15_v3_6_STABLE_package

Location:
branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

    r7454 r8408  
    2828   PUBLIC sbc_flx       ! routine called by step.F90 
    2929 
    30    INTEGER , PARAMETER ::   jpfld   = 6   ! maximum number of files to read  
     30   INTEGER , PARAMETER ::   jpfld   = 5   ! maximum number of files to read  
    3131   INTEGER , PARAMETER ::   jp_utau = 1   ! index of wind stress (i-component) file 
    3232   INTEGER , PARAMETER ::   jp_vtau = 2   ! index of wind stress (j-component) file 
     
    3434   INTEGER , PARAMETER ::   jp_qsr  = 4   ! index of solar heat file 
    3535   INTEGER , PARAMETER ::   jp_emp  = 5   ! index of evaporation-precipation file 
    36    INTEGER , PARAMETER ::   jp_press = 6  ! index of pressure for UKMO shelf fluxes 
    3736   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf    ! structure of input fields (file informations, fields read) 
    38    LOGICAL , PUBLIC    ::   ln_shelf_flx = .FALSE. ! UKMO SHELF specific flux flag 
    39    INTEGER             ::   jpfld_local   ! maximum number of files to read (locally modified depending on ln_shelf_flx)  
    4037 
    4138   !! * Substitutions 
     
    8582      REAL(wp) ::   zcdrag = 1.5e-3       ! drag coefficient 
    8683      REAL(wp) ::   ztx, zty, zmod, zcoef ! temporary variables 
    87       REAL     ::   cs                    ! UKMO SHELF: Friction co-efficient at surface 
    88       REAL     ::   totwindspd            ! UKMO SHELF: Magnitude of wind speed vector 
    89  
    90       REAL(wp) ::   rhoa  = 1.22         ! Air density kg/m3 
    91       REAL(wp) ::   cdrag = 1.5e-3       ! drag coefficient  
    9284      !! 
    9385      CHARACTER(len=100) ::  cn_dir                               ! Root directory for location of flx files 
    9486      TYPE(FLD_N), DIMENSION(jpfld) ::   slf_i                    ! array of namelist information structures 
    95       TYPE(FLD_N) ::   sn_utau, sn_vtau, sn_qtot, sn_qsr, sn_emp, sn_press  !  informations about the fields to be read 
    96       LOGICAL     ::   ln_foam_flx  = .FALSE.                     ! UKMO FOAM specific flux flag 
    97       NAMELIST/namsbc_flx/ cn_dir, sn_utau, sn_vtau, sn_qtot, sn_qsr, sn_emp,   & 
    98       &                    ln_foam_flx, sn_press, ln_shelf_flx 
     87      TYPE(FLD_N) ::   sn_utau, sn_vtau, sn_qtot, sn_qsr, sn_emp  !  informations about the fields to be read 
     88      NAMELIST/namsbc_flx/ cn_dir, sn_utau, sn_vtau, sn_qtot, sn_qsr, sn_emp 
    9989      !!--------------------------------------------------------------------- 
    10090      ! 
     
    119109         slf_i(jp_emp ) = sn_emp 
    120110         ! 
    121             ALLOCATE( sf(jpfld), STAT=ierror )        ! set sf structure 
    122             IF( ln_shelf_flx ) slf_i(jp_press) = sn_press 
    123     
    124             ! define local jpfld depending on shelf_flx logical 
    125             IF( ln_shelf_flx ) THEN 
    126                jpfld_local = jpfld 
    127             ELSE 
    128                jpfld_local = jpfld-1 
    129             ENDIF 
    130             ! 
     111         ALLOCATE( sf(jpfld), STAT=ierror )        ! set sf structure 
    131112         IF( ierror > 0 ) THEN    
    132113            CALL ctl_stop( 'sbc_flx: unable to allocate sf structure' )   ;   RETURN   
    133114         ENDIF 
    134          DO ji= 1, jpfld_local 
     115         DO ji= 1, jpfld 
    135116            ALLOCATE( sf(ji)%fnow(jpi,jpj,1) ) 
    136117            IF( slf_i(ji)%ln_tint ) ALLOCATE( sf(ji)%fdta(jpi,jpj,1,2) ) 
     
    151132         ENDIF 
    152133!CDIR COLLAPSE 
    153             !!UKMO SHELF effect of atmospheric pressure on SSH 
    154             ! If using ln_apr_dyn, this is done there so don't repeat here. 
    155             IF( ln_shelf_flx .AND. .NOT. ln_apr_dyn) THEN 
    156                DO jj = 1, jpjm1 
    157                   DO ji = 1, jpim1 
    158                      apgu(ji,jj) = (-1.0/rau0)*(sf(jp_press)%fnow(ji+1,jj,1)-sf(jp_press)%fnow(ji,jj,1))/e1u(ji,jj) 
    159                      apgv(ji,jj) = (-1.0/rau0)*(sf(jp_press)%fnow(ji,jj+1,1)-sf(jp_press)%fnow(ji,jj,1))/e2v(ji,jj) 
    160                   END DO 
    161                END DO 
    162             ENDIF ! ln_shelf_flx 
    163        
    164134         DO jj = 1, jpj                                           ! set the ocean fluxes from read fields 
    165135            DO ji = 1, jpi 
    166                    IF( ln_shelf_flx ) THEN 
    167                       !! UKMO SHELF - need atmospheric pressure to calculate Haney forcing 
    168                       pressnow(ji,jj) = sf(jp_press)%fnow(ji,jj,1) 
    169                       !! UKMO SHELF flux files contain wind speed not wind stress 
    170                       totwindspd = sqrt((sf(jp_utau)%fnow(ji,jj,1))**2.0 + (sf(jp_vtau)%fnow(ji,jj,1))**2.0) 
    171                       cs = 0.63 + (0.066 * totwindspd) 
    172                       utau(ji,jj) = cs * (rhoa/rau0) * sf(jp_utau)%fnow(ji,jj,1) * totwindspd 
    173                       vtau(ji,jj) = cs * (rhoa/rau0) * sf(jp_vtau)%fnow(ji,jj,1) * totwindspd 
    174                    ELSE 
    175                       utau(ji,jj) = sf(jp_utau)%fnow(ji,jj,1) 
    176                       vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1) 
    177                    ENDIF 
    178                    qsr (ji,jj) = sf(jp_qsr )%fnow(ji,jj,1) 
    179                    IF( ln_foam_flx .OR. ln_shelf_flx ) THEN 
    180                       !! UKMO FOAM flux files contain non-solar heat flux (qns) rather than total heat flux (qtot) 
    181                       qns (ji,jj) = sf(jp_qtot)%fnow(ji,jj,1) 
    182                       !! UKMO FOAM flux files contain the net DOWNWARD freshwater flux P-E rather then E-P 
    183                       emp (ji,jj) = -1. * sf(jp_emp )%fnow(ji,jj,1) 
    184                    ELSE 
    185                       qns (ji,jj) = sf(jp_qtot)%fnow(ji,jj,1) - sf(jp_qsr)%fnow(ji,jj,1) 
    186                       emp (ji,jj) = sf(jp_emp )%fnow(ji,jj,1) 
    187                    ENDIF 
     136               utau(ji,jj) = sf(jp_utau)%fnow(ji,jj,1) 
     137               vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1) 
     138               qns (ji,jj) = sf(jp_qtot)%fnow(ji,jj,1) - sf(jp_qsr)%fnow(ji,jj,1) 
     139               emp (ji,jj) = sf(jp_emp )%fnow(ji,jj,1) 
    188140            END DO 
    189141         END DO 
     
    191143         qns(:,:) = qns(:,:) - emp(:,:) * sst_m(:,:) * rcp        ! mass flux is at SST 
    192144         ! 
    193     
    194             !! UKMO FOAM wind fluxes need lbc_lnk calls owing to a bug in interp.exe 
    195             IF( ln_foam_flx ) THEN 
    196                CALL lbc_lnk( utau(:,:), 'U', -1. ) 
    197                CALL lbc_lnk( vtau(:,:), 'V', -1. ) 
    198             ENDIF 
    199      
    200145         !                                                        ! module of wind stress and wind speed at T-point 
    201146         zcoef = 1. / ( zrhoa * zcdrag ) 
     
    217162            WRITE(numout,*)  
    218163            WRITE(numout,*) '        read daily momentum, heat and freshwater fluxes OK' 
    219             DO jf = 1, jpfld_local 
     164            DO jf = 1, jpfld 
    220165               IF( jf == jp_utau .OR. jf == jp_vtau )   zfact =     1. 
    221166               IF( jf == jp_qtot .OR. jf == jp_qsr  )   zfact =     0.1 
  • branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

    r7454 r8408  
    4242   LOGICAL         ::   ln_sssr_bnd     ! flag to bound erp term  
    4343   REAL(wp)        ::   rn_sssr_bnd     ! ABS(Max./Min.) value of erp term [mm/day] 
    44    LOGICAL         ::   ln_UKMO_haney   ! UKMO specific flag to calculate Haney forcing   
    4544 
    4645   REAL(wp) , ALLOCATABLE, DIMENSION(:) ::   buffer   ! Temporary buffer for exchange 
     
    8079      INTEGER  ::   ierror   ! return error code 
    8180      !! 
    82       REAL(wp) ::   sst1,sst2                      ! sea surface temperature 
    83       REAL(wp) ::   e_sst1, e_sst2                 ! saturation vapour pressure 
    84       REAL(wp) ::   qs1,qs2                        ! specific humidity 
    85       REAL(wp) ::   pr_tmp                         ! temporary variable for pressure 
    86   
    87       REAL(wp), DIMENSION(jpi,jpj) ::  hny_frc1    ! Haney forcing for sensible heat, correction for latent heat    
    88       REAL(wp), DIMENSION(jpi,jpj) ::  hny_frc2    ! Haney forcing for sensible heat, correction for latent heat    
    89       !! 
    9081      CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files 
    9182      TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read 
     
    10495            ! 
    10596            IF( nn_sstr == 1 ) THEN                                   !* Temperature restoring term 
    106                   IF( ln_UKMO_haney ) THEN 
    107                      DO jj = 1, jpj 
    108                         DO ji = 1, jpi 
    109                            sst1   =  sst_m(ji,jj) 
    110                            sst2   =  sf_sst(1)%fnow(ji,jj,1)    
    111                            e_sst1 = 10**((0.7859+0.03477*sst1)/(1.+0.00412*sst1)) 
    112                            e_sst2 = 10**((0.7859+0.03477*sst2)/(1.+0.00412*sst2))          
    113                            pr_tmp = 0.01*pressnow(ji,jj)  !pr_tmp = 1012.0 
    114                            qs1    = (0.62197*e_sst1)/(pr_tmp-0.378*e_sst1) 
    115                            qs2    = (0.62197*e_sst2)/(pr_tmp-0.378*e_sst2) 
    116                            hny_frc1(ji,jj) = sst1-sst2                    
    117                            hny_frc2(ji,jj) = qs1-qs2                      
    118                           !Might need to mask off land points. 
    119                            hny_frc1(ji,jj)=-hny_frc1(ji,jj)*wndm(ji,jj)*1.42 
    120                            hny_frc2(ji,jj)=-hny_frc2(ji,jj)*wndm(ji,jj)*4688.0 
    121                            qns(ji,jj)=qns(ji,jj)+hny_frc1(ji,jj)+hny_frc2(ji,jj)    
    122                            qrp(ji,jj) = 0.e0 
    123                         END DO 
    124                      END DO 
    125                   ELSE 
    126                      DO jj = 1, jpj 
    127                         DO ji = 1, jpi 
    128                            zqrp = rn_dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) 
    129                            qns(ji,jj) = qns(ji,jj) + zqrp 
    130                            qrp(ji,jj) = zqrp 
    131                         END DO 
    132                      END DO 
    133                   ENDIF 
     97               DO jj = 1, jpj 
     98                  DO ji = 1, jpi 
     99                     zqrp = rn_dqdt * ( sst_m(ji,jj) - sf_sst(1)%fnow(ji,jj,1) ) 
     100                     qns(ji,jj) = qns(ji,jj) + zqrp 
     101                     qrp(ji,jj) = zqrp 
     102                  END DO 
     103               END DO 
    134104               CALL iom_put( "qrp", qrp )                             ! heat flux damping 
    135105            ENDIF 
     
    193163      CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files 
    194164      TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read 
    195       NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd, ln_UKMO_haney 
     165      NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd 
    196166      INTEGER     ::  ios 
    197167      !!---------------------------------------------------------------------- 
     
    219189         WRITE(numout,*) '      flag to bound erp term                 ln_sssr_bnd = ', ln_sssr_bnd 
    220190         WRITE(numout,*) '      ABS(Max./Min.) erp threshold           rn_sssr_bnd = ', rn_sssr_bnd, ' mm/day' 
    221          WRITE(numout,*) '      Haney forcing                          ln_UKMO_haney = ', ln_UKMO_haney 
    222191      ENDIF 
    223192      ! 
  • branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r7452 r8408  
    2525   USE trd_oce         ! trends: ocean variables 
    2626   USE trdtra          ! trends manager: tracers  
    27    USE tradwl          ! solar radiation penetration (downwell method) 
    2827   ! 
    2928   USE in_out_manager  ! I/O manager 
     
    139138 
    140139!!gm      IF( .NOT.ln_traqsr )   qsr(:,:) = 0.e0   ! no solar radiation penetration 
    141       IF( .NOT.ln_traqsr .and. .NOT.ln_tradwl ) THEN     ! no solar radiation penetration 
     140      IF( .NOT.ln_traqsr ) THEN     ! no solar radiation penetration 
    142141         qns(:,:) = qns(:,:) + qsr(:,:)      ! total heat flux in qns 
    143142         qsr(:,:) = 0.e0                     ! qsr set to zero 
  • branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r7454 r8408  
    2525   USE sbcrnf           ! surface boundary condition: runoff variables 
    2626   USE sbccpl           ! surface boundary condition: coupled formulation (call send at end of step) 
    27    USE sbcflx           ! surface boundary condition: Fluxes 
    2827   USE sbc_oce          ! surface boundary condition: ocean 
    2928   USE sbctide          ! Tide initialisation 
     
    3130 
    3231   USE traqsr           ! solar radiation penetration      (tra_qsr routine) 
    33    USE tradwl           ! POLCOMS style solar radiation    (tra_dwl routine)  
    3432   USE trasbc           ! surface boundary condition       (tra_sbc routine) 
    3533   USE trabbc           ! bottom boundary condition        (tra_bbc routine) 
  • branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

    r7454 r8408  
    2727   REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) ::   etot3         !: light absortion coefficient 
    2828   REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:,:) ::   facvol        !: volume for degraded regions 
    29    REAL(wp), PUBLIC, SAVE, ALLOCATABLE, DIMENSION(:,:)   ::   rlambda2      !: Lambda2 for downwell version of Short wave Radiation 
    30    REAL(wp), PUBLIC                                      ::   rlambda       !: Lambda  for downwell version of Short wave Radiation 
    3129 
    3230#if defined key_top  
     
    8078      !!                  ***  trc_oce_alloc  *** 
    8179      !!---------------------------------------------------------------------- 
    82       INTEGER ::   ierr(3)        ! Local variables 
     80      INTEGER ::   ierr(2)        ! Local variables 
    8381      !!---------------------------------------------------------------------- 
    8482      ierr(:) = 0 
    8583                     ALLOCATE( etot3 (jpi,jpj,jpk), STAT=ierr(1) ) 
    8684      IF( lk_degrad) ALLOCATE( facvol(jpi,jpj,jpk), STAT=ierr(2) ) 
    87                     ALLOCATE( rlambda2(jpi,jpj),   STAT=ierr(3) ) 
    8885      trc_oce_alloc  = MAXVAL( ierr ) 
    8986      ! 
    9087      IF( trc_oce_alloc /= 0 )   CALL ctl_warn('trc_oce_alloc: failed to allocate etot3 array') 
    91       IF( trc_oce_alloc /= 0 )   CALL ctl_warn('trc_oce_alloc: failed to allocate etot3, facvol or rlambda2 array') 
    9288   END FUNCTION trc_oce_alloc 
    9389 
Note: See TracChangeset for help on using the changeset viewer.