Ignore:
Timestamp:
Mar 20, 2024, 4:29:23 PM (9 months ago)
Author:
jleconte
Message:

Correction on moistadj_generic tendencies (Noe C)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r3277 r3278  
    16901690               pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)     &
    16911691                                                  + dqmoist(1:ngrid,1:nlayer,1:nq)
    1692                pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq)     &
    1693                                                   + dqmoist(1:ngrid,1:nlayer,1:nq)
    16941692               pdt(1:ngrid,1:nlayer) = pdt(1:ngrid,1:nlayer)+dtmoist(1:ngrid,1:nlayer)
    16951693
     
    17061704                 
    17071705                  if (is_master) then
    1708                      print*,'In moistadj atmospheric energy change   =',dEtot,' W m-2'
    1709                      print*,'In moistadj MAX atmospheric energy change   =',dtmoist_max*ptimestep,'K/step'
    1710                      print*,'In moistadj MIN atmospheric energy change   =',dtmoist_min*ptimestep,'K/step'
     1706                     print*,'In moistadj_generic atmospheric energy change   =',dEtot,' W m-2'
     1707                     print*,'In moistadj_generic MAX atmospheric energy change   =',dtmoist_max*ptimestep,'K/step'
     1708                     print*,'In moistadj_generic MIN atmospheric energy change   =',dtmoist_min*ptimestep,'K/step'
    17111709                  endif
    17121710
    17131711                  ! igcm_generic_vap & igcm_generic_ice are not declared below
    1714                   ! call planetwide_sumval(massarea(:,:)*dqmoist(:,:,igcm_generic_vap)*ptimestep/totarea_planet+        &
    1715                   !                        massarea(:,:)*dqmoist(:,:,igcm_generic_ice)*ptimestep/totarea_planet,dWtot)
    1716                   ! if (is_master) print*,'In moistadj atmospheric water change    =',dWtot,' kg m-2'
     1712                  call planetwide_sumval(massarea(:,:)*dqmoist(:,:,1)*ptimestep/totarea_planet+        &
     1713                                         massarea(:,:)*dqmoist(:,:,2)*ptimestep/totarea_planet,dWtot)
     1714                  if (is_master) print*,'In moistadj_generic atmospheric GCS change    =',dWtot,' kg m-2'
    17171715               
    17181716               endif ! end of 'enertest'
     
    18091807                     call planetwide_sumval(cell_area(:)*zdqssnow_generic(:,igcm_generic_ice)/totarea_planet*RLVTT/cpp,dItot)
    18101808                     dItot = dItot + dItot_tmp
    1811                      call planetwide_sumval(massarea(:,:)*dq_rain_generic_vap(:,:,igcm_generic_vap)*ptimestep/totarea_planet,dVtot_tmp)
     1809                     call planetwide_sumval(massarea(:,:)*dq_rain_generic_vap(:,:,igcm_generic_vap)/totarea_planet*RLVTT/cpp,dVtot_tmp)
    18121810                     call planetwide_sumval(cell_area(:)*zdqsrain_generic(:,igcm_generic_ice)/totarea_planet*RLVTT/cpp,dVtot)
    18131811                     dVtot = dVtot + dVtot_tmp
Note: See TracChangeset for help on using the changeset viewer.