Changeset 3558
- Timestamp:
- Dec 17, 2024, 2:17:56 PM (5 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
r3557 r3558 1567 1567 1568 1568 pdq(:,:,:) = pdq(:,:,:) + pdqmufi(:,:,:) 1569 1570 ! [BBT: Temporal tests]1571 ! >>>>>>>>>>>>>>>>>>>>>1572 !write(*,*) 'pdqmufi_prod', MAXVAL(pdqmufi_prod(1,:,11)), 'i', MAXLOC(pdqmufi_prod(1,:,11))1573 open(187,file='Mufi_tracers.out')1574 write(187,*) "Pressure (Pa) - Alt (m) - M0as (X/m3) - M3as (m3/m3) - M0af (X/m3) - M3af (m3/m3)"1575 do l=1,nlayer1576 write(187,*) pplay(1,l),zzlay(1,l), &1577 (pq(1,l,micro_indx(1))+pdq(1,l,micro_indx(1))*ptimestep)*int2ext(1,l), &1578 (pq(1,l,micro_indx(2))+pdq(1,l,micro_indx(2))*ptimestep)*int2ext(1,l), &1579 (pq(1,l,micro_indx(3))+pdq(1,l,micro_indx(3))*ptimestep)*int2ext(1,l), &1580 (pq(1,l,micro_indx(4))+pdq(1,l,micro_indx(4))*ptimestep)*int2ext(1,l)1581 enddo1582 close(187)1583 open(188,file='Mufi_tendencies.out')1584 write(188,*) "Pressure (Pa) - Alt (m) - dM0as (X/kg_air/s) - dM3as (m3/kg_air/s) - dM0af (X/kg_air/s) - dM3af (m3/kg_air/s)"1585 do l=1,nlayer1586 write(188,*) pplay(1,l),zzlay(1,l),pdqmufi(1,l,micro_indx(1)),pdqmufi(1,l,micro_indx(2)),pdqmufi(1,l,micro_indx(3)),pdqmufi(1,l,micro_indx(4))1587 enddo1588 close(188)1589 open(189,file='Mufi_diags.out')1590 write(189,*) "Pressure (Pa) - Alt (m) - rc_sph (m) - rc_fra (m) - aer_w_sph (m/s) - aer_w_fra (m/s) - aer_prec_sph (kg/m2/s) - aer_prec_fra (kg/m2/s)"1591 do l=1,nlayer1592 write(189,*) pplay(1,l),zzlay(1,l),mp2m_rc_sph(1,l),mp2m_rc_fra(1,l),mp2m_aer_s_w(1,l),mp2m_aer_f_w(1,l),mp2m_aer_s_prec(1),mp2m_aer_f_prec(1)1593 enddo1594 close(189)1595 ! <<<<<<<<<<<<<<<<<<<<<1596 1569 ENDIF ! end callmufi 1597 1570
Note: See TracChangeset
for help on using the changeset viewer.