Changeset 5953 for LMDZ6/trunk


Ignore:
Timestamp:
Dec 18, 2025, 6:54:06 PM (3 weeks ago)
Author:
yann meurdesoif
Message:

Desactivate call of ecumev6_flux, coare30_flux_cnrm, and coare_cp
need to be refactored to work properly.
YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/cdrag_mod.f90

    r5950 r5953  
    206206  LOGICAL               :: LPWG
    207207  real, dimension(1) :: z0m_1d, z_0h_1d, sqrt_zdu2_1d, zgeop1_rg_1d  ! Pour COARE CNRM
    208                                                                      ! convert scalar to 1D for call
     208  CHARACTER(len=80) :: abort_message
     209  CHARACTER(len=20) :: modname = 'pbl_surface_newfrac'                                                                     ! convert scalar to 1D for call
    209210
    210211
     
    339340         ENDIF
    340341!         call ini_csts
    341 !$gpum nocall
    342           call ecumev6_flux( knon, 1, z_0m,t1(i),tsurf(i),&
    343              q1(i),qsurf(i),sqrt(zdu2),zgeop1(i)/RG,PSSS,zgeop1(i)/RG,&
    344              psol(i),pat1(i), OPRECIP, OPWEBB,&
    345              PSFTH,PFSTQ,PUSTAR,PCD,PCDN,&
    346              PCH,PCE,PRI,PRESA,prain,&
    347              z_0h,OPERTFLUX,coeffs)
    348 
     342!!$gpum nocall
     343!ym : bricolage à la petite semaine, ne peut pas marcher sauf en 1D, et encore, incorrecte syntaxiquement (SCALAR<->ARRAY)
     344!ym desactivated
     345!          call ecumev6_flux( i, knon, 1, z_0m,t1(i),tsurf(i),&
     346!             q1(i),qsurf(i),sqrt(zdu2),zgeop1(i)/RG,PSSS,zgeop1(i)/RG,&
     347!             psol(i),pat1(i), OPRECIP, OPWEBB,&
     348!             PSFTH,PFSTQ,PUSTAR,PCD,PCDN,&
     349!             PCH,PCE,PRI,PRESA,prain,&
     350!             z_0h,OPERTFLUX,coeffs)
     351         abort_message="cannot call ecumev6_flux, routine mut be refactored"
     352         CALL abort_physic(modname,abort_message,1)
    349353         cdmm(i) = coeffs(1)
    350354         cdhh(i) = coeffs(2)
     
    365369           zgeop1_rg_1d=zgeop1(i)/RG
    366370
    367 !$gpum nocall
    368            call coare30_flux_cnrm(1, z0m_1d,t1(i),tsurf(i), q1(i),  &
    369                sqrt_zdu2_1d,zgeop1_rg_1d,zgeop1_rg_1d,psol(i),qsurf(i),PQSAT, &
    370                PSFTH,PFSTQ,PUSTAR,PCD,PCDN,PCH,PCE,PRI, &
    371                PRESA,prain,pat1(i),z_0h_1d, LPRECIP, LPWG, coeffs)
    372 
     371!!$gpum nocall
     372!ym : bricolage à la petite semaine, ne peut pas marcher sauf en 1D, et encore, incorrecte syntaxiquement (SCALAR<->ARRAY)
     373!ym desactivated
     374!           call coare30_flux_cnrm(1, z0m_1d,t1(i),tsurf(i), q1(i),  &
     375!               sqrt_zdu2_1d,zgeop1_rg_1d,zgeop1_rg_1d,psol(i),qsurf(i),PQSAT, &
     376!               PSFTH,PFSTQ,PUSTAR,PCD,PCDN,PCH,PCE,PRI, &
     377!               PRESA,prain,pat1(i),z_0h_1d, LPRECIP, LPWG, coeffs)
     378         abort_message="cannot call coare30_flux_cnrm, routine mut be refactored"
     379         CALL abort_physic(modname,abort_message,1)
    373380         cdmm(i) = coeffs(1)
    374381         cdhh(i) = coeffs(2)
     
    384391         write(*,*) "debug size",size(coeffs)
    385392
    386 !$gpum nocall
    387          call coare_cp(sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),&
    388                t1(i),q1(i),&
    389                zgeop1(i)/RG,zgeop1(i)/RG,zgeop1(i)/RG,&
    390                psol(i), pblh(i),&
    391                nit_bulk,&
    392                coeffs,z_0m,z_0h)
     393!!$gpum nocall
     394!ym : bricolage à la petite semaine, ne peut pas marcher sauf en 1D, et encore, incorrecte syntaxiquement (SCALAR<->ARRAY)
     395!ym desactivated
     396!ym         call coare_cp(sqrt(zdu2),t1(i)-tsurf(i),q1(i)-qsurf(i),&
     397!ym               t1(i),q1(i),&
     398!ym               zgeop1(i)/RG,zgeop1(i)/RG,zgeop1(i)/RG,&
     399!ym               psol(i), pblh(i),&
     400!ym               nit_bulk,&
     401!ym               coeffs,z_0m,z_0h)
     402         abort_message="cannot call coare_cp, routine mut be refactored"
     403         CALL abort_physic(modname,abort_message,1)
     404
    393405         cdmm(i) = coeffs(1)
    394406         cdhh(i) = coeffs(3)
Note: See TracChangeset for help on using the changeset viewer.