Changeset 6517


Ignore:
Timestamp:
2016-05-11T12:12:20+02:00 (6 years ago)
Author:
rfurner
Message:

small fixes to surge flux module

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcsurge.F90

    r6389 r6517  
    122122 
    123123      !                                            ! compute the surface ocean fluxes using CORE bulk formulea 
    124       IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL surge_oce( kt, sf, sst_m, ssu_m, ssv_m, rn_charn_const ) 
     124      IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL surge_oce( kt, sf, ssu_m, ssv_m, rn_charn_const ) 
    125125      ! 
    126126   END SUBROUTINE sbc_surge 
    127127    
    128128    
    129    SUBROUTINE surge_oce( kt, sf, pst, pu, pv, rn_charn_const ) 
     129   SUBROUTINE surge_oce( kt, sf, pu, pv, rn_charn_const ) 
    130130      !!--------------------------------------------------------------------- 
    131131      !!                     ***  ROUTINE surge_oce  *** 
     
    146146      INTEGER  , INTENT(in   )                 ::   kt    ! time step index 
    147147      TYPE(fld), INTENT(inout), DIMENSION(:)   ::   sf    ! input data 
    148       REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pst   ! surface temperature                      [Celcius] 
    149148      REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pu    ! surface current at U-point (i-component) [m/s] 
    150149      REAL(wp) , INTENT(in)   , DIMENSION(:,:) ::   pv    ! surface current at V-point (j-component) [m/s] 
     
    199198            !Iterate 
    200199            DO WHILE((abs(Cd(ji,jj)-z_Cd1))>1E-6) 
    201             z_Cd1=Cd(ji,jj) 
    202             z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav 
    203             Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2 
    204             zi=zi+1 
     200               z_Cd1=Cd(ji,jj) 
     201               z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav 
     202               Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2 
     203               zi=zi+1 
    205204            ENDDO 
    206205         ENDDO 
Note: See TracChangeset for help on using the changeset viewer.