Changeset 5748 for LMDZ6/trunk/libf/dyn3d
- Timestamp:
- Jul 2, 2025, 12:00:08 PM (5 weeks ago)
- Location:
- LMDZ6/trunk/libf/dyn3d
- Files:
-
- 3 edited
-
advtrac.f90 (modified) (1 diff)
-
check_isotopes.f90 (modified) (5 diffs)
-
dynetat0.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/advtrac.f90
r5324 r5748 15 15 USE control_mod, ONLY: iapp_tracvl, day_step 16 16 USE comconst_mod, ONLY: dtvr 17 USE strings_mod, ONLY: int2str18 17 USE dimensions_mod, ONLY: iim, jjm, llm, ndm 19 18 USE paramet_mod_h -
LMDZ6/trunk/libf/dyn3d/check_isotopes.f90
r5271 r5748 1 1 SUBROUTINE check_isotopes_seq(q, ip1jmp1, err_msg) 2 USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, int2str, real2str2 USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, num2str 3 3 USE infotrac, ONLY: nqtot, niso, nphas, isotope, isoCheck, iqIsoPha, isoSelect, & 4 4 ntiso, iH2O, nzone, tracers, isoName, itZonIso … … 80 80 CYCLE 81 81 END IF 82 CALL msg('ixt, iq = '//TRIM(strStack( int2str([ixt,iq]))), modname)83 msg1 = '('//TRIM(strStack( int2str([i,k])))//')'84 CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM( real2str(q1)), modname)85 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM( real2str(q2)), modname)82 CALL msg('ixt, iq = '//TRIM(strStack(num2str([ixt,iq]))), modname) 83 msg1 = '('//TRIM(strStack(num2str([i,k])))//')' 84 CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname) 85 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname) 86 86 CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1) 87 87 END DO … … 111 111 deltaD = (q2/q1/tnat(ixt)-1.)*1000. 112 112 IF(deltaD <= deltaDmax .AND. deltaD >= deltaDmin) CYCLE 113 CALL msg('ixt, iq = '//TRIM(strStack( int2str([ixt,iq]))), modname)114 msg1 = '('//TRIM(strStack( int2str([i,k])))//')'115 CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM( real2str(q1)), modname)116 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM( real2str(q2)), modname)117 CALL msg(TRIM(nm(iiso))//TRIM( real2str(deltaD)), modname)113 CALL msg('ixt, iq = '//TRIM(strStack(num2str([ixt,iq]))), modname) 114 msg1 = '('//TRIM(strStack(num2str([i,k])))//')' 115 CALL msg(TRIM(tracers(iqpar)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname) 116 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname) 117 CALL msg(TRIM(nm(iiso))//TRIM(num2str(deltaD)), modname) 118 118 CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1) 119 119 END DO … … 140 140 deltaD = (q2/q1/tnat(iso_HDO)-1.)*1000. 141 141 IF(deltaD <= deltaDmax .AND. deltaD >= deltaDmin) CYCLE 142 CALL msg('izon, ipha = '//TRIM(strStack( int2str([izon, ipha]))), modname)143 CALL msg( 'ixt, ieau = '//TRIM(strStack( int2str([ ixt, ieau]))), modname)144 msg1 = '('//TRIM(strStack( int2str([i,k])))//')'145 CALL msg(TRIM(tracers(iqeau)%name)//TRIM(msg1)//' = '//TRIM( real2str(q1)), modname)146 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM( real2str(q2)), modname)147 CALL msg('deltaD = '//TRIM( real2str(deltaD)), modname)142 CALL msg('izon, ipha = '//TRIM(strStack(num2str([izon, ipha]))), modname) 143 CALL msg( 'ixt, ieau = '//TRIM(strStack(num2str([ ixt, ieau]))), modname) 144 msg1 = '('//TRIM(strStack(num2str([i,k])))//')' 145 CALL msg(TRIM(tracers(iqeau)%name)//TRIM(msg1)//' = '//TRIM(num2str(q1)), modname) 146 CALL msg(TRIM(tracers(iq )%name)//TRIM(msg1)//' = '//TRIM(num2str(q2)), modname) 147 CALL msg('deltaD = '//TRIM(num2str(deltaD)), modname) 148 148 CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1) 149 149 END DO … … 163 163 IF(ABS(xtractot-xiiso) > errmax .AND. ABS(xtractot-xiiso)/MAX(MAX(ABS(xtractot),ABS(xiiso)),1e-18) > errmaxrel) THEN 164 164 CALL msg('Error in iso_verif_aberrant trac: '//TRIM(err_msg)) 165 CALL msg('iiso, ipha = '//TRIM(strStack( int2str([iiso, ipha]))), modname)166 CALL msg('q('//TRIM(strStack( int2str([i,k])))//',:) = '//TRIM(strStack(real2str(q(i,k,:)))), modname)165 CALL msg('iiso, ipha = '//TRIM(strStack(num2str([iiso, ipha]))), modname) 166 CALL msg('q('//TRIM(strStack(num2str([i,k])))//',:) = '//TRIM(strStack(num2str(q(i,k,:)))), modname) 167 167 CALL abort_gcm(modname, 'Error with isotopes: '//TRIM(err_msg), 1) 168 168 END IF -
LMDZ6/trunk/libf/dyn3d/dynetat0.f90
r5285 r5748 10 10 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, & 11 11 new2oldH2O, newHNO3, oldHNO3 12 USE strings_mod, ONLY: maxlen, msg, strStack, real2str, int2str12 USE strings_mod, ONLY: maxlen, msg, strStack, num2str 13 13 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, NF90_INQ_VARID, & 14 14 NF90_CLOSE, NF90_GET_VAR, NF90_NoErr … … 107 107 108 108 !------------------------------------------------------------------------------- 109 CALL msg('rad, omeg, g, cpp, kappa = '//TRIM(strStack( real2str([rad,omeg,g,cpp,kappa]))), modname)109 CALL msg('rad, omeg, g, cpp, kappa = '//TRIM(strStack(num2str([rad,omeg,g,cpp,kappa]))), modname) 110 110 CALL check_dim(im,iim,'im','im') 111 111 CALL check_dim(jm,jjm,'jm','jm') … … 210 210 CHARACTER(LEN=*), INTENT(IN) :: str1, str2 211 211 CHARACTER(LEN=maxlen) :: s1, s2 212 IF(n1/=n2) CALL abort_gcm(TRIM(modname), 'value of "'//TRIM(str1)//'" = '//TRIM( int2str(n1))// &213 ' read in starting file differs from gcm value of "'//TRIM(str2)//'" = '//TRIM( int2str(n2)), 1)212 IF(n1/=n2) CALL abort_gcm(TRIM(modname), 'value of "'//TRIM(str1)//'" = '//TRIM(num2str(n1))// & 213 ' read in starting file differs from gcm value of "'//TRIM(str2)//'" = '//TRIM(num2str(n2)), 1) 214 214 END SUBROUTINE check_dim 215 215
Note: See TracChangeset
for help on using the changeset viewer.
