New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3432 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 – NEMO

Ignore:
Timestamp:
2012-07-11T13:22:58+02:00 (12 years ago)
Author:
trackstand2
Message:

Merge branch 'ksection_partition'

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r3211 r3432  
    185185      CALL fld_read( kt, nn_fsbc, sf )        ! input fields provided at the current time-step 
    186186 
    187 #if defined key_lim3 
    188       tatm_ice(:,:) = sf(jp_tair)%fnow(:,:,1)                  ! LIM3: make Tair available in sea-ice 
    189 #endif 
    190187      !                                                        ! surface ocean fluxes computed with CLIO bulk formulea 
    191188      IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL blk_oce_core( sf, sst_m, ssu_m, ssv_m ) 
    … …  
    495492!!gm end 
    496493 
     494#if defined key_lim3 
     495      tatm_ice(:,:) = sf(jp_tair)%fnow(:,:,1)   ! LIM3: make Tair available in sea-ice. WARNING allocated after call to ice_init 
     496#endif 
    497497      ! ----------------------------------------------------------------------------- ! 
    498498      !    Wind components and module relative to the moving ocean ( U10m - U_ice )   ! 
    … …  
    790790      !!---------------------------------------------------------------------- 
    791791      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released, iwrk_in_use, iwrk_not_released 
    792       USE wrk_nemo, ONLY: dU10 => wrk_2d_1        ! dU                             [m/s] 
    793       USE wrk_nemo, ONLY: dT => wrk_2d_2          ! air/sea temperature difference   [K] 
    794       USE wrk_nemo, ONLY: dq => wrk_2d_3          ! air/sea humidity difference      [K] 
    795       USE wrk_nemo, ONLY: Cd_n10 => wrk_2d_4      ! 10m neutral drag coefficient 
    796       USE wrk_nemo, ONLY: Ce_n10 => wrk_2d_5      ! 10m neutral latent coefficient 
    797       USE wrk_nemo, ONLY: Ch_n10 => wrk_2d_6      ! 10m neutral sensible coefficient 
    798       USE wrk_nemo, ONLY: sqrt_Cd_n10 => wrk_2d_7 ! root square of Cd_n10 
    799       USE wrk_nemo, ONLY: sqrt_Cd => wrk_2d_8     ! root square of Cd 
    800       USE wrk_nemo, ONLY: T_vpot => wrk_2d_9      ! virtual potential temperature    [K] 
    801       USE wrk_nemo, ONLY: T_star => wrk_2d_10     ! turbulent scale of tem. fluct. 
    802       USE wrk_nemo, ONLY: q_star => wrk_2d_11     ! turbulent humidity of temp. fluct. 
    803       USE wrk_nemo, ONLY: U_star => wrk_2d_12     ! turb. scale of velocity fluct. 
    804       USE wrk_nemo, ONLY: L => wrk_2d_13          ! Monin-Obukov length              [m] 
    805       USE wrk_nemo, ONLY: zeta_u => wrk_2d_14     ! stability parameter at height zu 
    806       USE wrk_nemo, ONLY: zeta_t => wrk_2d_15     ! stability parameter at height zt 
    807       USE wrk_nemo, ONLY: U_n10 => wrk_2d_16      ! neutral wind velocity at 10m     [m] 
    808       USE wrk_nemo, ONLY: xlogt => wrk_2d_17, xct => wrk_2d_18, zpsi_hu => wrk_2d_19, zpsi_ht => wrk_2d_20, zpsi_m => wrk_2d_21 
     792      USE wrk_nemo, ONLY: dU10 => wrk_2d_14        ! dU                             [m/s] 
     793      USE wrk_nemo, ONLY: dT => wrk_2d_15          ! air/sea temperature difference   [K] 
     794      USE wrk_nemo, ONLY: dq => wrk_2d_16          ! air/sea humidity difference      [K] 
     795      USE wrk_nemo, ONLY: Cd_n10 => wrk_2d_17      ! 10m neutral drag coefficient 
     796      USE wrk_nemo, ONLY: Ce_n10 => wrk_2d_18      ! 10m neutral latent coefficient 
     797      USE wrk_nemo, ONLY: Ch_n10 => wrk_2d_19      ! 10m neutral sensible coefficient 
     798      USE wrk_nemo, ONLY: sqrt_Cd_n10 => wrk_2d_20 ! root square of Cd_n10 
     799      USE wrk_nemo, ONLY: sqrt_Cd => wrk_2d_21     ! root square of Cd 
     800      USE wrk_nemo, ONLY: T_vpot => wrk_2d_22      ! virtual potential temperature    [K] 
     801      USE wrk_nemo, ONLY: T_star => wrk_2d_23     ! turbulent scale of tem. fluct. 
     802      USE wrk_nemo, ONLY: q_star => wrk_2d_24     ! turbulent humidity of temp. fluct. 
     803      USE wrk_nemo, ONLY: U_star => wrk_2d_25     ! turb. scale of velocity fluct. 
     804      USE wrk_nemo, ONLY: L => wrk_2d_26          ! Monin-Obukov length              [m] 
     805      USE wrk_nemo, ONLY: zeta_u => wrk_2d_27     ! stability parameter at height zu 
     806      USE wrk_nemo, ONLY: zeta_t => wrk_2d_28     ! stability parameter at height zt 
     807      USE wrk_nemo, ONLY: U_n10 => wrk_2d_29      ! neutral wind velocity at 10m     [m] 
     808      USE wrk_nemo, ONLY: xlogt => wrk_2d_30, xct => wrk_2d_31, zpsi_hu => wrk_2d_32, zpsi_ht => wrk_2d_33, zpsi_m => wrk_2d_34 
    809809      USE wrk_nemo, ONLY: stab => iwrk_2d_1      ! 1st guess stability test integer 
    810810      !! 
    … …  
    833833      !!  * Start 
    834834 
    835       IF(  wrk_in_use(2, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21) .OR. & 
    836           iwrk_in_use(2, 1) ) THEN 
     835      IF(  wrk_in_use(2,             14,15,16,17,18,19,        & 
     836                         20,21,22,23,24,25,26,27,28,29,        &          
     837                         30,31,32,33,34)                .OR.   & 
     838          iwrk_in_use(2, 1)                               ) THEN 
    837839         CALL ctl_stop('TURB_CORE_2Z: requested workspace arrays unavailable')   ;   RETURN 
    838       END IF 
     840      ENDIF 
    839841 
    840842      !! Initial air/sea differences 
    … …  
    911913      END DO 
    912914      !! 
    913       IF(  wrk_not_released(2, 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21) .OR.   & 
    914           iwrk_not_released(2, 1)    )   CALL ctl_stop('TURB_CORE_2Z: requested workspace arrays unavailable') 
     915     IF( wrk_not_released(2,              14,15,16,17,18,19,          & 
     916         &                    20,21,22,23,24,25,26,27,28,29,          & 
     917         &                    30,31,32,33,34                )   .OR.  &   
     918         iwrk_not_released(2, 1)                                  )   & 
     919         CALL ctl_stop('TURB_CORE_2Z: failed to release workspace arrays') 
    915920      ! 
    916921    END SUBROUTINE TURB_CORE_2Z 
    … …  
    920925      !------------------------------------------------------------------------------- 
    921926      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 
    922       USE wrk_nemo, ONLY:     X2 => wrk_2d_33 
    923       USE wrk_nemo, ONLY:     X  => wrk_2d_34 
    924       USE wrk_nemo, ONLY: stabit => wrk_2d_35 
     927      USE wrk_nemo, ONLY:     X2 => wrk_2d_35 
     928      USE wrk_nemo, ONLY:     X  => wrk_2d_36 
     929      USE wrk_nemo, ONLY: stabit => wrk_2d_37 
    925930      !! 
    926931      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: zta 
    … …  
    930935      !------------------------------------------------------------------------------- 
    931936 
    932       IF( wrk_in_use(2, 33,34,35) ) THEN 
     937      IF( wrk_in_use(2, 35,36,37) ) THEN 
    933938         CALL ctl_stop('psi_m: requested workspace arrays unavailable')   ;   RETURN 
    934939      ENDIF 
    … …  
    939944         &    + (1. - stabit)*(2*log((1. + X)/2) + log((1. + X2)/2) - 2*atan(X) + pi/2)  ! Unstable  
    940945 
    941       IF( wrk_not_released(2, 33,34,35) )   CALL ctl_stop('psi_m: failed to release workspace arrays') 
     946      IF( wrk_not_released(2, 35,36,37) )   CALL ctl_stop('psi_m: failed to release workspace arrays') 
    942947      ! 
    943948    END FUNCTION psi_m 
    … …  
    947952      !------------------------------------------------------------------------------- 
    948953      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 
    949       USE wrk_nemo, ONLY:     X2 => wrk_2d_33 
    950       USE wrk_nemo, ONLY:     X  => wrk_2d_34 
    951       USE wrk_nemo, ONLY: stabit => wrk_2d_35 
     954      USE wrk_nemo, ONLY:     X2 => wrk_2d_35 
     955      USE wrk_nemo, ONLY:     X  => wrk_2d_36 
     956      USE wrk_nemo, ONLY: stabit => wrk_2d_37 
    952957      ! 
    953958      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   zta 
    … …  
    956961      !------------------------------------------------------------------------------- 
    957962 
    958       IF( wrk_in_use(2, 33,34,35) ) THEN 
     963      IF( wrk_in_use(2, 35,36,37) ) THEN 
    959964         CALL ctl_stop('psi_h: requested workspace arrays unavailable')   ;   RETURN 
    960965      ENDIF 
    … …  
    965970         &    + (1. - stabit)*(2.*log( (1. + X2)/2. ))                 ! Unstable 
    966971 
    967       IF( wrk_not_released(2, 33,34,35) )   CALL ctl_stop('psi_h: failed to release workspace arrays') 
     972      IF( wrk_not_released(2, 35,36,37) )   CALL ctl_stop('psi_h: failed to release workspace arrays') 
    968973      ! 
    969974    END FUNCTION psi_h 
Note: See TracChangeset for help on using the changeset viewer.