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 7351 for branches/2016/dev_INGV_UKMO_2016/NEMOGCM/TOOLS/SIREN/src/interp_linear.f90 – NEMO

Ignore:
Timestamp:
2016-11-28T17:04:10+01:00 (7 years ago)
Author:
emanuelaclementi
Message:

ticket #1805 step 3: /2016/dev_INGV_UKMO_2016 aligned to the trunk at revision 7161

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/TOOLS/SIREN/src/interp_linear.f90

    r5617 r7351  
    627627       
    628628      IF( ld_even(jp_I) )THEN 
    629          dl_dx=1./REAL(id_rho(jp_I)-1) 
     629         dl_dx=1._dp/REAL(id_rho(jp_I)-1,dp) 
    630630      ELSE ! odd refinement 
    631          dl_dx=1./REAL(id_rho(jp_I)) 
     631         dl_dx=1._dp/REAL(id_rho(jp_I),dp) 
    632632      ENDIF 
    633633 
    634634      IF( ld_even(jp_J) )THEN 
    635          dl_dy=1./REAL(id_rho(jp_J)-1) 
     635         dl_dy=1._dp/REAL(id_rho(jp_J)-1,dp) 
    636636      ELSE ! odd refinement 
    637          dl_dy=1./REAL(id_rho(jp_J)) 
     637         dl_dy=1._dp/REAL(id_rho(jp_J),dp) 
    638638      ENDIF 
    639639 
    … …  
    642642 
    643643         IF( ld_even(jp_J) )THEN 
    644             dl_y=(jj-1)*dl_dy - dl_dy*0.5  
     644            dl_y=REAL(jj-1,dp)*dl_dy - dl_dy*0.5_dp 
    645645         ELSE ! odd refinement 
    646             dl_y=(jj-1)*dl_dy  
     646            dl_y=REAL(jj-1,dp)*dl_dy  
    647647         ENDIF 
    648648 
    … …  
    653653 
    654654            IF( ld_even(jp_I) )THEN 
    655                dl_x=(ji-1)*dl_dx - dl_dx*0.5  
     655               dl_x=REAL(ji-1,dp)*dl_dx - dl_dx*0.5_dp  
    656656            ELSE ! odd refinement 
    657                dl_x=(ji-1)*dl_dx  
     657               dl_x=REAL(ji-1,dp)*dl_dx  
    658658            ENDIF 
    659659 
    … …  
    692692       
    693693      IF( ld_even )THEN 
    694          dl_dx=1./REAL(id_rho-1) 
     694         dl_dx=1._dp/REAL(id_rho-1,dp) 
    695695      ELSE ! odd refinement 
    696          dl_dx=1./REAL(id_rho) 
     696         dl_dx=1._dp/REAL(id_rho,dp) 
    697697      ENDIF 
    698698 
    699699      DO ji=1,id_rho+1 
    700700         IF( ld_even )THEN 
    701             dl_x=(ji-1)*dl_dx - dl_dx*0.5  
     701            dl_x=REAL(ji-1,dp)*dl_dx - dl_dx*0.5_dp  
    702702         ELSE ! odd refinement 
    703             dl_x=(ji-1)*dl_dx  
     703            dl_x=REAL(ji-1,dp)*dl_dx  
    704704         ENDIF 
    705705 
Note: See TracChangeset for help on using the changeset viewer.