Ignore:
Timestamp:
Jun 20, 2014, 12:27:40 PM (10 years ago)
Author:
Laurent Fairhead
Message:

Modifications pour avoir 1+1=2 avec la nouvelle physique (testé avec NPv4.12)


Modifications needed to ensure 1+1=2 with the new physics (tested with NPv4.12)

Location:
LMDZ5/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phyetat0.F90

    r2057 r2069  
    1515       solsw, t_ancien, u_ancien, v_ancien, w01, wake_cstar, wake_deltaq, &
    1616       wake_deltat, wake_fip, wake_pe, wake_s, zgam, zmax0, zmea, zpic, zsig, &
    17        zstd, zthe, zval
     17       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
    1818  USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
    1919  USE infotrac, only: nbtr, type_trac, tname, niadv
     
    997997  PRINT*, '(ecart-type) detr_therm:', xmin, xmax
    998998
     999  CALL get_field("ALE_BL", ale_bl, found)
     1000  IF (.NOT. found) THEN
     1001     PRINT*, "phyetat0: Le champ <ale_bl> est absent"
     1002     PRINT*, "Depart legerement fausse. Mais je continue"
     1003     ale_bl=0.
     1004  ENDIF
     1005  xmin = 1.0E+20
     1006  xmax = -1.0E+20
     1007  xmin = MINval(ale_bl)
     1008  xmax = MAXval(ale_bl)
     1009  PRINT*, '(ecart-type) ale_bl:', xmin, xmax
     1010
     1011  CALL get_field("ALE_BL_TRIG", ale_bl_trig, found)
     1012  IF (.NOT. found) THEN
     1013     PRINT*, "phyetat0: Le champ <ale_bl_trig> est absent"
     1014     PRINT*, "Depart legerement fausse. Mais je continue"
     1015     ale_bl_trig=0.
     1016  ENDIF
     1017  xmin = 1.0E+20
     1018  xmax = -1.0E+20
     1019  xmin = MINval(ale_bl_trig)
     1020  xmax = MAXval(ale_bl_trig)
     1021  PRINT*, '(ecart-type) ale_bl_trig:', xmin, xmax
     1022
     1023  CALL get_field("ALP_BL", alp_bl, found)
     1024  IF (.NOT. found) THEN
     1025     PRINT*, "phyetat0: Le champ <alp_bl> est absent"
     1026     PRINT*, "Depart legerement fausse. Mais je continue"
     1027     alp_bl=0.
     1028  ENDIF
     1029  xmin = 1.0E+20
     1030  xmax = -1.0E+20
     1031  xmin = MINval(alp_bl)
     1032  xmax = MAXval(alp_bl)
     1033  PRINT*, '(ecart-type) alp_bl:', xmin, xmax
     1034
    9991035  ! Read and send field trs to traclmdz
    10001036
  • LMDZ5/trunk/libf/phylmd/phyredem.F90

    r2057 r2069  
    323323
    324324  CALL put_field("DETR_THERM", "DETR_THERM", detr_therm)
     325
     326  CALL put_field("ALE_BL", "ALE_BL", Ale_bl)
     327
     328  CALL put_field("ALE_BL_TRIG", "ALE_BL_TRIG", Ale_bl_trig)
     329
     330  CALL put_field("ALP_BL", "ALP_BL", Alp_bl)
    325331
    326332  ! trs from traclmdz_mod
  • LMDZ5/trunk/libf/phylmd/physiq.F90

    r2068 r2069  
    11211121        endif
    11221122
    1123         do i = 1,klon
    1124            Ale_bl(i)=0.
    1125            Alp_bl(i)=0.
    1126         enddo
     1123!        do i = 1,klon
     1124!           Ale_bl(i)=0.
     1125!           Alp_bl(i)=0.
     1126!        enddo
    11271127
    11281128        !================================================================================
     
    19861986        enddo
    19871987!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     1988
    19881989
    19891990     endif
Note: See TracChangeset for help on using the changeset viewer.