Dear all Eventually the same kind of problem appears with standard Type 33.a. *TRACE* UNIT 3 TYPE 33 AT TIME 1.3750000000000000E+00 ITERATION 1 CALL 25 POST-CONVERGENCE CALL P 6.0000000000000000E+00 1.0000000000000000E+00 2.0000000000000000E+00 I 8.5248000000000025E-03 5.1599257000000001E+01 1.0000000000000000E+00 O 8.5248000000000025E-03 8.5248000000000025E-03 5.1599257000000001E+01 1.2856158812905980E+00 1.2747489018520892E+00 O 1.0000000000000000E+02 8.5248000000000025E-03 8.5248000000000025E-03 0.0000000000000000E+00 1.0000000000000000E+00 T = O[7] = 8.5E-3 [°C] x = O[1] = 8.5E-3 [] h = O[3] = 51.6 [kJ/kg] (No warning printed in lst file.) I am not sure if it is related with the issue of PSYCHROMETICS routine below. Thank you for your opinion. Regards -- Antoine Gautier antoine.gautier@solamen.fr De : Antoine Gautier [mailto:antoine.gautier@solamen.fr] Dear all Same problem illustrated in a more simplified way. source code MODE=6 PSYDAT(6)=0.8524754597715760D-02 PSYDAT(7)=0.5159925701758463D+02 CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE, & PSYDAT,EMODE,STATUS) write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D, & " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7) CALL MESSAGES(-1,trim(myMessage), & 'Notice',info(1),info(2)) simulation response *** Notice at time : 9.125000 Generated by Unit : 2 Generated by Type : 206 Message : PSYDAT(2)= 0.1732587158021238D+02 PSYDAT(6)= 0.8524754597715760D-02 PSYDAT(7)= 0.5159925701758463D+02 DB Temperature = 17.3°C instead of ~ 30°C for that state point. Is there any explanation? Thanks a lot. Regards -- Antoine Gautier antoine.gautier@solamen.fr De : Antoine Gautier [mailto:antoine.gautier@solamen.fr] Dear all Please have a look below at the short excerpt of a source code. MODE=6 PSYDAT(6)=wAirEnt PSYDAT(7)=hAirEnt CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE, & PSYDAT,EMODE,STATUS) write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D, & " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7) CALL MESSAGES(-1,trim(myMessage), & 'Notice',info(1),info(2)) PSYDAT(2)=29.82D0 CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE, & PSYDAT,EMODE,STATUS) write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D, & " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7) CALL MESSAGES(-1,trim(myMessage), & 'Notice',info(1),info(2)) I know PSYDAT(2) is not used in MODE 6 but please have a look now at the response I get… *** Notice at time : 16.375000 Generated by Unit : 2 Generated by Type : 206 Message : PSYDAT(2)= 0.1732587158021238D+02 PSYDAT(6)= 0.8524754597715760D-02 PSYDAT(7)= 0.5159925701758463D+02 *** Notice at time : 16.375000 Generated by Unit : 2 Generated by Type : 206 Message : PSYDAT(2)= 0.2981999969482422D+02 PSYDAT(6)= 0.8524754597715760D-02 PSYDAT(7)= 0.5159925701758463D+02 Same humidity ratio, same enthalpy and quite different dry bulb temperature. Why the second PSYCHROMETRICS statement does not rewrite on PSYDAT(2)? Any help would be quite appreciated as I am slowly but surely getting crazy! Thanks a lot. Regards -- Antoine Gautier antoine.gautier@solamen.fr |