- Timestamp:
- Jul 24, 2024, 4:39:59 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writedynav.F90
r5117 r5118 9 9 USE temps_mod, ONLY: itau_dyn 10 10 USE lmdz_description, ONLY: descript 11 USE lmdz_iniprint, ONLY: lunout, prt_level 11 12 12 13 IMPLICIT NONE … … 33 34 include "paramet.h" 34 35 include "comgeom.h" 35 include "iniprint.h"36 36 37 37 ! Arguments 38 38 39 REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) 40 REAL teta(ip1jmp1 *llm), phi(ip1jmp1, llm), ppk(ip1jmp1*llm)41 REAL ps(ip1jmp1), masse(ip1jmp1, llm) 42 REAL phis(ip1jmp1) 39 REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) 40 REAL teta(ip1jmp1 * llm), phi(ip1jmp1, llm), ppk(ip1jmp1 * llm) 41 REAL ps(ip1jmp1), masse(ip1jmp1, llm) 42 REAL phis(ip1jmp1) 43 43 REAL q(ip1jmp1, llm, nqtot) 44 44 INTEGER time … … 47 47 ! Variables locales 48 48 49 INTEGER ndex2d(ip1jmp1), ndexu(ip1jmp1 *llm), ndexv(ip1jm*llm)49 INTEGER ndex2d(ip1jmp1), ndexu(ip1jmp1 * llm), ndexv(ip1jm * llm) 50 50 INTEGER iq, ii, ll 51 REAL tm(ip1jmp1 *llm)51 REAL tm(ip1jmp1 * llm) 52 52 REAL vnat(ip1jm, llm), unat(ip1jmp1, llm) 53 53 LOGICAL ok_sync … … 74 74 ! Vents U 75 75 76 CALL histwrite(histuaveid, 'u', itau_w, unat, 77 iip1*jjp1*llm, ndexu)76 CALL histwrite(histuaveid, 'u', itau_w, unat, & 77 iip1 * jjp1 * llm, ndexu) 78 78 79 79 ! Vents V 80 80 81 CALL histwrite(histvaveid, 'v', itau_w, vnat, 82 iip1*jjm*llm, ndexv)81 CALL histwrite(histvaveid, 'v', itau_w, vnat, & 82 iip1 * jjm * llm, ndexv) 83 83 84 84 ! Temperature potentielle moyennee 85 85 86 CALL histwrite(histaveid, 'theta', itau_w, teta, 87 iip1*jjp1*llm, ndexu)86 CALL histwrite(histaveid, 'theta', itau_w, teta, & 87 iip1 * jjp1 * llm, ndexu) 88 88 89 89 ! Temperature moyennee 90 90 91 91 do ii = 1, ijp1llm 92 tm(ii) = teta(ii) * ppk(ii)/cpp92 tm(ii) = teta(ii) * ppk(ii) / cpp 93 93 enddo 94 CALL histwrite(histaveid, 'temp', itau_w, tm, 95 iip1*jjp1*llm, ndexu)94 CALL histwrite(histaveid, 'temp', itau_w, tm, & 95 iip1 * jjp1 * llm, ndexu) 96 96 97 97 ! Geopotentiel 98 98 99 CALL histwrite(histaveid, 'phi', itau_w, phi, 100 iip1*jjp1*llm, ndexu)99 CALL histwrite(histaveid, 'phi', itau_w, phi, & 100 iip1 * jjp1 * llm, ndexu) 101 101 102 102 ! Traceurs … … 109 109 ! Masse 110 110 111 CALL histwrite(histaveid, 'masse', itau_w, masse, 112 iip1*jjp1*llm, ndexu)111 CALL histwrite(histaveid, 'masse', itau_w, masse, & 112 iip1 * jjp1 * llm, ndexu) 113 113 114 114 ! Pression au sol 115 115 116 CALL histwrite(histaveid, 'ps', itau_w, ps, iip1 *jjp1, ndex2d)116 CALL histwrite(histaveid, 'ps', itau_w, ps, iip1 * jjp1, ndex2d) 117 117 118 118 ! Geopotentiel au sol … … 121 121 122 122 IF (ok_sync) THEN 123 124 125 123 CALL histsync(histaveid) 124 CALL histsync(histvaveid) 125 CALL histsync(histuaveid) 126 126 ENDIF 127 127
Note: See TracChangeset
for help on using the changeset viewer.