Changeset 11325


Ignore:
Timestamp:
2019-07-23T10:19:17+02:00 (13 months ago)
Author:
smasson
Message:

dev_r10984_HPC-13 : bugfixes for sette with debug options

Files:
5 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdydta.F90

    r11317 r11325  
    415415               bf_alias => bf(jp_bdyssh,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    416416               bn_alias => bn_ssh                                          ! alias for ssh structure of nambdy_dta  
     417               iszdim = idx_bdy(jbdy)%nblenrim(igrd)                       ! length of this bdy on this MPI processus : only rim 
    417418            ENDIF 
    418419            ! ===================== 
     
    427428               bf_alias => bf(jp_bdyu2d,jbdy:jbdy)                         ! alias for u2d structure of bdy number jbdy 
    428429               bn_alias => bn_u2d                                          ! alias for u2d structure of nambdy_dta  
     430               iszdim = idx_bdy(jbdy)%nblenrim(igrd)                       ! length of this bdy on this MPI processus : only rim 
    429431            ENDIF 
    430432            IF( jfld == jp_bdyv2d ) THEN 
     
    436438               bf_alias => bf(jp_bdyv2d,jbdy:jbdy)                         ! alias for v2d structure of bdy number jbdy 
    437439               bn_alias => bn_v2d                                          ! alias for v2d structure of nambdy_dta  
     440               iszdim = idx_bdy(jbdy)%nblenrim(igrd)                       ! length of this bdy on this MPI processus : only rim 
    438441            ENDIF 
    439442            ! ===================== 
     
    449452               bf_alias => bf(jp_bdyu3d,jbdy:jbdy)                         ! alias for u3d structure of bdy number jbdy 
    450453               bn_alias => bn_u3d                                          ! alias for u3d structure of nambdy_dta  
    451             ENDIF 
     454               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
     455           ENDIF 
    452456            IF( jfld == jp_bdyv3d ) THEN 
    453457               cl3 = 'v3d' 
     
    459463               bf_alias => bf(jp_bdyv3d,jbdy:jbdy)                         ! alias for v3d structure of bdy number jbdy 
    460464               bn_alias => bn_v3d                                          ! alias for v3d structure of nambdy_dta  
    461             ENDIF 
     465               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
     466           ENDIF 
    462467 
    463468            ! ===================== 
     
    472477               bf_alias => bf(jp_bdytem,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    473478               bn_alias => bn_tem                                          ! alias for ssh structure of nambdy_dta  
     479               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
    474480            ENDIF 
    475481            IF( jfld == jp_bdysal ) THEN 
     
    481487               bf_alias => bf(jp_bdysal,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    482488               bn_alias => bn_sal                                          ! alias for ssh structure of nambdy_dta  
     489               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
    483490            ENDIF 
    484491 
     
    494501               bf_alias => bf(jp_bdya_i,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    495502               bn_alias => bn_a_i                                          ! alias for ssh structure of nambdy_dta  
    496             ENDIF 
     503               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
     504           ENDIF 
    497505            IF( jfld == jp_bdyh_i ) THEN 
    498506               cl3 = 'h_i' 
     
    503511               bf_alias => bf(jp_bdyh_i,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    504512               bn_alias => bn_h_i                                          ! alias for ssh structure of nambdy_dta  
     513               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
    505514            ENDIF 
    506515            IF( jfld == jp_bdyh_s ) THEN 
     
    512521               bf_alias => bf(jp_bdyh_s,jbdy:jbdy)                         ! alias for ssh structure of bdy number jbdy 
    513522               bn_alias => bn_h_s                                          ! alias for ssh structure of nambdy_dta  
    514             ENDIF 
    515  
     523               iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
     524            ENDIF 
    516525 
    517526            IF( llneed ) THEN                                              ! dta_bdy(jbdy)%xxx will be needed 
    518527               !                                                           !   -> must be associated with an allocated target 
    519                iszdim = idx_bdy(jbdy)%nblen(igrd)                          ! length of this bdy on this MPI processus 
    520528               ALLOCATE( bf_alias(1)%fnow( iszdim, 1, ipk ) )              ! allocate the target 
    521529               ! 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DIA/diawri.F90

    r10425 r11325  
    657657         ENDIF 
    658658 
    659          IF( .NOT. ln_cpl ) THEN 
     659         IF( ln_ssr ) THEN 
    660660            CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping"         , "W/m2"   ,   &  ! qrp 
    661661               &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
     
    665665               &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    666666         ENDIF 
    667  
    668          IF( ln_cpl .AND. nn_ice <= 1 ) THEN 
    669             CALL histdef( nid_T, "sohefldp", "Surface Heat Flux: Damping"         , "W/m2"   ,   &  ! qrp 
    670                &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    671             CALL histdef( nid_T, "sowafldp", "Surface Water Flux: Damping"        , "Kg/m2/s",   &  ! erp 
    672                &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    673             CALL histdef( nid_T, "sosafldp", "Surface salt flux: Damping"         , "Kg/m2/s",   &  ! erp * sn 
    674                &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    675          ENDIF 
    676           
     667        
    677668         clmx ="l_max(only(x))"    ! max index on a period 
    678669!         CALL histdef( nid_T, "sobowlin", "Bowl Index"                         , "W-point",   &  ! bowl INDEX  
     
    814805      ENDIF 
    815806 
    816       IF( .NOT. ln_cpl ) THEN 
     807      IF( ln_ssr ) THEN 
    817808         CALL histwrite( nid_T, "sohefldp", it, qrp           , ndim_hT, ndex_hT )   ! heat flux damping 
    818809         CALL histwrite( nid_T, "sowafldp", it, erp           , ndim_hT, ndex_hT )   ! freshwater flux damping 
    819          IF( ln_ssr ) zw2d(:,:) = erp(:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1) 
    820          CALL histwrite( nid_T, "sosafldp", it, zw2d          , ndim_hT, ndex_hT )   ! salt flux damping 
    821       ENDIF 
    822       IF( ln_cpl .AND. nn_ice <= 1 ) THEN 
    823          CALL histwrite( nid_T, "sohefldp", it, qrp           , ndim_hT, ndex_hT )   ! heat flux damping 
    824          CALL histwrite( nid_T, "sowafldp", it, erp           , ndim_hT, ndex_hT )   ! freshwater flux damping 
    825          IF( ln_ssr ) zw2d(:,:) = erp(:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1) 
     810         zw2d(:,:) = erp(:,:) * tsn(:,:,1,jp_sal) * tmask(:,:,1) 
    826811         CALL histwrite( nid_T, "sosafldp", it, zw2d          , ndim_hT, ndex_hT )   ! salt flux damping 
    827812      ENDIF 
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/TRA/eosbn2.F90

    r11317 r11325  
    3030   !!   eos_insitu_2d : Compute the in situ density for 2d fields 
    3131   !!   bn2           : Compute the Brunt-Vaisala frequency 
     32   !!   bn2           : compute the Brunt-Vaisala frequency 
     33   !!   eos_pt_from_ct: compute the potential temperature from the Conservative Temperature 
    3234   !!   eos_rab       : generic interface of in situ thermal/haline expansion ratio  
    3335   !!   eos_rab_3d    : compute in situ thermal/haline expansion ratio 
     
    7476 
    7577   !                               !!** Namelist nameos ** 
    76    LOGICAL , PUBLIC ::   ln_TEOS10   ! determine if eos_pt_from_ct is used to compute sst_m 
    77    LOGICAL , PUBLIC ::   ln_EOS80   ! determine if eos_pt_from_ct is used to compute sst_m 
    78    LOGICAL , PUBLIC ::   ln_SEOS   ! determine if eos_pt_from_ct is used to compute sst_m 
     78   LOGICAL , PUBLIC ::   ln_TEOS10 
     79   LOGICAL , PUBLIC ::   ln_EOS80 
     80   LOGICAL , PUBLIC ::   ln_SEOS 
    7981 
    8082   ! Parameters 
     
    16471649         ! 
    16481650      CASE( np_seos )                        !==  Simplified EOS     ==! 
     1651 
     1652         r1_S0  = 0.875_wp/35.16504_wp   ! Used to convert CT in potential temperature when using bulk formulae (eos_pt_from_ct) 
     1653          
    16491654         IF(lwp) THEN 
    16501655            WRITE(numout,*) 
  • utils/build/arch/CNRS/arch-X64_IRENE_DEBUG.fcm

    r10339 r11325  
    4444%CPP               cpp 
    4545%FC                  mpif90 -c -cpp 
    46 %FCFLAGS             -i4 -r8 -g -O0 -debug all -traceback -fp-model precise -ftrapuv -fpe0 
     46%FCFLAGS             -i4 -r8 -g -O0 -debug all -traceback -fp-model precise -ftrapuv -fpe0 -check bounds 
    4747%FFLAGS              %FCFLAGS 
    4848%LD                  mpif90 
  • utils/build/arch/arch-osx_gfortran_debug.fcm

    r10533 r11325  
    4242%NCDF_HOME           /opt/local 
    4343%HDF5_HOME           /opt/local 
    44 %XIOS_HOME           /Users/$( whoami )/xios-2.0 
     44%XIOS_HOME           /Users/$( whoami )/xios-2.5 
    4545%OASIS_HOME          /not/defined 
    4646 
     
    5454%CPP               cpp-mp-4.8 -Dkey_nosignedzero  
    5555%FC                mpif90  
    56 #%FCFLAGS             -fdefault-real-8 -O3 -funroll-all-loops -fcray-pointer -ffree-line-length-none -fbounds-check 
    57 %FCFLAGS             -fdefault-real-8 -O0 -funroll-all-loops -fcray-pointer -ffree-line-length-none -fcheck=bounds -finit-real=nan 
     56%FCFLAGS             -fdefault-real-8 -O0 -g -fbacktrace -funroll-all-loops -fcray-pointer -ffree-line-length-none -fcheck=bounds -finit-real=nan 
    5857%FFLAGS              %FCFLAGS 
    5958%LD                  %FC 
Note: See TracChangeset for help on using the changeset viewer.