Changeset 4708


Ignore:
Timestamp:
2014-07-10T17:17:46+02:00 (6 years ago)
Author:
rfurner
Message:

changes to remove tracer modules and produce barotropic model

Location:
branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r4570 r4708  
    147147 
    148148      IF( lk_vvl ) THEN 
    149          z3d(:,:,:) = tsn(:,:,:,jp_tem) * fse3t_n(:,:,:) 
     149         z3d(:,:,:) = tsn(:,:,:,jp_tem) !* fse3t_n(:,:,:) 
    150150         CALL iom_put( "toce" , z3d                        )   ! heat content 
    151151         CALL iom_put( "sst"  , z3d(:,:,1)                 )   ! sea surface heat content 
    152152         z3d(:,:,1) = tsn(:,:,1,jp_tem) * z3d(:,:,1) 
    153153         CALL iom_put( "sst2" , z3d(:,:,1)                 )   ! sea surface content of squared temperature 
    154          z3d(:,:,:) = tsn(:,:,:,jp_sal) * fse3t_n(:,:,:)             
     154         z3d(:,:,:) = tsn(:,:,:,jp_sal) !* fse3t_n(:,:,:)             
    155155         CALL iom_put( "soce" , z3d                        )   ! salinity content 
    156156         CALL iom_put( "sss"  , z3d(:,:,1)                 )   ! sea surface salinity content 
     
    166166      END IF 
    167167      IF( lk_vvl .AND. (.NOT. ln_dynadv_vec) ) THEN 
    168          CALL iom_put( "uoce" , un(:,:,:) * fse3u_n(:,:,:) )    ! i-transport 
    169          CALL iom_put( "voce" , vn(:,:,:) * fse3v_n(:,:,:) )    ! j-transport 
     168         CALL iom_put( "uoce" , un(:,:,:) )!* fse3u_n(:,:,:) )    ! i-transport 
     169         CALL iom_put( "voce" , vn(:,:,:) )!* fse3v_n(:,:,:) )    ! j-transport 
    170170      ELSE 
    171171         CALL iom_put( "uoce" , un                         )    ! i-current 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r4370 r4708  
    203203      INTEGER  :: ji, jj, jk 
    204204      REAL(wp) ::   zsal = 35.50 
     205#if defined key_surge 
     206      REAL(wp) ::   ztem = 10.0 
     207#endif 
    205208      !!---------------------------------------------------------------------- 
    206209      ! 
     
    209212      IF(lwp) WRITE(numout,*) '~~~~~~~~~~   and constant salinity (',zsal,' psu)' 
    210213      ! 
     214#if defined key_surge 
     215      tsn(:,:,:,jp_tem) = ztem * tmask(:,:,:) 
     216      tsb(:,:,:,jp_tem) = tsn(:,:,:,jp_tem) 
     217#else 
    211218      DO jk = 1, jpk 
    212219         tsn(:,:,jk,jp_tem) = (  ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH((fsdept(:,:,jk)-80.)/30.) )   & 
     
    214221         tsb(:,:,jk,jp_tem) = tsn(:,:,jk,jp_tem) 
    215222      END DO 
     223#endif 
    216224      tsn(:,:,:,jp_sal) = zsal * tmask(:,:,:) 
    217225      tsb(:,:,:,jp_sal) = tsn(:,:,:,jp_sal) 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r4624 r4708  
    303303      ! ocean albedo assumed to be constant + modify now Qsr to include the diurnal cycle                    ! Short Wave 
    304304      zztmp = 1. - albo 
     305#if defined key_surge 
     306      qsr(:,:)=0._wp 
     307      zqlw(:,:) = 0._wp 
     308#else 
    305309      IF( ln_dm2dc ) THEN   ;   qsr(:,:) = zztmp * sbc_dcy( sf(jp_qsr)%fnow(:,:,1) ) * tmask(:,:,1) 
    306310      ELSE                  ;   qsr(:,:) = zztmp *          sf(jp_qsr)%fnow(:,:,1)   * tmask(:,:,1) 
     
    308312!CDIR COLLAPSE 
    309313      zqlw(:,:) = (  sf(jp_qlw)%fnow(:,:,1) - Stef * zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:)  ) * tmask(:,:,1)   ! Long  Wave 
     314#endif 
    310315      ! ----------------------------------------------------------------------------- ! 
    311316      !     II    Turbulent FLUXES                                                    ! 
     
    382387      !  Turbulent fluxes over ocean 
    383388      ! ----------------------------- 
     389#if ! defined key_surge 
    384390      IF( ln_2m .OR. ( ln_bulk2z .AND. rn_zqt /= rn_zu ) ) THEN 
    385391         ! Values of temp. and hum. adjusted to height of wind must be used 
     
    394400!CDIR COLLAPSE 
    395401      zqla (:,:) = Lv * zevap(:,:)                                                              ! Latent Heat 
     402#endif 
    396403 
    397404      IF(ln_ctl) THEN 
     
    410417      ! ----------------------------------------------------------------------------- ! 
    411418      
     419#if defined key_surge 
     420      emp (:,:) = 0._wp 
     421      qns(:,:)  = 0._wp 
     422#else 
    412423!CDIR COLLAPSE 
    413424      emp (:,:) = (  zevap(:,:)                                          &   ! mass flux (evap. - precip.) 
     
    421432         &     + sf(jp_snow)%fnow(:,:,1) * rn_pfac                                &   ! add solid  precip heat content at min(Tair,Tsnow) 
    422433         &     * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic  
     434#endif 
    423435      ! 
    424436      CALL iom_put( "qlw_oce",   zqlw )                 ! output downward longwave heat over the ocean 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r4624 r4708  
    252252         nstop = nstop + 1 
    253253      ENDIF 
     254      IF ( lk_surge .and. .not. ( ln_blk_core .or. ln_ana ) ) & 
     255         &   CALL ctl_stop( ' surge model only compatible with analytical fluxes or core formulae' )      
    254256      IF(lwp) THEN 
    255257         WRITE(numout,*) 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r4645 r4708  
    419419         &                  CALL zdf_ddm_init      ! double diffusive mixing 
    420420      !                                         ! Lateral physics 
     421#if ! defined key_surge 
    421422                            CALL ldf_tra_init      ! Lateral ocean tracer physics 
     423#endif 
    422424                            CALL ldf_dyn_init      ! Lateral ocean momentum physics 
    423425      IF( lk_ldfslp     )   CALL ldf_slp_init      ! slope of lateral mixing 
    424426 
     427#if ! defined key_surge 
    425428      !                                     ! Active tracers 
    426429                            CALL tra_qsr_init   ! penetrative solar radiation qsr 
     
    431434                            CALL tra_ldf_init   ! lateral mixing 
    432435                            CALL tra_zdf_init   ! vertical mixing and after tracer fields 
     436#endif 
    433437 
    434438      !                                     ! Dynamics 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/par_oce.F90

    r4205 r4708  
    9999   LOGICAL, PUBLIC, PARAMETER ::   lk_vopt_loop = .FALSE.  !: vector optimization flag 
    100100#endif 
     101#if defined key_surge 
     102   LOGICAL, PUBLIC, PARAMETER ::   lk_surge     = .TRUE.   ! flag for 2d baratropic modelling 
     103#else 
     104   LOGICAL, PUBLIC, PARAMETER ::   lk_surge     = .FALSE.  ! flag for 2d baratropic modelling 
     105#endif 
    101106 
    102107   !!---------------------------------------------------------------------- 
  • branches/UKMO/2014_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/step.F90

    r4624 r4708  
    223223      ! Active tracers                              (ua, va used as workspace) 
    224224      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     225#if ! defined key_surge 
    225226                             tsa(:,:,:,:) = 0.e0            ! set tracer trends to zero 
    226227 
     
    257258                             CALL tra_nxt( kstp )                ! tracer fields at next time step 
    258259      ENDIF 
    259  
     260#endif 
    260261      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    261262      ! Dynamics                                    (tsa used as workspace) 
Note: See TracChangeset for help on using the changeset viewer.