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 4313 – NEMO

Changeset 4313


Ignore:
Timestamp:
2013-11-22T13:12:59+01:00 (10 years ago)
Author:
cetlod
Message:

dev_MERGE_2013 : fix on argument call of eos

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r4230 r4313  
    243243      tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:)  * tmask(:,:,:)    ! salinity 
    244244      ! 
    245       CALL eos    ( tsn, rhd, rhop )                                       ! In any case, we need rhop 
     245      CALL eos    ( tsn, rhd, rhop, gdept_0(:,:,:) )                                       ! In any case, we need rhop 
    246246      CALL zdf_mxl( kt )                                                   ! In any case, we need mxl  
    247247      ! 
     
    521521      !!--------------------------------------------------------------------- 
    522522#if defined key_ldfslp && ! defined key_c1d 
    523       CALL eos( pts, rhd, rhop )   ! Time-filtered in situ density  
     523      CALL eos( pts, rhd, rhop, gdept_0(:,:,:) )   ! Time-filtered in situ density  
    524524      CALL bn2( pts, rn2 )         ! before Brunt-Vaisala frequency 
    525525      IF( ln_zps )   & 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r4245 r4313  
    761761            tsb(:,:,:,:) = tsn(:,:,:,:)               ! Update before fields 
    762762 
    763             CALL eos( tsb, rhd, rhop )                ! Before potential and in situ densities 
     763            CALL eos( tsb, rhd, rhop, gdept_0(:,:,:) )                ! Before potential and in situ densities 
    764764          
    765765            IF( ln_zps .AND. .NOT. lk_c1d ) & 
     
    769769 
    770770#if defined key_zdfkpp 
    771             CALL eos( tsn, rhd )                      ! Compute rhd 
     771            CALL eos( tsn, rhd, fsdept_n(:,:,:) )                      ! Compute rhd 
    772772#endif 
    773773 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/C1D/step_c1d.F90

    r4247 r4313  
    130130                             CALL tra_nxt    ( kstp )        ! tracer fields at next time step 
    131131      IF( ln_zdfnpc      )   CALL tra_npc    ( kstp )        ! applied non penetrative convective adjustment on (t,s) 
    132                              CALL eos( tsb, rhd, rhop )      ! now (swap=before) in situ density for dynhpg module 
     132                             CALL eos( tsb, rhd, rhop, gdept_0(:,:,:) )      ! now (swap=before) in situ density for dynhpg module 
    133133 
    134134      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r4292 r4313  
    9898      ztsn(:,:,:,jp_tem) = tsn(:,:,:,jp_tem)                    ! thermosteric ssh 
    9999      ztsn(:,:,:,jp_sal) = sn0(:,:,:) 
    100       CALL eos( ztsn, zrhd )                       ! now in situ density using initial salinity 
     100      CALL eos( ztsn, zrhd, fsdept_n(:,:,:) )                       ! now in situ density using initial salinity 
    101101      ! 
    102102      zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
     
    112112       
    113113      !                                         ! steric sea surface height 
    114       CALL eos( tsn, zrhd, zrhop )                 ! now in situ and potential density 
     114      CALL eos( tsn, zrhd, zrhop, fsdept_n(:,:,:) )                 ! now in situ and potential density 
    115115      zrhop(:,:,jpk) = 0._wp 
    116116      CALL iom_put( 'rhop', zrhop ) 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r4292 r4313  
    130130         ENDIF 
    131131         ! 
    132          CALL eos( tsb, rhd, rhop )        ! before potential and in situ densities 
     132         CALL eos( tsb, rhd, rhop, gdept_0(:,:,:) )        ! before potential and in situ densities 
    133133#if ! defined key_c1d 
    134134         IF( ln_zps )   CALL zps_hde( nit000, jpts, tsb, gtsu, gtsv,  & ! zps: before hor. gradient 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

    r4292 r4313  
    228228         CALL iom_get( numror, jpdom_autoglo, 'rhop'   , rhop    )   ! now    potential density 
    229229      ELSE 
    230          CALL eos    ( tsn, rhd, rhop )    
     230         CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) )    
    231231      ENDIF 
    232232#if defined key_zdfkpp 
     
    234234         CALL iom_get( numror, jpdom_autoglo, 'rhd'    , rhd     )   ! now    in situ density anomaly 
    235235      ELSE 
    236          CALL eos( tsn, rhd )   ! compute rhd 
     236         CALL eos( tsn, rhd, fsdept_n(:,:,:) )   ! compute rhd 
    237237      ENDIF 
    238238#endif 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/tranpc.F90

    r3294 r4313  
    8181         inpci = 0 
    8282 
    83          CALL eos( tsa, rhd, zrhop )         ! Potential density 
     83         CALL eos( tsa, rhd, zrhop, fsdept_n(:,:,:) )         ! Potential density 
    8484 
    8585         IF( l_trdtra )   THEN                    !* Save ta and sa trends 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRD/trdicp.F90

    r3294 r4313  
    276276         zkepe(:,:,:) = 0._wp 
    277277    
    278          CALL eos( tsn, rhd, rhop )       ! now potential and in situ densities 
     278         CALL eos( tsn, rhd, rhop, fsdept_n(:,:,:) )       ! now potential and in situ densities 
    279279 
    280280         zcof = 0.5_wp / rau0             ! Density flux at w-point 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/step.F90

    r4312 r4313  
    112112          ! Note that the computation of vertical velocity above, hence "after" sea level 
    113113          ! is necessary to compute momentum advection for the rhs of barotropic loop: 
    114                                   CALL eos    ( tsn, rhd, rhop )                 ! now in situ density for hpg computation 
     114                                  CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) )                 ! now in situ density for hpg computation 
    115115          IF( ln_zps      )       CALL zps_hde( kstp, jpts, tsn, gtsu, gtsv,  &  ! zps: now hor. derivative 
    116116                &                                          rhd, gru , grv  )     ! of t, s, rd at the last ocean level 
Note: See TracChangeset for help on using the changeset viewer.